# Experimental and Numerical Study of Texture Evolution and Anisotropic Plastic Deformation of Pure Magnesium under Various Strain Paths

Authors: Hamad F. Alharbi; Monis Luqman; Ehab El-Danaf; Nabeel H. Alharthi
Journal: Advances in Materials Science and Engineering
 (2018)

---

## Abstract

The deformation behavior and texture evolution of pure magnesium were investigated during plane strain compression, simple compression, and uniaxial tension at room temperature. The distinctive stages in the measured anisotropic stress-strain responses and numerically computed strain-hardening rates were correlated with texture and deformation mechanisms. More specifically, in plane strain compression and simple compression, the onset of tensile twins and the accompanying texture-hardening effect were associated with the initial high strain-hardening rates observed in specimens loaded in directions perpendicular to the crystallographicc-axis in most of the grains. The subsequent drop in strain-hardening rates in these samples was correlated with the exhaustion of tensile twins and the activation of pyramidal <c+a> slip systems. The falling strain-hardening rates were observed in simple compression and plane strain compression with loading directions parallel to the c-axis where the second pyramidal <c+a> slip systems were the only slip families that can accommodate deformation. For uniaxial tension with the basal plane parallel to the tensile axis, the prismatic <a> and second pyramidal <c+a> slips are the main deformation mechanisms. The predicted relative slip and twin activities from the crystal plasticity simulations clearly showed the effect of texture on the type of activated deformation mechanisms.


---

## Body

## 1. Introduction

Magnesium and its alloys have potential applications in automobile and aerospace industries due to their low density, high specific strength, and good castability [1, 2]. However, magnesium materials have low ductility at room temperature compared to other ductile metals [3, 4]. The low ductility in magnesium can be generally attributed to the lack of independent easy-activated slip systems that can accommodate large plastic strains.Magnesium has a close-packed (hcp) crystal structure and exhibits a high degree of plastic anisotropy due to the possible activation of various slip and twin families and the resulting pronounced texture. The slip systems include the basal <a>, prism <a>, pyramidal <a>, and second pyramidal <c+a> slip systems [5–7]. The basal <a> family is the easiest slip system to activate in magnesium with lower slip resistance compared to the other systems. The second pyramidal <c+a> slip system is the only slip family that can accommodate plastic strain in the basal normal direction (all the other slip systems operate in the basal <a> direction only), but it is also very hard to activate. Twinning is an alternative deformation mechanism in magnesium that can also accommodate plastic strains in the nonbasal directions. The {101¯2}check space<1¯011> tensile twins and the {101¯1}check space<101¯2¯> contraction twins are the most common twin systems observed in magnesium and its alloys [8–11]. The extension twins give extension along the crystallographic c-axis in magnesium and reorient the crystal by 86° around <12¯10> directions, while the contraction twins supply contraction along the crystallographic c-axis and reorient the crystal by 56° around <12¯10> directions [8, 9].The extension twins are readily observed in magnesium and its alloys during simple compression tests when the loading direction is perpendicular to the crystallographicc-axis of the crystal. Such deformation modes are commonly triggered in textured (polycrystalline) magnesium when compressed parallel to the extrusion direction of an extruded bar or compressed perpendicular to the normal direction of a rolled sheet [10–12]. The extension twins can also be activated in magnesium during a uniaxial tension if the tensile loading direction is parallel to the crystallographic c-axis of the lattice crystals [13, 14]. The activation of extension twins in magnesium has important implications on the texture evolution and strain-hardening rates. The initial high strain-hardening rates observed in magnesium are commonly linked to the activation of extension twins in the early stage of deformation. The widely accepted hypothesis for the physical mechanism by which deformation twinning contributes to the unusually high strain-hardening rates is the texture hardening, where tensile twins usually reorient the grains into hard (unfavored slip) orientations. An alternative hypothesis for explaining the high strain-hardening rates associated with the onset of deformation twinning is the Hall–Petch-like effect, where twinning divides the grains thereby reducing the effective slip distance [15–18]. However, the latter explanation is not likely to be effective due to the fast growth (thickening) of tensile twins in magnesium until it consumes the entire grain.The current study presents detailed investigations of the mechanisms that govern the anisotropic plastic response of pure magnesium under various strain paths at room temperature. We have measured the anisotropic stress-strain response in different mechanical tests including plane strain compression, simple compression, and uniaxial tensile tests in different sample directions that were initially processed by extrusion and rolling. We have also measured the texture of the initial and deformed samples using electron backscattered diffraction (EBSD) technique. The differences in the strain-hardening responses and texture evolution between uniaxial loading and plane strain compression tests under certain loading directions were discussed. The highly observed plastic anisotropy in pure Mg was related to the type of predominant slip and/or twin deformation modes in each sample orientation and supported by predictions from crystal plasticity simulation. Such information is important for developing more accurate crystal plasticity models for Mg alloys.


## 2. Experimental Procedure

The material used in this study was an as-cast billet of high-purity magnesium (99.9%). A cylindrical bar with 35 mm (diameter) × 50 mm (long) and a plate of 70 mm × 58 mm × 19 mm were cut from the billet for extrusion and rolling processing, respectively. The cylindrical bar and plate were first annealed at 250°C in a vacuum tube furnace for 4 hrs. The cylindrical bar was then hot extruded at 250°C from 35 mm to 15 mm (∼5.4 extrusion ratio), and the plate was hot rolled at about 200°C from 19 mm to 8.6 mm (∼55% reduction in thickness) in multiple passes. The extruded bar and rolled plate were subsequently used to fabricate five sets of samples with different orientations to investigate the effect of strain path change and the associated texture on the anisotropic mechanical response of pure magnesium under various mechanical tests.Five sets of samples with different orientations were machined from the cylindrical extruded bar and rolled plate for subsequent mechanical tests as depicted in Figure1. The orientations of the five set of samples were designed such that the loading directions in the mechanical tests align either perpendicular or parallel to the typical crystallographic orientations of the c-axis resulting from extrusion and rolling processes. These samples were subjected to plane strain compression (PSC) (ND90ED0 and ND0ED0 samples), simple compression (CA0ED0 and CA0ND0 samples), and uniaxial tensile (TA90ND0 sample) tests at room temperature. The orientations of these samples can be distinguished by their labels, where the number in the label illustrates the angle between the loading direction in the mechanical test and the original processing direction (extrusion direction ED0 of the extruded bar or normal direction ND0 of the rolled plate). As an example, the label of ND90ED0 sample indicates that the loading (normal) direction ND in the PSC test makes an angle of 90° with respect to the original extrusion direction ED0. The labels also show the rolling (RD), transverse (TD), and normal (ND) directions in plane strain compression, compression axis (CA) in simple compression, and tensile axis (TA) in uniaxial tensile tests. Two specimens for each set of orientations were mechanically tested, and the variance in the measured stress-strain responses between the two samples tested for each group was found to be very small.Figure 1
Schematic illustration of the orientations of the samples machined from the cylindrical extruded bar and rolled plate and used in plane strain compression (ND90ED0 and ND0ED0), simple compression (CA0ED0 and CA0ND0), and uniaxial tensile (TA90ND0) tests. The typical crystallographic orientations of the c-axis resulting from extrusion and rolling processes are also shown.The plane strain compression (PSC) tests were carried out in a channel die that allows an elongation of the specimens in RD direction only, whereas deformation is suppressed in TD direction. Teflon sheet and high-pressure grease were used between the sample and the die surfaces to reduce frictional effect. The dimensions of the specimens for the PSC tests were 11 mm (RD) × 8 mm (TD) × 7 mm (ND). As the labels of the PSC specimens indicate, the loading direction ND is perpendicular to the original extrusion direction ED0 in ND90ED0 samples, whereas it is parallel to ED0 in ND0ED0 samples. For simple compression test, cylindrical samples were cut with dimensions of 13 mm (diameter) × 20 mm (height) and 5 mm (diameter) × 8 mm (height) for CA0ED0 and CA0ND0, respectively. The CA0ED0 samples were machined from the cylindrical extruded bar with the compression axis CA parallel to ED0, while CA0ND0 samples were cut from the rolled plate with CA parallel to the plate’s normal direction ND0. The tensile specimens TA90ND0 were machined from the rolled plate in the TD0-RD0 plane with the tensile axis TA parallel to the plate’s rolling direction RD0. The tensile samples were machined as per ASTM E8 with a gauge length of 25 mm. All mechanical tests were performed at room temperature with a constant initial strain rate of about 10−3 s−1, using a screw-driven INSTRON machine. The raw data of the PSC and simple compression tests were corrected for machine compliance. The change of gauge length in the uniaxial tensile tests was measured using an axial extensometer (Epsilon 3542-0125M).The texture and microstructure evolutions of all samples were characterized using a field-emission scanning electron microscope (JEOL JSM-7600F) equipped with an electron backscattered diffraction (EBSD) camera supplied by Oxford Instruments. The specimens were prepared for EBSD scanning by mechanical polishing using subsequent oil-based diamond suspensions of 9, 3, and 1µm, respectively. This was followed by 4–6 s chemical etching using a solution of 60 ml ethanol, 20 ml distilled water, 15 ml acetic acid, and 5 ml nitric acid. The EBSD scan areas for all samples were taken as large as possible (>700 µm × 700 µm) to provide a more representative texture information, but with a slightly large step size of 1 µm to be able to complete the scanning in a reasonable time.


## 3. Experimental Results

### 3.1. True Stress-Strain Curves and Strain-Hardening Rates

The true stress-strain curves measured on different pure magnesium samples under various mechanical tests including PSC, simple compression, and uniaxial tensile are shown in Figure2. In order to easily compare the mechanical responses between the different tests, the flow curves are plotted based on the von Mises definition of the equivalent stress and equivalent strain [19]. The strong anisotropic plastic deformation responses shown in Figure 2 have been reported in prior literature for magnesium and its alloys (e.g. [11, 13, 14, 20, 21]). Figure 3 shows the corresponding strain-hardening rates computed numerically from the stress-strain curves. The strain-hardening rates were normalized by the shear modulus of magnesium. This figure clearly indicates the major differences in the strain-hardening behavior between the different specimens tested in this study.Figure 2
Flow curves of pure magnesium at room temperature under different mechanical tests including PSC (ND90ED0 and ND0ED0), simple compression (CA0ED0 and CA0ND0), and uniaxial tensile (TA90ND0). The loading directions with respect to the typical crystallographic c-axis in most grains for each specimen are shown between parentheses.Figure 3
Normalized strain-hardening responses of pure magnesium at room temperature under different mechanical tests including PSC (ND90ED0 and ND0ED0), simple compression (CA0ED0 and CA0ND0), and uniaxial tensile (TA90ND0).


### 3.2. Microstructure and Texture Evolution

Microstructure and texture evolutions of all specimens were investigated by EBSD scanning which reveals both the structure and crystallographic orientations of the grains. Figure4 shows the EBSD inverse pole figure (IPF) maps and pole figures of the cylindrical extruded bar and rolled plate that were used to fabricate all five sets of samples shown in Figure 1. The IPF maps are collected on planes perpendicular to the extrusion direction (ED0) and normal direction (ND0) of the cylindrical extruded bar and rolled plate, respectively. Therefore, the IPF maps reflect the orientations of the extruded bar’s extrusion direction (ED0) and rolled plate’s normal direction (ND0) with respect to the local crystal frames. The measured pole figures clearly indicate that the extruded bar exhibited a typical prismatic 101¯0 fiber texture, where basal planes oriented parallel to the extrusion direction ED0 (or the crystallographic c-axis in most grains oriented in a radial direction that is perpendicular to ED0) (e.g. [22–27]). The rolled plate also showed a typical fiber texture with the crystallographic c-axis in most grains aligned parallel to plate’s normal direction ND0 (e.g. [21, 28, 29]). Based on this texture, the orientations of the five sets of samples were designed such that the loading directions in the mechanical tests align either perpendicular or parallel to the typical crystallographic c-axis in most of the grains. More specifically, the loading direction was selected to be (1) parallel to the crystallographic c-axis in most grains in the samples ND90ED0 and CA0ND0 and (2) perpendicular to the crystallographic c-axis in most grains in the samples ND0ED0, CA0ED0, and TA90ND0. In order to make this clear, we show in Figure 5 the {0001} and {101¯0}check space pole figures of the initial texture for all samples using the new reference frame for each mechanical test which were obtained by appropriate rotations of the pole figures in Figure 4.Figure 4
EBSD inverse pole figure (IPF) maps (top) and pole figures (bottom) showing microstructure and texture of the (a) cylindrical extruded bar and (b) rolled plate of the pure magnesium. The IPF maps reflect the orientations of the bar’s extrusion direction (ED0) and rolled plate’s normal direction (ND0) with respect to the local crystal frames. The measured pole figures show the initial texture of the processed (extruded and rolled) material. The reference frames (X0, Y0, ED0) and (RD0, TD0, ND0) are the same as in Figure 1.
(a)
(b)Figure 5
{
0001
} and {101¯0}check space pole figures of the initial texture for the five set of samples shown in Figure 1 using the reference frame for each mechanical test: (RD, TD, ND) for PSC, (1, 2, CA) for simple compression and (1, 2, TA) for uniaxial tension, where 1 and 2 refer to the lateral directions of the compression or tension test. The pole figures were obtained by appropriate rotations of the pole figures shown in Figure 4 with respect to the new reference frame. (a) ND90ED0. (b) ND0ED0. (c) CA0ED0. (d) CA0ND0. (e) TA90ND0.
(a)
(b)
(c)
(d)
(e)The texture and microstructure of all specimens after each mechanical test are collected by EBSD scanning. For the purpose of simplicity and clarity, the texture will be described in terms of the orientation of the crystallographicc-axis in the crystal with respect to the loading direction in the mechanical test. The texture of ND90ED0 and ND0ED0 samples after PSC tests at room temperature is shown in Figure 6. The EBSD scans were acquired in the RD-TD plane in the midsection of the sample. It is clear that the crystallographic c-axis in most grains aligns parallel to ND direction in both ND90ED0 and ND0ED0 samples after PSC tests (the poles reappear at the center of the pole figures for both samples). It should be noted that even though the starting texture for ND90ED0 and ND0ED0 samples was different (Figures 5(a) and 5(b)), the final texture after PSC tests for the two set of samples was found to be similar. The reason for such texture change is related to the type of deformation mechanisms, which will be explained in more detail in Discussion. Figures 7(a) and 7(b) present the texture after simple compression test at room temperature for CA0ED0 and CA0ND0 samples, respectively. The pole figures are collected in a plane normal to the compression axis CA in the midsection of the sample, with CA lying at the center of the pole figures. It is seen once again that even though the starting texture for CA0ED0 and CA0ND0 samples is different (Figures 5(c) and 5(d)), the final texture after simple compression test is similar to the crystallographic c-axis in most grains that align parallel to the loading direction CA in both CA0ED0 and CA0ND0 samples. Lastly, the texture after uniaxial tensile test for the TA90ND0 sample is illustrated in Figure 7(c) which is collected in the tensile plane about 3 mm apart from the fracture surface in the gauge length region. It is noted that the texture of TA90ND0 sample after the tensile test has little change compared to the initial texture.Figure 6
EBSD inverse pole figure (IPF) maps (left) and pole figures (right) showing microstructure and texture in RD-TD plane of the (a) ND0ED0 and (b) ND90ED0 pure magnesium samples after plane strain compression. The color legend is the same as in Figure 4.
(a)
(b)Figure 7
EBSD inverse pole figure (IPF) maps (left) and pole figures (right) showing texture of (a) CA0ED0 and (b) CA0ND0 pure magnesium samples after simple compression collected in a plane normal to the compression axis CA. The IPF and pole figures in (c) illustrate the texture of r sample after uniaxial tensile test collected in the tensile plane. The labels 1 and 2 in (c) represent the lateral directions of the tension test. The color legend is the same as in Figure 4.
(a)
(b)
(c)


## 3.1. True Stress-Strain Curves and Strain-Hardening Rates

The true stress-strain curves measured on different pure magnesium samples under various mechanical tests including PSC, simple compression, and uniaxial tensile are shown in Figure2. In order to easily compare the mechanical responses between the different tests, the flow curves are plotted based on the von Mises definition of the equivalent stress and equivalent strain [19]. The strong anisotropic plastic deformation responses shown in Figure 2 have been reported in prior literature for magnesium and its alloys (e.g. [11, 13, 14, 20, 21]). Figure 3 shows the corresponding strain-hardening rates computed numerically from the stress-strain curves. The strain-hardening rates were normalized by the shear modulus of magnesium. This figure clearly indicates the major differences in the strain-hardening behavior between the different specimens tested in this study.Figure 2
Flow curves of pure magnesium at room temperature under different mechanical tests including PSC (ND90ED0 and ND0ED0), simple compression (CA0ED0 and CA0ND0), and uniaxial tensile (TA90ND0). The loading directions with respect to the typical crystallographic c-axis in most grains for each specimen are shown between parentheses.Figure 3
Normalized strain-hardening responses of pure magnesium at room temperature under different mechanical tests including PSC (ND90ED0 and ND0ED0), simple compression (CA0ED0 and CA0ND0), and uniaxial tensile (TA90ND0).


## 3.2. Microstructure and Texture Evolution

Microstructure and texture evolutions of all specimens were investigated by EBSD scanning which reveals both the structure and crystallographic orientations of the grains. Figure4 shows the EBSD inverse pole figure (IPF) maps and pole figures of the cylindrical extruded bar and rolled plate that were used to fabricate all five sets of samples shown in Figure 1. The IPF maps are collected on planes perpendicular to the extrusion direction (ED0) and normal direction (ND0) of the cylindrical extruded bar and rolled plate, respectively. Therefore, the IPF maps reflect the orientations of the extruded bar’s extrusion direction (ED0) and rolled plate’s normal direction (ND0) with respect to the local crystal frames. The measured pole figures clearly indicate that the extruded bar exhibited a typical prismatic 101¯0 fiber texture, where basal planes oriented parallel to the extrusion direction ED0 (or the crystallographic c-axis in most grains oriented in a radial direction that is perpendicular to ED0) (e.g. [22–27]). The rolled plate also showed a typical fiber texture with the crystallographic c-axis in most grains aligned parallel to plate’s normal direction ND0 (e.g. [21, 28, 29]). Based on this texture, the orientations of the five sets of samples were designed such that the loading directions in the mechanical tests align either perpendicular or parallel to the typical crystallographic c-axis in most of the grains. More specifically, the loading direction was selected to be (1) parallel to the crystallographic c-axis in most grains in the samples ND90ED0 and CA0ND0 and (2) perpendicular to the crystallographic c-axis in most grains in the samples ND0ED0, CA0ED0, and TA90ND0. In order to make this clear, we show in Figure 5 the {0001} and {101¯0}check space pole figures of the initial texture for all samples using the new reference frame for each mechanical test which were obtained by appropriate rotations of the pole figures in Figure 4.Figure 4
EBSD inverse pole figure (IPF) maps (top) and pole figures (bottom) showing microstructure and texture of the (a) cylindrical extruded bar and (b) rolled plate of the pure magnesium. The IPF maps reflect the orientations of the bar’s extrusion direction (ED0) and rolled plate’s normal direction (ND0) with respect to the local crystal frames. The measured pole figures show the initial texture of the processed (extruded and rolled) material. The reference frames (X0, Y0, ED0) and (RD0, TD0, ND0) are the same as in Figure 1.
(a)
(b)Figure 5
{
0001
} and {101¯0}check space pole figures of the initial texture for the five set of samples shown in Figure 1 using the reference frame for each mechanical test: (RD, TD, ND) for PSC, (1, 2, CA) for simple compression and (1, 2, TA) for uniaxial tension, where 1 and 2 refer to the lateral directions of the compression or tension test. The pole figures were obtained by appropriate rotations of the pole figures shown in Figure 4 with respect to the new reference frame. (a) ND90ED0. (b) ND0ED0. (c) CA0ED0. (d) CA0ND0. (e) TA90ND0.
(a)
(b)
(c)
(d)
(e)The texture and microstructure of all specimens after each mechanical test are collected by EBSD scanning. For the purpose of simplicity and clarity, the texture will be described in terms of the orientation of the crystallographicc-axis in the crystal with respect to the loading direction in the mechanical test. The texture of ND90ED0 and ND0ED0 samples after PSC tests at room temperature is shown in Figure 6. The EBSD scans were acquired in the RD-TD plane in the midsection of the sample. It is clear that the crystallographic c-axis in most grains aligns parallel to ND direction in both ND90ED0 and ND0ED0 samples after PSC tests (the poles reappear at the center of the pole figures for both samples). It should be noted that even though the starting texture for ND90ED0 and ND0ED0 samples was different (Figures 5(a) and 5(b)), the final texture after PSC tests for the two set of samples was found to be similar. The reason for such texture change is related to the type of deformation mechanisms, which will be explained in more detail in Discussion. Figures 7(a) and 7(b) present the texture after simple compression test at room temperature for CA0ED0 and CA0ND0 samples, respectively. The pole figures are collected in a plane normal to the compression axis CA in the midsection of the sample, with CA lying at the center of the pole figures. It is seen once again that even though the starting texture for CA0ED0 and CA0ND0 samples is different (Figures 5(c) and 5(d)), the final texture after simple compression test is similar to the crystallographic c-axis in most grains that align parallel to the loading direction CA in both CA0ED0 and CA0ND0 samples. Lastly, the texture after uniaxial tensile test for the TA90ND0 sample is illustrated in Figure 7(c) which is collected in the tensile plane about 3 mm apart from the fracture surface in the gauge length region. It is noted that the texture of TA90ND0 sample after the tensile test has little change compared to the initial texture.Figure 6
EBSD inverse pole figure (IPF) maps (left) and pole figures (right) showing microstructure and texture in RD-TD plane of the (a) ND0ED0 and (b) ND90ED0 pure magnesium samples after plane strain compression. The color legend is the same as in Figure 4.
(a)
(b)Figure 7
EBSD inverse pole figure (IPF) maps (left) and pole figures (right) showing texture of (a) CA0ED0 and (b) CA0ND0 pure magnesium samples after simple compression collected in a plane normal to the compression axis CA. The IPF and pole figures in (c) illustrate the texture of r sample after uniaxial tensile test collected in the tensile plane. The labels 1 and 2 in (c) represent the lateral directions of the tension test. The color legend is the same as in Figure 4.
(a)
(b)
(c)


## 4. Discussions

### 4.1. Deformation Mechanisms and Texture Evolution

The deformation behavior of magnesium and its alloys can be very different depending on the initial texture and the type of plastic deformation mechanisms. The mechanical deformation behavior of the pure magnesium tested in this study as described by the flow curves shown in Figure2 can be divided into two groups each having a special characteristic in common. The first group (ND0ED0 and CA0ED0 samples) has a concave stress-strain curve consisting of an initial drop in hardening rate followed by a marked increase up to a true strain of 0.1 where another drop is noticed. On the other hand, the second group (ND90ED0, CA0ND0, and TA90ND0 samples) has the typical falling strain-hardening rates similar to the saturation-type hardening observed in high and medium staking fault energy face-centered cubic metals. The offset yield point for the samples in the first group is about 75 MPa, whereas it is about 130 MPa for the samples in the second group. All flow curves in both groups have a sharp loss of strength at the end of the mechanical tests accompanied by the occurrence of inhomogeneous deformation in terms of shear cracks. It should be clear that this classification does not assume the same deformation mechanisms for all samples within the same group. The deformation mechanisms and texture evolution for the samples in each group are discussed next.For the first group consisting of ND0ED0 and CA0ED0 samples, the loading (compression) direction in the mechanical tests (PSC for ND0ED0 and simple compression for CA0ED0) aligns nearly perpendicular to the crystallographic c-axis in most of the grains as the texture indicated in Figures 5(b) and 5(c). This type of deformation mode is known to activate the {101¯2}check space<1¯011> tensile twins extensively at the early stage of deformation [21, 30]. This twin mode reorients the crystal by 86° around <12¯10> directions bringing the basal plane to be nearly perpendicular to the loading direction. The formation of this basal texture ({0001}-fiber texture) leads to texture hardening where most grains become hard to deform with the c-axis nearly parallel to the loading axis. The {101¯2}check space tensile twinning in our samples appears to start at an axial true strain of about 0.03 in ND0ED0 and CA0ED0 samples. This strain value corresponds to the points at which the strain-hardening rates start to increase as illustrated in Figures 2 and 3 (marked using arrows, note that Figure 2 shows the equivalent strain). This observation is in agreement with previously reported experiments in Mg alloys that detected the initial appearance of tensile twins at strains in the range of 0.02 to 0.04 [21, 31]. It should be clear that, in this study, we have not tried to detect the initial appearance of tensile twins experimentally at this small-strain level since this would require interrupting the mechanical tests frequently. The production of extension twins and the accompanied texture hardening lead to the sharp increase in flow stresses only after a small amount of true strain, where the flow stresses increased from about 80 MPa to about 160 MPa by a true strain of about 0.07 in both ND0ED0 and CA0ED0 samples. Therefore, the very high strain-hardening rates observed in samples ND0ED0 and CA0ED0 can be attributed to the activation of extension twins and the associated texture-hardening effect. Furthermore, the exhaustion of extension twins in these samples is estimated to occur at about 0.08 strain. This strain value corresponds to the peak strain-hardening rates shown in Figure 3. It should be noted that this value is very close to the maximum shear strain that can be produced by extension twins, estimated to be 0.1289 shear strain (or about 0.065 axial strain). In other words, most of the grains in ND0ED0 and CA0ED0 samples are expected to form a basal texture due to extension twins at about 0.08 strain. The plastic strain after the exhaustion of extension twins beyond an equivalent strain of about 0.08 is most likely accommodated by the second pyramidal <c+a> slip. Once slip and/or twinning cannot accommodate the shape change, a sudden loss of strength occurred due to the development of an inhomogeneous deformation in terms of shear bands that leads to softening and eventually cracking.It is interesting to note that the flow responses and strain-hardening rates in ND0ED0 and CA0ED0 samples are quite similar in spite of the difference in the applied boundary conditions for the two mechanical tests. To explain this, let us first consider the geometrical constraints in the two mechanical tests. In simple compression with CA normal to the c-axis, the sample has no lateral constraint that would suppress the extension of crystals along the c-axis as a result of extension twins. However, in plane strain compression with ND normal to the c-axis, the crystals with c-axis parallel to TD are not expected to twin because twinning would cause broadening in the constrained direction, TD. Thus, only the crystals with c-axis not perfectly aligned with TD can twin (especially those with c-axis parallel to the free direction, RD). In other words, the amount of twin activities is expected to be higher in simple compression compared to plane strain compression when the loading direction aligns perpendicular to the c-axis. This observation has been numerically investigated using crystal plasticity models for ideal orientations in prior literature (e.g. [20, 30, 32]). However, in the current study, it is noted that the flow responses and strain-hardening rates of the plane strain compressed sample (ND0ED0) are quite similar to the simple compressed sample (CA0ED0). This can be attributed to the existence of many grains in ND0ED0 sample with enough scatter of the c-axis about TD and therefore the easy activation of extension twins compared to other slip systems. This is clearly evidenced by the fact that the final texture of ND0ED0 sample at the end of the mechanical test had a strong basal texture and was also supported by the rapid activation of extension twins as predicted by crystal plasticity simulations. It is worth mentioning that the extension twins in magnesium usually grow quickly and encompass the entire grains. Thus, no evidence of a high volume fraction of extension twins was observed in the final microstructure of ND0ED0 and CA0ED0 samples at the end of the mechanical test.The samples in the second group (ND90ED0, CA0ND0, and TA90ND0) have similar flow stress with falling strain-hardening rates as shown in Figure 2. Unlike the samples in the first group, the extension twins are not expected to occur in these samples [10, 13, 14, 33]. For ND90ED0 and CA0ND0 samples, the compression direction aligns parallel to the crystallographic c-axis in most grains (Figures 5(a) and 5(d)), and therefore the <a> slips in basal or prismatic slip systems would not be activated. In this case, the second hard pyramidal <c+a> system is the dominant slip mode that can accommodate the plastic strain. Hence, it is not surprising to see high strain-hardening rates in the early portion of the stress-strain curves. The activation of compression twins is an additional possible deformation mechanism reported in the literature for such deformation modes [11, 20, 21]. For TA90ND0 sample in the second group, since the tensile axis TA is parallel to the basal plane of most grains (Figure 5(e)), the prismatic <a> and second pyramidal <c+a> slips have more favorable orientations compared to extension twin and basal slip [34, 35].It is known that the activation of extension twins in magnesium has significant effects on the texture evolution since tensile twins reorient the crystal by 86° around12¯10 directions. This can be clearly seen in the texture evolution of the samples in the first group (ND0ED0 and CA0ED0). These samples begin the mechanical tests with the loading direction aligned nearly perpendicular to the crystallographic c-axis in most grains (Figures 5(b) and 5(c)). However, the activation of tensile twins after about 0.03 axial strains alters the orientations of the crystallographic c-axis in these samples to become nearly parallel to the loading direction forming the strong basal texture ({0001}-fiber texture) shown in Figures 6(a) and 7(a). On the other hand, the texture of ND90ED0 and CA0ND0 samples in the second group with no active tensile twins did not change significantly (Figures 6(b) and 7(b)) compared to their initial basal texture shown in Figures 5(a) and 5(d). These samples were expected to deform mainly by slip in the second pyramidal <c+a> system. It is interesting to note that even though the initial texture of the samples in the two groups before the mechanical tests is different (the loading direction is perpendicular to the c-axis in ND0ED0 and CA0ED0 samples and parallel to the c-axis in ND90ED0 and CA0ND0 samples), the final texture for these samples is relatively similar. It is also seen that the flow response and strain-hardening rates of the samples in the first group (ND0ED0 and CA0ED0) beyond a small strain of 0.08 (the point at which extension twins deplete or saturate) are roughly similar to the behavior of the ND90ED0 and CA0ND0 samples in the second group. This can be related to the similarity in texture and deformation mechanism, where the texture of the ND0ED0 and CA0ED0 samples at 0.08 strain transforms mainly to a basal texture, which is similar to the initial texture of the ND90ED0 and CA0ND0. The samples with such a fiber texture deform essentially by slip in the second pyramidal <c+a> system.


### 4.2. Crystal Plasticity Simulation

In order to examine the effect of texture on the activation of various slip and twin systems, the different experimental deformation modes tested in this study were modeled using crystal plasticity constitutive equations. In this model, the plastic velocity gradient tensor is assumed to have contributions from both slip and deformation twinning [36, 37]:(1)LP=∑αNsγ˙αm0α⊗n0α+∑βNtwf˙βγtwcheck spacem0β⊗n0β,where γ˙α represents the shearing rate on the slip system α, γtw describes the constant shear associated with the twin system, and fβ denotes its volume fraction in the given crystal. m0α(β) and n0α(β) denote the slip (twin shear) direction and the slip (twin) plane normal to the slip (twin) system α(β), respectively, in the initial configuration. Ns and Ntw refer to the total number of available slip and twin systems, respectively. The evolution of the plastic shearing rate γ˙α and the deformation twin volume fraction f˙β are quantified using power-law relations [38, 39]:(2)γ˙α=γ˙0ταssα1/msignτα,(3)f˙β=γ˙0γtwτβstwβ1/m,where γ˙0 is a reference value of slip rate taken as 0.001 s−1 and m is a rate-sensitivity parameter taken as 0.02 and assumed to be the same for both slip and twinning. τα⋅(τβ) and ssα⋅(stwβ) denote the resolved shear stress and slip (twin) resistance of the α-slip (β-twin) system, respectively. An extended version of the saturation-type hardening laws has been employed in this work to describe the slip-hardening evolution [36–38]:(4)S˙sα=hs1+C∑βNtwfβ1−ssαssatαacheck space∑kNsγ˙k,ssatα=ss0+sprcheck space∑βNtwfβ0.5,where the parameters hs, C, ssatα, a, and spr aim to capture the phenomenological slip-hardening behavior due to slip and twin interactions. More details about the physical interpretations of these parameters can be found in [36–38]. The twin resistance was assumed to be constant in the current simulation.A total of 18 slip systems (3 basal <a>, 3 prism <a>, 6 pyramidal <a>, and 6{102¯2}check space second pyramidal <c+a>) and 6 {101¯2}check space1¯011 tensile twin systems were considered as listed in Table 1. The simulations were conducted assuming that the slip resistance values for the basal, prismatic, pyramidal <a>, second pyramidal <c+a>, and tensile twin modes were in the ratio of 1 : 50 : 50 : 80 : 7 [30, 40, 41]. The values of the various slip and twin hardening parameters in (4) were established by fitting the measured stress-strain curves of CA0ND0, CA0ED0, and TA90ND0 samples. These mechanical tests were selected because they produced the most distinct mechanical behavior (Figure 2) since the loading directions align either parallel or perpendicular to the crystallographic c-axis in most grains in these samples. The initial texture in the samples was captured using a set of 1000 discrete crystal orientations selected randomly from the initial measured texture. Five loading conditions were simulated corresponding to plane strain compression, simple compression, and uniaxial tensile tests with the loading direction aligned either along the crystallographic c-axis or along a direction normal to the c-axis. These loading conditions cover the five sets of experiments performed in this study: namely, ND90ED0 (ND ∥c-axis), ND0ED0 (ND ⊥c-axis), CA0ND0 (CA ∥c-axis), CA0ED0 (CA ⊥c-axis), and TA90ND0 (TA ⊥c-axis). Figure 8 shows the predicted final textures after an effective strain of about 0.15 for the five set of mechanical tests performed in this study using Taylor-type crystal plasticity simulations. It can be seen that the predicted textures are in reasonable agreement with the measured-final textures shown in Figures 6 and 7. The deviation between the simulated and measured textures appears in the pole figure of the plane strain compressed samples where the simulated results do not produce the strong-basal texture obtained in the measured texture (Figure 6).Table 1
Various slip and twin systems used in the crystal plasticity simulations for Mg.
System
Shear direction,m0
Plane normal,n0
Basal <a> slip
0
0
0
1
−2
1
1
0
0
0
0
1
1
−2
1
0
0
0
0
1
1
1
−2
0
Prism <a> slip
0
1
−1
0
−2
1
1
0
1
0
−1
0
1
−2
1
0
1
−1
0
0
1
1
−2
0
Pyramidal <a> slip
0
1
−1
1
−2
1
1
0
0
−1
1
1
−2
1
1
0
1
0
−1
1
1
−2
1
0
−1
0
1
1
1
−2
1
0
1
−1
0
1
1
1
−2
0
−1
1
0
1
1
1
−2
0
Second pyramidal <c+a> slip
1
1
−2
2
−1
−1
2
3
−1
2
−1
2
1
−2
1
3
−2
1
1
2
2
−1
−1
3
−1
−1
2
2
1
1
−2
3
1
−2
1
2
−1
2
−1
3
2
−1
−1
2
−2
1
1
3
Extension twin
0
−1
1
2
0
1
−1
1
1
0
−1
2
−1
0
1
1
−1
1
0
2
1
−1
0
1
0
1
−1
2
0
−1
1
1
−1
0
1
2
1
0
−1
1
1
−1
0
2
−1
1
0
1Figure 8
Predicted final texture ({0001} pole figures) using crystal plasticity simulations under the five set of mechanical tests performed in this study: plane strain compression (ND0ED0 and ND90ED0), simple compression (CA0ED0 and CA0ND0), and uniaxial tension (TA90ND0).Our main interest is to capture the most favored (active) slip or twin modes for the various loading conditions performed in this study using crystal plasticity simulation. Figure9 shows the relative slip and twin activities for plane strain compression of pure Mg with the loading direction aligned either parallel (ND90ED0 sample) or perpendicular (ND0ED0 sample) to the crystallographic c-axis in most grains. The predictions from the crystal plasticity model suggest that much of the plastic strain in ND90ED0 sample (ND ∥ c-axis) is accommodated by the second pyramidal <c+a> slip. However, the plastic strain in ND0ED0 sample (ND ⊥c-axis) is initially accommodated by the activation of tensile twins up to a strain of about 0.07, a point at which most grains are oriented such that the c-axis becomes nearly parallel to the loading, beyond that the second hard pyramidal <c+a> slip is activated to accommodate the further deformation. In the case of simple compression, Figure 10 presents the most favored deformation mechanisms predicted by the crystal plasticity simulation. These results are somehow similar to PSC for the two sample orientations tested experimentally in this study: namely, CA0ND0 (CA ∥c-axis) and CA0ED0 (CA ⊥c-axis). More specifically, Figure 10(a) shows that the second pyramidal <c+a> slip is the predominant active deformation mechanism in CA0ND0 sample (CA ∥c-axis) for the entire deformation. However, for CA0ED0 sample (Figure 10(b)), it is seen once again that the tensile twins are initially operated and then followed by the activation of the second pyramidal <c+a> slip (after a strain of about 0.07). Finally, for the uniaxial tensile test (TA90ND0 sample), the crystal plasticity predictions shown in Figure 11 suggest that the plastic strain is most likely accommodated by the prismatic <a> and second pyramidal <c+a> slips. These results are expected for uniaxial tension when the basal plane is parallel to the tensile axis [34, 35].Figure 9
Relative slip and twin activities predicted by crystal plasticity simulation during plane strain compression of pure Mg for the two sample orientations tested in this study: (a) ND90ED0 sample (ND ∥ c-axis) and (b) ND0ED0 sample (ND ⊥c-axis).
(a)
(b)Figure 10
Relative slip and twin activities predicted by crystal plasticity simulation during simple compression of pure Mg for the two sample orientations tested in this study: (a) CA0ND0 (CA ∥c-axis) and (b) CA0ED0 (CA ⊥c-axis).
(a)
(b)Figure 11
Relative slip and twin activities predicted by crystal plasticity simulation during uniaxial tension of pure Mg with the basal plane parallel to the tensile axis, TA90ND0 sample.


## 4.1. Deformation Mechanisms and Texture Evolution

The deformation behavior of magnesium and its alloys can be very different depending on the initial texture and the type of plastic deformation mechanisms. The mechanical deformation behavior of the pure magnesium tested in this study as described by the flow curves shown in Figure2 can be divided into two groups each having a special characteristic in common. The first group (ND0ED0 and CA0ED0 samples) has a concave stress-strain curve consisting of an initial drop in hardening rate followed by a marked increase up to a true strain of 0.1 where another drop is noticed. On the other hand, the second group (ND90ED0, CA0ND0, and TA90ND0 samples) has the typical falling strain-hardening rates similar to the saturation-type hardening observed in high and medium staking fault energy face-centered cubic metals. The offset yield point for the samples in the first group is about 75 MPa, whereas it is about 130 MPa for the samples in the second group. All flow curves in both groups have a sharp loss of strength at the end of the mechanical tests accompanied by the occurrence of inhomogeneous deformation in terms of shear cracks. It should be clear that this classification does not assume the same deformation mechanisms for all samples within the same group. The deformation mechanisms and texture evolution for the samples in each group are discussed next.For the first group consisting of ND0ED0 and CA0ED0 samples, the loading (compression) direction in the mechanical tests (PSC for ND0ED0 and simple compression for CA0ED0) aligns nearly perpendicular to the crystallographic c-axis in most of the grains as the texture indicated in Figures 5(b) and 5(c). This type of deformation mode is known to activate the {101¯2}check space<1¯011> tensile twins extensively at the early stage of deformation [21, 30]. This twin mode reorients the crystal by 86° around <12¯10> directions bringing the basal plane to be nearly perpendicular to the loading direction. The formation of this basal texture ({0001}-fiber texture) leads to texture hardening where most grains become hard to deform with the c-axis nearly parallel to the loading axis. The {101¯2}check space tensile twinning in our samples appears to start at an axial true strain of about 0.03 in ND0ED0 and CA0ED0 samples. This strain value corresponds to the points at which the strain-hardening rates start to increase as illustrated in Figures 2 and 3 (marked using arrows, note that Figure 2 shows the equivalent strain). This observation is in agreement with previously reported experiments in Mg alloys that detected the initial appearance of tensile twins at strains in the range of 0.02 to 0.04 [21, 31]. It should be clear that, in this study, we have not tried to detect the initial appearance of tensile twins experimentally at this small-strain level since this would require interrupting the mechanical tests frequently. The production of extension twins and the accompanied texture hardening lead to the sharp increase in flow stresses only after a small amount of true strain, where the flow stresses increased from about 80 MPa to about 160 MPa by a true strain of about 0.07 in both ND0ED0 and CA0ED0 samples. Therefore, the very high strain-hardening rates observed in samples ND0ED0 and CA0ED0 can be attributed to the activation of extension twins and the associated texture-hardening effect. Furthermore, the exhaustion of extension twins in these samples is estimated to occur at about 0.08 strain. This strain value corresponds to the peak strain-hardening rates shown in Figure 3. It should be noted that this value is very close to the maximum shear strain that can be produced by extension twins, estimated to be 0.1289 shear strain (or about 0.065 axial strain). In other words, most of the grains in ND0ED0 and CA0ED0 samples are expected to form a basal texture due to extension twins at about 0.08 strain. The plastic strain after the exhaustion of extension twins beyond an equivalent strain of about 0.08 is most likely accommodated by the second pyramidal <c+a> slip. Once slip and/or twinning cannot accommodate the shape change, a sudden loss of strength occurred due to the development of an inhomogeneous deformation in terms of shear bands that leads to softening and eventually cracking.It is interesting to note that the flow responses and strain-hardening rates in ND0ED0 and CA0ED0 samples are quite similar in spite of the difference in the applied boundary conditions for the two mechanical tests. To explain this, let us first consider the geometrical constraints in the two mechanical tests. In simple compression with CA normal to the c-axis, the sample has no lateral constraint that would suppress the extension of crystals along the c-axis as a result of extension twins. However, in plane strain compression with ND normal to the c-axis, the crystals with c-axis parallel to TD are not expected to twin because twinning would cause broadening in the constrained direction, TD. Thus, only the crystals with c-axis not perfectly aligned with TD can twin (especially those with c-axis parallel to the free direction, RD). In other words, the amount of twin activities is expected to be higher in simple compression compared to plane strain compression when the loading direction aligns perpendicular to the c-axis. This observation has been numerically investigated using crystal plasticity models for ideal orientations in prior literature (e.g. [20, 30, 32]). However, in the current study, it is noted that the flow responses and strain-hardening rates of the plane strain compressed sample (ND0ED0) are quite similar to the simple compressed sample (CA0ED0). This can be attributed to the existence of many grains in ND0ED0 sample with enough scatter of the c-axis about TD and therefore the easy activation of extension twins compared to other slip systems. This is clearly evidenced by the fact that the final texture of ND0ED0 sample at the end of the mechanical test had a strong basal texture and was also supported by the rapid activation of extension twins as predicted by crystal plasticity simulations. It is worth mentioning that the extension twins in magnesium usually grow quickly and encompass the entire grains. Thus, no evidence of a high volume fraction of extension twins was observed in the final microstructure of ND0ED0 and CA0ED0 samples at the end of the mechanical test.The samples in the second group (ND90ED0, CA0ND0, and TA90ND0) have similar flow stress with falling strain-hardening rates as shown in Figure 2. Unlike the samples in the first group, the extension twins are not expected to occur in these samples [10, 13, 14, 33]. For ND90ED0 and CA0ND0 samples, the compression direction aligns parallel to the crystallographic c-axis in most grains (Figures 5(a) and 5(d)), and therefore the <a> slips in basal or prismatic slip systems would not be activated. In this case, the second hard pyramidal <c+a> system is the dominant slip mode that can accommodate the plastic strain. Hence, it is not surprising to see high strain-hardening rates in the early portion of the stress-strain curves. The activation of compression twins is an additional possible deformation mechanism reported in the literature for such deformation modes [11, 20, 21]. For TA90ND0 sample in the second group, since the tensile axis TA is parallel to the basal plane of most grains (Figure 5(e)), the prismatic <a> and second pyramidal <c+a> slips have more favorable orientations compared to extension twin and basal slip [34, 35].It is known that the activation of extension twins in magnesium has significant effects on the texture evolution since tensile twins reorient the crystal by 86° around12¯10 directions. This can be clearly seen in the texture evolution of the samples in the first group (ND0ED0 and CA0ED0). These samples begin the mechanical tests with the loading direction aligned nearly perpendicular to the crystallographic c-axis in most grains (Figures 5(b) and 5(c)). However, the activation of tensile twins after about 0.03 axial strains alters the orientations of the crystallographic c-axis in these samples to become nearly parallel to the loading direction forming the strong basal texture ({0001}-fiber texture) shown in Figures 6(a) and 7(a). On the other hand, the texture of ND90ED0 and CA0ND0 samples in the second group with no active tensile twins did not change significantly (Figures 6(b) and 7(b)) compared to their initial basal texture shown in Figures 5(a) and 5(d). These samples were expected to deform mainly by slip in the second pyramidal <c+a> system. It is interesting to note that even though the initial texture of the samples in the two groups before the mechanical tests is different (the loading direction is perpendicular to the c-axis in ND0ED0 and CA0ED0 samples and parallel to the c-axis in ND90ED0 and CA0ND0 samples), the final texture for these samples is relatively similar. It is also seen that the flow response and strain-hardening rates of the samples in the first group (ND0ED0 and CA0ED0) beyond a small strain of 0.08 (the point at which extension twins deplete or saturate) are roughly similar to the behavior of the ND90ED0 and CA0ND0 samples in the second group. This can be related to the similarity in texture and deformation mechanism, where the texture of the ND0ED0 and CA0ED0 samples at 0.08 strain transforms mainly to a basal texture, which is similar to the initial texture of the ND90ED0 and CA0ND0. The samples with such a fiber texture deform essentially by slip in the second pyramidal <c+a> system.


## 4.2. Crystal Plasticity Simulation

In order to examine the effect of texture on the activation of various slip and twin systems, the different experimental deformation modes tested in this study were modeled using crystal plasticity constitutive equations. In this model, the plastic velocity gradient tensor is assumed to have contributions from both slip and deformation twinning [36, 37]:(1)LP=∑αNsγ˙αm0α⊗n0α+∑βNtwf˙βγtwcheck spacem0β⊗n0β,where γ˙α represents the shearing rate on the slip system α, γtw describes the constant shear associated with the twin system, and fβ denotes its volume fraction in the given crystal. m0α(β) and n0α(β) denote the slip (twin shear) direction and the slip (twin) plane normal to the slip (twin) system α(β), respectively, in the initial configuration. Ns and Ntw refer to the total number of available slip and twin systems, respectively. The evolution of the plastic shearing rate γ˙α and the deformation twin volume fraction f˙β are quantified using power-law relations [38, 39]:(2)γ˙α=γ˙0ταssα1/msignτα,(3)f˙β=γ˙0γtwτβstwβ1/m,where γ˙0 is a reference value of slip rate taken as 0.001 s−1 and m is a rate-sensitivity parameter taken as 0.02 and assumed to be the same for both slip and twinning. τα⋅(τβ) and ssα⋅(stwβ) denote the resolved shear stress and slip (twin) resistance of the α-slip (β-twin) system, respectively. An extended version of the saturation-type hardening laws has been employed in this work to describe the slip-hardening evolution [36–38]:(4)S˙sα=hs1+C∑βNtwfβ1−ssαssatαacheck space∑kNsγ˙k,ssatα=ss0+sprcheck space∑βNtwfβ0.5,where the parameters hs, C, ssatα, a, and spr aim to capture the phenomenological slip-hardening behavior due to slip and twin interactions. More details about the physical interpretations of these parameters can be found in [36–38]. The twin resistance was assumed to be constant in the current simulation.A total of 18 slip systems (3 basal <a>, 3 prism <a>, 6 pyramidal <a>, and 6{102¯2}check space second pyramidal <c+a>) and 6 {101¯2}check space1¯011 tensile twin systems were considered as listed in Table 1. The simulations were conducted assuming that the slip resistance values for the basal, prismatic, pyramidal <a>, second pyramidal <c+a>, and tensile twin modes were in the ratio of 1 : 50 : 50 : 80 : 7 [30, 40, 41]. The values of the various slip and twin hardening parameters in (4) were established by fitting the measured stress-strain curves of CA0ND0, CA0ED0, and TA90ND0 samples. These mechanical tests were selected because they produced the most distinct mechanical behavior (Figure 2) since the loading directions align either parallel or perpendicular to the crystallographic c-axis in most grains in these samples. The initial texture in the samples was captured using a set of 1000 discrete crystal orientations selected randomly from the initial measured texture. Five loading conditions were simulated corresponding to plane strain compression, simple compression, and uniaxial tensile tests with the loading direction aligned either along the crystallographic c-axis or along a direction normal to the c-axis. These loading conditions cover the five sets of experiments performed in this study: namely, ND90ED0 (ND ∥c-axis), ND0ED0 (ND ⊥c-axis), CA0ND0 (CA ∥c-axis), CA0ED0 (CA ⊥c-axis), and TA90ND0 (TA ⊥c-axis). Figure 8 shows the predicted final textures after an effective strain of about 0.15 for the five set of mechanical tests performed in this study using Taylor-type crystal plasticity simulations. It can be seen that the predicted textures are in reasonable agreement with the measured-final textures shown in Figures 6 and 7. The deviation between the simulated and measured textures appears in the pole figure of the plane strain compressed samples where the simulated results do not produce the strong-basal texture obtained in the measured texture (Figure 6).Table 1
Various slip and twin systems used in the crystal plasticity simulations for Mg.
System
Shear direction,m0
Plane normal,n0
Basal <a> slip
0
0
0
1
−2
1
1
0
0
0
0
1
1
−2
1
0
0
0
0
1
1
1
−2
0
Prism <a> slip
0
1
−1
0
−2
1
1
0
1
0
−1
0
1
−2
1
0
1
−1
0
0
1
1
−2
0
Pyramidal <a> slip
0
1
−1
1
−2
1
1
0
0
−1
1
1
−2
1
1
0
1
0
−1
1
1
−2
1
0
−1
0
1
1
1
−2
1
0
1
−1
0
1
1
1
−2
0
−1
1
0
1
1
1
−2
0
Second pyramidal <c+a> slip
1
1
−2
2
−1
−1
2
3
−1
2
−1
2
1
−2
1
3
−2
1
1
2
2
−1
−1
3
−1
−1
2
2
1
1
−2
3
1
−2
1
2
−1
2
−1
3
2
−1
−1
2
−2
1
1
3
Extension twin
0
−1
1
2
0
1
−1
1
1
0
−1
2
−1
0
1
1
−1
1
0
2
1
−1
0
1
0
1
−1
2
0
−1
1
1
−1
0
1
2
1
0
−1
1
1
−1
0
2
−1
1
0
1Figure 8
Predicted final texture ({0001} pole figures) using crystal plasticity simulations under the five set of mechanical tests performed in this study: plane strain compression (ND0ED0 and ND90ED0), simple compression (CA0ED0 and CA0ND0), and uniaxial tension (TA90ND0).Our main interest is to capture the most favored (active) slip or twin modes for the various loading conditions performed in this study using crystal plasticity simulation. Figure9 shows the relative slip and twin activities for plane strain compression of pure Mg with the loading direction aligned either parallel (ND90ED0 sample) or perpendicular (ND0ED0 sample) to the crystallographic c-axis in most grains. The predictions from the crystal plasticity model suggest that much of the plastic strain in ND90ED0 sample (ND ∥ c-axis) is accommodated by the second pyramidal <c+a> slip. However, the plastic strain in ND0ED0 sample (ND ⊥c-axis) is initially accommodated by the activation of tensile twins up to a strain of about 0.07, a point at which most grains are oriented such that the c-axis becomes nearly parallel to the loading, beyond that the second hard pyramidal <c+a> slip is activated to accommodate the further deformation. In the case of simple compression, Figure 10 presents the most favored deformation mechanisms predicted by the crystal plasticity simulation. These results are somehow similar to PSC for the two sample orientations tested experimentally in this study: namely, CA0ND0 (CA ∥c-axis) and CA0ED0 (CA ⊥c-axis). More specifically, Figure 10(a) shows that the second pyramidal <c+a> slip is the predominant active deformation mechanism in CA0ND0 sample (CA ∥c-axis) for the entire deformation. However, for CA0ED0 sample (Figure 10(b)), it is seen once again that the tensile twins are initially operated and then followed by the activation of the second pyramidal <c+a> slip (after a strain of about 0.07). Finally, for the uniaxial tensile test (TA90ND0 sample), the crystal plasticity predictions shown in Figure 11 suggest that the plastic strain is most likely accommodated by the prismatic <a> and second pyramidal <c+a> slips. These results are expected for uniaxial tension when the basal plane is parallel to the tensile axis [34, 35].Figure 9
Relative slip and twin activities predicted by crystal plasticity simulation during plane strain compression of pure Mg for the two sample orientations tested in this study: (a) ND90ED0 sample (ND ∥ c-axis) and (b) ND0ED0 sample (ND ⊥c-axis).
(a)
(b)Figure 10
Relative slip and twin activities predicted by crystal plasticity simulation during simple compression of pure Mg for the two sample orientations tested in this study: (a) CA0ND0 (CA ∥c-axis) and (b) CA0ED0 (CA ⊥c-axis).
(a)
(b)Figure 11
Relative slip and twin activities predicted by crystal plasticity simulation during uniaxial tension of pure Mg with the basal plane parallel to the tensile axis, TA90ND0 sample.


## 5. Conclusions

The anisotropic stress-strain responses of strongly textured pure magnesium were investigated at room temperature under three different monotonic deformation paths: (1) plane strain compression; (2) simple compression; and (3) uniaxial tensile tests. The samples in these mechanical tests were oriented such that the loading direction aligned either parallel or normal to the crystallographicc-axis in most of the grains. The main conclusions of this study are summarized below based on our experimental observations and supported by the predicted relative slip and twin activities using crystal plasticity simulations:(a)
The activation of{101¯2}check space tensile twins and the accompanying texture-hardening effect were associated with the initial increase in strain-hardening rate observed in samples subjected to plane strain compression or simple compression with the loading direction normal to the crystallographic c-axis.(b)
The transition in the stress-strain response from the initial increase in strain-hardening rates to the falling hardening rates was correlated with the exhaustion of deformation twins and the activation of slips along the second pyramidal <c+a> systems.(c)
The slip on the second hard pyramidal <c+a> system was found to be the dominant deformation mechanism that can accommodate the plastic strains in samples subjected to plane strain compression or simple compression with the loading direction parallel to the crystallographicc-axis.(d)
The prismatic <a> and second pyramidal <c+a> slips are the main deformation mechanisms during uniaxial tension when the basal plane is parallel to the tensile axis.



---

*Source: 2867281-2018-01-30.xml*