Creuziger et al. Integrating Materials and Manufacturing Innovation 2014, 3:1

http://www.immijournal.com/content/3/1/1

RESEARCH Open Access

Crystallographic texture evolution in 1008 steel

sheet during multi-axial tensile strain paths

Adam Creuziger1*, Lin Hu2,3, Thomas Gnäupel-Herold1 and Anthony D Rollett2

*Correspondence:

adam.creuziger@nist.gov

1National Institute of Standards and

Technology, 100 Bureau Dr., 20899,

Gaithersburg, MD, USA

Full list of author information is

available at the end of the article

Abstract

This paper considers the crystallographic texture evolution in a 1008 low carbon steel.

The texture evolution along uniaxial, plane strain and balanced biaxial strain states

were measured. For uniaxial testing, grains tend to rotate such that the {111}110 ̄  slip

directions are aligned with the loading axis. For plane strain and balanced biaxial strain

states, the majority of grains are distributed with the {111} plane parallel to the sample

normal direction. Accompanying visco-plastic self consistent (VPSC) predictions of the

texture evolution were made along same strain paths and strain increments.

Comparing between the measured texture evolution and computational texture

evolution indicate that the VPSC model qualitatively predicts the measured texture

evolution, but the rate at which the texture evolution occurs is over predicted.

Keywords: Crystallographic texture; Steel; Neutron diffraction; Metal forming

Introduction

Steel has been in use for nearly 4000 years, and remains one of the main materials (along

with concrete, glass and asphalt) used around the world. A large majority of the steels

created are carbon steels, with a large proportion of these carbon steels strained after production to create a required shape. Accurate prediction of the final shape and mechanical

properties after straining continues to be an important industrial goal. One of the methods applied to this problem is measurement of the crystallographic texture (or preferred

orientation) as a function of strain path and strain increment. Accurate prediction of the

crystallographic texture and how the texture evolves will allow better estimates of the

flow stress and plastic contraction, which in turn provides a more accurate prediction of

the final shape and failure during forming. The particular case of deformation of sheet

stock material is discussed in this paper. The majority of the literature discusses the high

strain increments relevant for rolling sheet material. Unfortunately existing [1-6] experimental data on the texture evolution along multiaxial strain states at intermediate strain

increments () relevant for forming are incomplete, so assessment of model accuracy is

difficult in this regime.

This paper focuses on the measurement of the crystallographic texture evolution in

1008 AISI/SAE low carbon steel sheets after strain increments typical in forming structures from the rolled sheet. The texture is measured by neutron diffraction to interrogate a

large number of grains. Five different strain states were investigated: uniaxial strain along

© 2014 Creuziger et al.; licensee Springer. This is an Open Access article distributed under the terms of the Creative Commons

Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction

in any medium, provided the original work is properly cited.
Creuziger et al. Integrating Materials and Manufacturing Innovation 2014, 3:1 Page 2 of 19

http://www.immijournal.com/content/3/1/1

the rolling direction (URD), strain along the rolling direction with the transverse direction

fixed near zero strain (a plane strain condition, or PSRD), balanced biaxial strain (BB) in

the rolling and transverse directions, strain along the transverse direction with the rolling

direction fixed near zero strain (PSTD), and uniaxial strain along the transverse direction (UTD). In the plane of the sheet (normal direction or ND) the sample is assumed to

be in a plane stress condition. In addition, the texture evolution is predicted from the as

received texture using the visco-plastic self consistent (VPSC) model and compared with

the experimental data. Data on the texture evolution of low carbon steels will also provide

a well characterized reference for comparison to more complicated materials, such as high

strength steels (HSS) like high strength low alloy (HSLA) steel or advanced high strength

steels (AHSS) such as dual phase (DP) or transformation induced plasticity (TRIP) steels.

The data described in this paper has been submitted to the NIST/MatDL DSpace repository as part of the IMMI General Articles Collection [7]. The repository contains files

that include a description of the data contained in the repository, a schematic of the the

analysis and computational processes, some of the programs used to analyze the data, a

collection of the texture data files at various points in the analysis, and input files to the

VPSC model.

Background

Using the American Iron and Steel Institute (AISI) and Society for Automotive Engineers

(SAE) grade designation, 1008 steel is defined as an iron alloy possessing a mass fraction

of carbon ≤ 0.0010, a mass fraction of manganese between 0.0030 and 0.0050, a mass

fraction of phosphorus ≤ 0.0004, and a mass fraction of sulfur ≤ 0.0005 [8]. Owing to the

low carbon content and absence of austenitic stabilizers, the structure is almost entirely

ferritic (body centered cubic crystal structure) with some carbides. This allows for high

ductility (approximately 40% elongation), but low strength (tensile strength approximately

350 MPa) [9].

The texture evolution of steels has been widely studied in the past, particularly with

respect to the complete processing path from the melt to hot and cold rolling. These topics are of interest to create desirable textures for deep drawing. Some articles on these

topics include [2,3,10,11], which discuss steel textures and their origin from deformation

and recrystalization of the austenite phase and subsequent transformation from austenite

to ferrite. The complex details of the austenitic transformation are beyond the scope of

this paper, but it is important to note that deformation and recrystallization in the austenite phase are one of the primary causes of crystallographic texture in the ferritic phase.

The review by Inagaki [12] summarizes work on the inter-relationship between texture,

grain boundaries and dislocation structures in low carbon steels. Hutchinson [13] discusses texture control in low carbon steels as a function of processing parameters. Bodin

et al. [2] shows the strong effects of processing path and temperature on the texture.

The predicted plastic strain ratio (R value) has been used to correlate the effect of a particular texture component to the macroscopic properties. For deep drawing applications,

sheet material that has a large volume of grains aligned with the γ -fiber ({111} plane parallel to the ND) are preferred as they have a high average R value (R ̄) and a planar anisotropy

(R) value near 0 [10,13]. Orientations with the {100} plane parallel to the ND direction

such as ‘cube’ and ‘shear’ orientations are not desirable, since they have low average R values [10,13]. This is due to the orientation of slip systems for body centered cubic (BCC)
Creuziger et al. Integrating Materials and Manufacturing Innovation 2014, 3:1 Page 3 of 19

http://www.immijournal.com/content/3/1/1

materials. Conversely, for high strength materials where high ductility is not required,

orientations along the α-fiber (110 crystal direction parallel to rolling direction) are preferred for high strength [14] but may produce anisotropy and brittleness, particularly for

the {113}110 orientation [10].

There is also a significant amount of investigation into the textures caused by cold

rolling in the ferrite phase. This data is particularly relevant to this paper as cold rolling,

apart from a hydrostatic term, produces a stress state similar to plane strain along the

rolling direction (PSRD). As shown in a series of review papers, cold rolling ferrite

strengthens the α and γ fibers, particularly the α fiber at higher reduction ratios [3,10,11].

Figure thirteen in Kestens and Jonas [11] shows relative amounts of texture intensity as a

function of thickness reduction. Inagaki [12] discusses some of the microstructural reasons for this, and theorizes that grain pairs with {111}[uvw] type grain boundaries that

allow the grain to rotate easily, keeping orientations along the γ fiber stable.

Kestens and Jonas [11] proposes that {211}011 along the α fiber is the stable final orientation in plane strain compression, while {111}110 and {554}225 orientations are

caused by ferritic recrystallization after plane strain compression. Inagaki [12] proposes

that the {223}011 orientation is the final stable orientation in plane strain compression.

There is additional literature on cross-rolling of low carbon steel sheets [5], a strain state

similar to plane strain along the transverse direction (PSTD). Results from Wronski et al.

[5] show that the shear orientation and γ fiber are stable orientations after near 80%

thickness reduction.

There has been limited investigation into the texture evolution along other strain paths

relevant for sheet metal forming. Yerra et al. [6] is one of the most complete, following five

different strain paths from equal biaxial to uniaxial tension, all with the largest strain values in the rolling direction. The orientation distribution functions (ODFs) are measured

at a single level of strain for each strain path. In Yerra et al. [6], the experimental textures were not found to change greatly between the different strain modes nor the volume

fractions of specific orientations. However a calculation using the Lamel model by the

same authors show a significant difference between the strain paths. Other notable works

include Luzin et al. [15], which measured texture after balanced biaxial straining. The γ

fiber in a balanced biaxial strain mode is intensified, in agreement with Mirshams et al.

[16] and theoretical calculations that predict the θ (100 parallel to ND direction) and

γ fibers to be the stable orientations in uniaxial compression along the ND. Other strain

modes, such as torsion [17] and simple shear [18] are outside the scope of this paper.

As steel is such a relevant industrial material, there have been attempts to predict the

texture evolution. Savoie and Jonas [19] used a rate sensitive analysis to model the texture

along a flange during deep drawing. These simulations include a rotation rate map that

shows the stable end orientations for different constraints imposed on the models. They

found that all initial orientations tend to rotate toward the {223}110 orientation. Toth

et al. [20] used a rate dependent slip theory that contains models for pancake, lath and

Taylor models to calculate texture evolution in low carbon steels during rolling. Starting

with the ‘RGoss’ orientation, {111}110, and shear orientations; they found the {111}112

orientation also forms, but becomes unstable at higher strains. Many more stable orientations were found using this method, but most were observed along the α fiber and γ

fiber as well as {554}225. Raabe [21] used a Taylor type model to predict the rolling

texture in low carbon deep drawing steel as a function of reduction ratio. The Taylor
Creuziger et al. Integrating Materials and Manufacturing Innovation 2014, 3:1 Page 4 of 19

http://www.immijournal.com/content/3/1/1

model showed increases in the intensity of the α fiber and γ fiber in good agreement with

experiment [21].

Methods

Experimental

Chemical analysis

Chemical analysis was performed following ASTM E415 [22] for all elements except sulfur which followed ASTM E1019 [23]. The resulting elemental mass fractions are given

in Table 1. These are within the required parameters for AISI/SAE 1008 steel with the

exception of manganese which is slightly lower than specified.

Grain size measurements

The average grain size was measured according to ASTM E112 [24]. Linear intercept

methods along the rolling direction (RD), transverse direction (TD) and normal direction (ND) were performed along five lines in three images over two different samples

for all six cross-sectional views (N = 60 per direction). These averages were summed

together to find the overall average grain size of 18.7 μm and aspect ratios of 1.0:0.78:0.46

for RD:TD:ND, normalized to RD. Complete information on the grain size is shown in

Table 2.

Uniaxial sample testing

Samples were rough cut from 1.17 mm (0.046 in, 18 gauge) thick sheet stock into sample blanks approximately 19 mm (0.75 in) × 200 mm (8 in) using a hydraulic shear. A

TENSILKUTa machine was used to create dogbone samples that conform to ASTM E8

[25] standards for sheet specimens with a 50.8 mm (2 inch) gauge length and 12.7 mm

(0.5 inch) width. After machining, samples were loaded into a screw driven tensile frame

moving at a crosshead speed of 0.1 mm

sec resulting in a nominal strain rate of 0.002 strainsec .

Axial and transverse strains were measured using extensometers. The strains imparted to

the samples are listed in Table 3. The von Mises effective strain was calculated assuming

volume conservation and neglecting shear strains:

eff =

√2

3



(11 − 22)

2 + (22 − 33)2 + (33 − 11)2 1

2 (1)

Multi-axial sample testing

Samples for non-uniaxial stain paths were loaded in a custom testing machine [26] using

a sample geometry developed by Raghavan [27]. A clamp hold down force of -450 kN

(force control) and a ram rate of 0.1 mm

sec was applied resulting in a nominal strain rate of

0.2 strain

sec . Digital Image Correlation (DIC) was used to measure the plastic strains from a

circular area with an undeformed diameter of 32 mm (1.25 inches). Strains are listed in

Table 3.

Table 1 Mass fraction analysis performed by an independent lab per [22] and [23] for sulfur

Element Measured 1008 AISI-SAE carbon steel grade [8]

Carbon 0.0007 0.0010 max

Manganese 0.0029 0.0030 to 0.0050

Phosphorus 0.00010 0.00040 max

Sulfur 0.00014 0.00040 max

Iron Balance
Creuziger et al. Integrating Materials and Manufacturing Innovation 2014, 3:1 Page 5 of 19

http://www.immijournal.com/content/3/1/1

Table 2 Average grain sizes for 1008 steel using the linear intercept method

Average Standard Confidence Relative

Plane Direction grain size [μm] deviation [μm] interval [μm] accuracy [%]

ND RD 24.15 2.94 1.10 4.55

TD RD 28.72 3.25 1.21 4.22

ND TD 17.85 1.50 0.56 3.14

RD TD 23.40 2.27 0.85 3.62

RD ND 11.79 0.93 0.35 2.96

TD ND 12.59 0.78 0.29 2.32

Images were recorded after etching with Nital (2%) for 45 seconds. Images were recorded at 200 × magnification on all

three planes on two different samples each. n = 30 for each row.

Texture measurements

Crystallographic texture measurements were performed using neutron diffraction at the

BT-8 Residual Stress Diffractometer at NIST [28]. Strained and as received samples were

cut into 8 mm squares and stacked and glued together until a cube was formed. For the

average grain sizes calculated in Table 2, there are on the order of 2 × 107 grains in each

sample. Three complete pole figures from the (200), (220), and (211) reflections were measured. A 5° hexagonal grid [29] of the pole figure space was used, resulting in 469 separate

measurements, with the exception of the URD and UTD samples strained to 10% and 20%.

In these four samples a hexagonal grid of 6° spacing was used, resulting in 331 separate

measurements. The raw data was then processed in the program PF [30] and converted

to pole figures in the popLA format [31]. The as received samples exhibited orthotropic

Table 3 Applied strains, measured after unloading the samples, expressed as true strain

Name Nominal strain (%) RD TD ND eff

URD10 10% 0.100 -0.063 -0.037 0.101

URD20 20% 0.200 -0.127 -0.072 0.202

URD30 30% 0.302 -0.186 -0.116 0.305

PSRD10 10% 0.10 -0.01 -0.09 0.11

PSRD10 10% 0.10 -0.01 -0.09 0.11

PSRD20 20% 0.20 -0.02 -0.17 0.21

PSRD20 20% 0.20 -0.02 -0.17 0.21

PSRD30 30% 0.33 -0.03 -0.29 0.36

PSRD30 30% 0.33 -0.03 -0.30 0.37

BB10 10% 0.10 0.10 -0.19 0.19

BB20 20% 0.20 0.20 -0.40 0.41

BB30 30% 0.32 0.32 -0.63 0.63

PSTD10 10% -0.01 0.10 -0.09 0.11

PSTD10 10% -0.01 0.10 -0.09 0.11

PSTD20 20% -0.02 0.22 -0.20 0.24

PSTD20 20% -0.02 0.22 -0.20 0.24

PSTD30 30% -0.03 0.34 -0.31 0.37

PSTD30 30% -0.03 0.33 -0.30 0.37

UTD10 10% -0.063 0.100 -0.037 0.101

UTD20 20% -0.125 0.201 -0.076 0.203

UTD30 30% -0.167 0.297 -0.130 0.297

For uniaxial testing, strains are measured along the RD and TD by extensometers, for BB and PS testing RD and TD strains are

measured by digital image correlation. ND strains are calculated from volume conservation. The von Mises effective strain is

calculated as shown in equation 1.
Creuziger et al. Integrating Materials and Manufacturing Innovation 2014, 3:1 Page 6 of 19

http://www.immijournal.com/content/3/1/1

sample symmetry from the rolling process. All sample strains occurred aligned with the

initial orthotropic symmetry resulting in orthotropic symmetry in the strained samples.

For these reasons, orthotropic sample symmetry was applied to the pole figures and the

orientation distribution functions (ODFs).

Orientation Distribution Function (ODF) calculation

From the series of pole figures, orientation distribution functions (ODFs) were calculated

to determine the distribution of grains as a function of crystallographic orientation in

the bulk material. Two different ODF solvers were investigated, the WIMV method [32]

in popLA and mTex [33]. The calculation of the ODF is an inverse problem, which are

known to not necessarily have a unique solution. For each solver a range of different ODFs

can be calculated that match the pole figures depending on the, number of iterations of

the ODF solving loop, resolution of the ODF, and other input parameters. To facilitate

comparison with the VPSC ODF calculations, the WIMV method was used. Uncertainty

in ODF measurements is an ongoing measurement challenge currently under investigation [34]. For the purposes of this paper, an uncertainty of 10% in the ODF values is

assumed [35].

Good agreement (2% calculation error) with the pole figures was found after 12 iterations. Experimental and recalculated pole figures were qualitatively inspected to check

agreement as well. The ODF from the as received sample was used to calculate an input

file for the VPSC analysis. The ODF is expressed in a crystal orientation distribution file

that has a regular 5° grid in the Bunge Euler angles φ1, , φ2. Due to the crystal and sample symmetries, the ODF can be expressed in a subspace with the range 0° ≤ φ1 ≤ 90°,

0° ≤  ≤ 90°, 0° ≤ φ2 ≤ 90°, albeit with an additional three fold symmetry in this

subspace. A value for intensity is listed at each point in this grid, expressed in units of

‘multiples of a uniform distribution’ or MUD (often referred to as ‘multiples of a random

distribution’ MRD).

Excerpts from the ODF were also extracted along fiber textures. The fibers of interest

for this paper are the α fiber, γ fiber and θ fiber. The α fiber runs along the  axis, with

φ1 = 0° and φ2 = 45°. The γ fiber runs along the φ1 axis, with  = 55.4° and φ2 = 45°.

The γ1, γ3 and γ2, γ4 pairs are equivalent points by symmetry.The θ fiber runs along the

φ1 axis, with  = 0° and φ2 = 45°. Specific orientations of interest are listed in Table 4,

with a shorthand name, Miller indices, and Bunge Euler angles. Additionally, difference

ODF plots were created between adjacent strain steps. These difference ODFs show the

locations that grains are rotating to or away from. However, as a consequence of the nonlinear nature of the intensity scale of the ODF, caution should be applied to ascribing

quantitative significance to the values in the difference plots. To aid in the idenification of

the orientations listed in Table 4, a key showing the location of each orientation located

in the φ2 = 45° cross section is shown in Figure 1.

Computational

The viscoplastic self-consistent (VPSC) model developed by Lebensohn and Tomé [36]

was used to simulate the texture development. The VPSC model treats each grain as a viscoplastic ellipsoidal inclusion embedded in and interacting with a homogeneous effective

matrix (HEM). Under the HEM assumption, individual grain response varies as a function of its orientation and the stiffness of the grain-to-matrix interaction, depending on
Creuziger et al. Integrating Materials and Manufacturing Innovation 2014, 3:1 Page 7 of 19

http://www.immijournal.com/content/3/1/1

Table 4 Relevant orientations investigated

Name {hkl} uvw φ1  φ2

Cube {001} 100 45° 0° 45°

Shear {001} 110 ̄  0° 0° 45°

Goss {011} 100 90° 90° 45°

RGoss {011} 011 ̄  0° 90° 45°

Copper {112} 1 ̄11 ̄  90° 35.26° 45°

Brass {011} 211 ̄  54.74° 90° 45°

S {213} 3 ̄64 ̄  58.98° 36.70° 63.44°

{554} {554} 2 ̄25 ̄  90° 60.50° 45°

γ1 {111} 110 ̄  0° 54.74° 45°

γ2 {111} 121 ̄  30° 54.74° 45°

γ3 {111} 011 ̄  60° 54.74° 45°

γ4 {111} 1 ̄12 ̄  90° 54.74° 45°

α1 {115} 110 ̄  0° 15.79° 45°

α2 {113} 110 ̄  0° 25.24° 45°

α3 {112} 110 ̄  0° 35.26° 45°

α4 {223} 110 ̄  0° 43.31° 45°

All orientations are described by Miller Indices and Bunge Euler angles.

the choice of homogenization methods at grain level. Details of the model can be found

in [36].

A file containing 24 938 discrete orientations was created to approximate the calculated

ODF from the as received sample. This list of orientations was used as an input file to the

VPSC algorithm. Three other input files were also required: a file containing the applied

velocity gradients and applied Cauchy stresses; a file containing the unit cell parameters,

elastic moduli of the single crystal, slip planes, and hardening behavior of slip families; and

an input file containing the path and name of the subfiles, grain aspect ratio, hardening

model, and interaction model.

Figure 1 A key showing the Bunge Euler angle location of orientations described in Table 4 for the

φ2 = 45°cross section of ODF space. All orientations except the S orientation appear in this cross section.

The same tick mark intervals were used in the key shown above as Figures 2, 3, 4, and 5.
Creuziger et al. Integrating Materials and Manufacturing Innovation 2014, 3:1 Page 8 of 19

http://www.immijournal.com/content/3/1/1

For all simulations a velocity gradient is imposed on all components of the velocity

gradient tensor matching the macroscopic strains in Table 3 in the following manner.

Imposed velocity tensors were the RD, TD and ND strains normalized such that the

highest positive strain was set equal to 1. The strain increment was the value of strain

along the highest positive strain direction. All stresses calculated by the VPSC model

are therefore deviatoric stresses. Only the {110}111 and {112}111 families of slip

planes and directions were considered in this model as initial experimentation including the {123}111 family of planes did not show a significant effect on the texture.

The ellipsoid ratios used were based on the experimental data listed in Table 2 normalized along the rolling direction and having values of 1.0, 0.78 and 0.46 for the RD, TD

and ND respectively. Secant linearization and first nearest neighbor interaction models were chosen. The outputs used here from the VPSC algorithm include a list of

discrete orientations that have been rotated to new orientations due to the imposed

velocity gradient. These discrete orientations were binned into a regular grid to create

an ODF.

Results

Experimental textures

Three different techniques have been used to show how the experimentally measured

texture evolves as a function of applied strain. Figure 2 shows the intensities of the ODF

Figure 2 Experimentally determined orientation distributions for 1008 steel as a function of applied

strain path and strain value, including the as received (AR) data. Bunge Euler angle convention was

used, with only the φ2 = 45° sections shown. Tick marks are shown in 15° increments. The color key shows

the value of intensity in terms of multiples of a uniform distribution (MUD). The strains paths are abbreviated

as follows: URD - uniaxial along the rolling direction, PSRD - plane strain along the rolling direction,

BB - balanced biaxial strain, PSTD - plane strain along the transverse direction, UTD - uniaxial along the

transverse direction. Strain values are the nominal strains of the largest strain, exact values can be found in

Table 3. The key in Figure 1 shows the locations of particular orientations in the same style plot.
Creuziger et al. Integrating Materials and Manufacturing Innovation 2014, 3:1 Page 9 of 19

http://www.immijournal.com/content/3/1/1

Figure 3 Differences between experimentally determined orientation distributions. Differences are

taken between adjacent strain steps (10% - as received (AR), 20% - 10%, 30% - 20%). Red indicates increasing

intensity, blue indicates decreasing intensity.

Figure 4 Orientation distributions calculated by VPSC for 1008 steel as a function of applied strain

path and strain value. The conventions described in Figure 2 also apply here.
Creuziger et al. Integrating Materials and Manufacturing Innovation 2014, 3:1 Page 10 of 19

http://www.immijournal.com/content/3/1/1

Figure 5 Differences in orientation distributions calculated by VPSC. Adjacent strain states were

subtracted (10% - as received (AR), 20%-10%, 30%-20%). Red indicates increasing intensity, blue indicates

decreasing intensity.

on φ2 = 45° cross sections. Figure 3 is a plot of the differences between adjacent strain

steps (10% - as received, 20% - 10%, 30% - 20%) on the same φ2 = 45° cross sections

shown in Figure 2. The smaller strain step is subtracted from the larger strain step, so

positive values indicate an increase in intensity at the given orientation. The third technique extracts line profiles along a path with two fixed Euler angles. These line profiles

correspond to the α fiber (Figure 6 top), γ fiber (Figure 7 top), and θ fiber (Figure 8 top).

For brevity, the discussion of the URD deformation contains the relevant figure numbers for discussion of the fiber components, they will not be repeated for each strain

mode.

Figure 2 includes the as received texture, as well as the three nominal strain increments

along the five strain paths explored. The as received texture of the 1008 steel is typical for

a rolled body centered cubic (BCC) material for deep drawing applications, with strong

intensity along the γ fiber, and peaks along the 110 RD directions of the (111) plane. The

intensity values for α fiber orientations and the shear orientation are also larger than in

a uniform texture distribution. The intensities of cube, shear, Goss, RGoss, copper, brass,

and S orientations are low, as expected in a BCC material.

When samples are strained in uniaxial tension along the rolling direction (URD), the

ODF intensity along the γ1 and γ3 orientations increases, as shown in the first row in

Figure 2, Figure 3, and Figure 7. The ODF intensity along the {554}2 ̄25 ̄ , γ2, and γ4

decrease. The largest change in the ODF intensities are along these orientations occurs

while straining from 0% to 10%. The ODF intensity along all of the α fiber orientations
Creuziger et al. Integrating Materials and Manufacturing Innovation 2014, 3:1 Page 11 of 19

http://www.immijournal.com/content/3/1/1

Figure 6 Top: experimentally measured texture along the α fiber, only highest strain levels shown.

The width of the line corresponds to the assumed standard deviation of 10% of the value. Bottom:

computationally predicted by VPSC texture along the α fiber, only highest strain levels shown. The width of

the line corresponds to the assumed standard deviation of 10% of the value. The ordinate range between the

two figures is dissimilar, a thin black line was drawn in the computational data set to show the maximum

scale (MUD = 20) for the experimental data.

increase, shown in Figure 2, Figure 3, and Figure 6. Along the θ fiber the grains along

the shear orientation increase, shown in Figure 2, Figure 3, and Figure 8, but the ODF

intensity along the cube orientation remains the same.

The PSRD strain state is notable in that there are very few changes from the as received

structure, depicted in the second row in Figures 2 and 3. For all of the orientations
Creuziger et al. Integrating Materials and Manufacturing Innovation 2014, 3:1 Page 12 of 19

http://www.immijournal.com/content/3/1/1

Figure 7 Top: experimentally measured texture along the γ fiber, only highest strain levels shown.

The width of the line corresponds to the assumed standard deviation of 10% of the value. Bottom:

computationally predicted by VPSC texture along the γ fiber, only highest strain levels shown. The width of

the line corresponds to the assumed standard deviation of 10% of the value. The ordinate range between the

two figures is dissimilar, a thin black line was drawn in the computational data set to show the maximum

scale (MUD = 20) for the experimental data.

investigated, the ODF intensity changes only very slightly at each strain level. The ODF

intensity of orientations along the α fiber are basically unchanged. There are some modest

increases in the intensity along the γ fiber, particularly filling in the γ2 and γ4 orientations to make the intensity along the γ fiber slightly higher and more uniform in

intensity. The θ fiber shows the same trend as the samples deformed with uniaxial tension
Creuziger et al. Integrating Materials and Manufacturing Innovation 2014, 3:1 Page 13 of 19

http://www.immijournal.com/content/3/1/1

Figure 8 Top: experimentally measured texture along the θ fiber, only highest strain levels shown.

The width of the line corresponds to the assumed standard deviation of 10% of the value. Bottom:

computationally predicted by VPSC texture along the θ fiber, only highest strain levels shown. The width of

the line corresponds to the assumed standard deviation of 10% of the value. The ordinate range between the

two figures is dissimilar, a thin black line was drawn in the computational data set to show the maximum

scale (MUD = 6) for the experimental data. The inset figures show the details near φ1 = 45° where the data

plots overlap.

along the rolling direction with the shear orientation increasing and the cube orientation

unchanged.

The dominant texture evolution in the BB strain state is intensification of the γ fiber.

The 20% sample still shows the highest intensity along the 110 RD directions, but by

30% strain all orientations along the γ fiber are approximately equal and have an intensity
Creuziger et al. Integrating Materials and Manufacturing Innovation 2014, 3:1 Page 14 of 19

http://www.immijournal.com/content/3/1/1

of 15 to 16 times higher than a uniform distribution. Grains oriented along the α fiber

have decreased sharply outside the range of the γ fiber and θ fiber. The BB strain state

also shows the strongest changes along the θ fiber. The cube orientation increases from

a value of 0.5 MUD in the as received strain state to over 5 MUD after 30% BB straining.

Similarly, the shear orientation increases dramatically.

The PSTD strain state results in a texture evolution that is similar to the BB strain

state. Grains oriented along the γ fiber increase, particularly the γ2 and γ4 orientations,

to become nearly constant in intensity at 10 MUD by 30% strain. Grains oriented along

the α fiber decrease, except for those near the γ fiber. The ODF intensity along the lower

part of the α fiber (φ1 ≥ 60°) increases more than observed in the RD, PSRD or BB strain

states. There is an increase in the number of grains oriented along the shear orientation,

while the cube orientation remains unchanged.

Due to the symmetries in this projection of Euler space, the texture evolution of UTD

samples look like the corresponding URD ODF mirrored about the φ1 = 90° axis. Grains

oriented along the γ2 and γ4 orientations increase, while γ1 and γ3 intensity decreases.

Orientations along the α fiber decrease. There is some increase in the copper orientation,

which can be related by symmetry operations to an orientation near the α3 orientation.

Along the θ fiber the fraction of grains oriented along these directions are similar to

measurements of the URD and PSRD samples.

Computational textures

Similar to the experimental data shown in Figure 2, Figure 4 shows the calculated intensities of the ODF in φ2 = 45° cross sections from VPSC simulations. Figure 5 is a plot of

the differences between adjacent strain steps (10% - as received, 20% - 10%, 30% - 20%) in

the ODFs calculated with VPSC. Line profiles along the α fiber (Figure 6 bottom), γ fiber

(Figure 7 bottom) and θ fiber (Figure 8 bottom) are also shown. Similar to the experimental data, the discussion of the URD deformation contains the relevant figure numbers for

discussion of the fiber components, they will not be repeated for each strain mode.

The intensity scale for each of the computational figures was adjusted to accommodate

the higher intensity values calculated in the VPSC code. However the color scale in the

range of the experimental intensities are unchanged for Figure 4 and Figure 5 to allow

for comparison, additional colors were added for intensities outside of the experimental

range. For Figures 6, 7, and 8 a thin black line is used to depict the maximum in the

experimental line profile figures.

The VPSC as received texture shown in Figure 4 is similar to, but less smooth than

the experimental as received texture shown in Figure 2. Regions around γ2 and γ4 are

slightly depleted in intensity compared to the experimental data. Similarly, areas around

the copper and brass orientations are also lower than observed experimentally.

The main intensity changes for URD occur along γ1 and γ3. There are corresponding

decreases in the γ2 and γ4 orientations. As the amount of strain increases, the intensity

around the γ3 orientation increases in a very narrow band in the φ1 axis while broadening in the  axis at higher strains. The α fiber increases from the as received sample,

but with larger increases at  > 55° than at  < 55° observed in the experimental

data. Along the θ fiber, the shear orientation shows a sharp drop, observed in all VPSC

data sets. These observations are likely computational artifacts due to the low number of grains (2.5 × 105) compared to the experimental data (≈ 2 × 107). The θ fiber
Creuziger et al. Integrating Materials and Manufacturing Innovation 2014, 3:1 Page 15 of 19

http://www.immijournal.com/content/3/1/1

orientations show modest increases in the shear texture, but the cube texture is essentially

unchanged.

When the strain state changes to plane strain along the rolling direction (PSRD), the

texture evolution is similar to the VPSC URD sample. The γ1 through γ3 orientations

increase but the γ2 and γ4 orientations decrease. However, the {554} orientation, which is

close to the γ4 orientation, increases like the γ1 and γ3 orientations. Unlike experimental

data, orientations near the α3 and α4 orientations show large increases, particularly at 20%

strain, but the other α fiber orientations only show slight increases in ODF intensity. Similar to the α4 orientations, the shear orientation shows a large increase at 20% strain, and

then decreases back to the as received intensity at 30%. The cube and shear orientations

in the θ fiber are also largely unchanged.

Balanced biaxial (BB) strains show the sharpest texture evolution. The texture sharpening is primarily along the γ fiber. All orientations along γ fiber increase by similar

amounts for each strain step, such that by 30% strain, intensities of 80 MUD are observed.

The α fiber shows a sharp peak where it intersects with the γ fiber, but decreases in all

other orientations. The θ fiber shows a large increase as well, predominantly along the

shear orientation.

Plane strain along the transverse direction (PSTD) shows increases in the ODF intensity

at all orientations along the γ fiber after 10% strain, and then increases in the γ2 and γ4

thereafter. Similar to the URD sample, the α fiber increases more at  > 55° than at

 < 55°. The shear orientation increases in intensity, but the cube orientation remains

nearly constant.

Finally, the UTD samples show intensity increases along γ2 and γ4, the 110 TD directions of the (111) plane, but of lower magnitude than the PSTD strain state. The intensity

of orientations γ1 and γ3 decrease. Like the URD case, as strain increases the highest

intensity bands tend to align themselves with the  axis. The intensity along the α fiber

orientations generally decrease, with only a slight increase at  = 65° to 70°. Orientations

near the shear texture also show large increases.

Discussion

Qualitatively there is good agreement between the experimental and computational texture evolution models. Most changes from the texture evolution have the correct ‘sense’

to them, experimental increases are matched to computational increases and vice versa.

Both the computational and experimental textures show significant changes even at

relatively low (10%) values of strain. There are some subtle differences between the Experimental and VPSC as received textures. As discussed in the results and methods sections,

the VPSC as received texture is less smooth, likely due to the three order of magnitude

reduction in the number of grains explored in the VPSC model (24 938 orientations) than

in the experimental data (estimated at ≈ 2 × 107). For orientations with φ2 = 45° and

 < 20° there are less than 10 grains for each of the orientations in the 5° bin and at  = 0

most of the orientations only have 1 grain associated with each bin. This may also be the

cause of oddities around the shear orientation in Figure 8, where large peaks are observed

5° from the orientation that are not observed experimentally in Figure 8.

There are a few places where there is qualitative disagreement between the trends in

the VPSC predictions and experimental data. In the experimental URD and PSRD strain

states there are broad increases in the α fiber, particularly at  < 55°, shown in Figure 6.
Creuziger et al. Integrating Materials and Manufacturing Innovation 2014, 3:1 Page 16 of 19

http://www.immijournal.com/content/3/1/1

The VPSC data in Figure 6 shows some increase in the intensity at  < 55°, but a larger

increase in intensity at  > 55°. Similarly the increase in cube orientation over the shear

orientation observed experimentally in Figure 8 under BB straining was not observed

in the VPSC predictions in Figure 8. At higher strain levels the VPSC model seemed to

rotate grains away from the γ fiber to orientations ±5° to ±10° away from  = 55°

(γ fiber), as shown by the decreases in intensity along the γ fiber in the difference map in

Figure 5. This may also affect the PSTD results data, where the intensity pattern observed

experimentally in Figure 7 does not match the VPSC predictions shown in Figure 7.

Quantitatively, the VPSC model tends to over predict the rate of texture evolution leading to much higher intensities and sharper orientation distributions. Overall, the intensity

values for experimental data at 30% strain better match the VPSC data at 10% strain. The

BB 30% strain state is a particularly good example of the difference in rate of texture evolution, with the intensity predicted by the VPSC over-predicted by a factor of 5. Another

example is the PSRD strain state, where the α4 and {554} orientations are predicted to

have increases larger than are observed experimentally.

Comparing this data with prior computational and analytical work, the uniaxial samples

show the 110 slip directions of the {111} plane as stable orientations, and instability of

slip directions 211 in the {111} plane. The PSRD strain state is indeed similar to the

predictions for cold rolling the ferrite phase. The experimental data indicates that while

the α4 ({223}110 ̄ ) orientation is the most stable, there seems to be a broad distribution

of orientations that are increasing at the strain levels investigated. The stability of the α

fiber changes with more transverse direction strain, such that in the BB and PSTD strain

states this fiber is not as stable. This result could be used in production, for applications

where α fiber is not desirable, as more lubrication on the rollers or cross rolling could be

used to remove some of the α fiber and increase the number of grains with orientations

aligned with the γ fiber; albeit with some increase in the cube and shear components as

well. The increase in γ fiber agrees well with observations by [15] and [16] but in our work

more cube orientated grains are observed than by Luzin et al. [15]. The experimental and

VPSC data are more consistent with the Lamel predictions than experimental data from

Yerra et al. [6].

A similar study 5754 aluminum compared Los Alamos polycrystal plasticity (LApp),

VPSC and visco-plastic fast Fourier transform (VPFFT) models for multiaxial strain states

[37]. Hu et al. [37] showed similar trends to those observed here, with the rate of texture evolution over predicted in the VPSC model, particularly for the BB strain state. The

VPFFT model calculated a rate of texture evolution closer to those observed experimentally. In the VPFFT model, the inter-grain and intra-grain interactions are included in the

model, likely improving the accuracy. However, the VPFFT model described by Hu et al.

[37] is discretized as voxel elements. Approximately 500 voxels were required for each

grain, as opposed to the single ellipsoid per grain in the VPSC model, resulting in the

VPFFT model containing one fifth as many grains as the VPSC model for a similar computational time. The reduction in number of grains resulted in a lower fidelity to the as

received texture. The texture in the 5754 aluminum studied was much closer to a uniform

distribution than the steel texture measured here, so a larger number of grains would

be required to accurately model the steel microstructure. Following work by Dawson

et al. [38], the addition of the elastic component to the viscoplastic models (EVPSC or

EVPFFT approaches) is not expected to improve the model as the strains in this study are
Creuziger et al. Integrating Materials and Manufacturing Innovation 2014, 3:1 Page 17 of 19

http://www.immijournal.com/content/3/1/1

relatively large strains along monotonically increasing strain paths. Additionally, the construction of an accurate microstructure with respect to the spatial correlation of adjacent

grains required as an input to a VPFFT or crystal plasticity based finite element analysis

(CP-FEA) models is a topic of current research.

However, there is some concern that even accurate microstructure may be insufficient

to model the deformation of polycrystalline materials. Pokharel et al. [39] compare a

modeled microstructure in VPFFT that is taken from experimental data from near field

high-energy X-ray diffraction microscopy (nf-HEDM) [40]. Even with a microstructure

model that closely matches the experimental data, the agreement between experiment

and model are not as close as desired. Similar to the results shown in this paper, there

is reasonable qualitative agreement, but there are significant qualitative and quantitative

variations. Pokharel et al. [39] conclude that the “empirically derived constitutive relations used in continuum scale simulations” (i. e. the Voce model in VPFFT and VPSC)

as the likely cause of the disagreement and recommend further work on the development of a constitutive relation based on dislocation interactions. Pokharel et al. [39] also

touches on CP-FEA, but states that the empirically derived continuum level constitutive relationships used in FEA may suffer similar problems as observed in the VPFFT

model.

Conclusions

To summarize, the crystallographic texture evolution in a 1008 low carbon steel along

five different strain states have been experimentally measured and computationally

predicted using VPSC. Comparisons between measured texture evolution and computational texture evolution indicates that the VPSC model is generally qualitatively accurate

in predicting the orientations to which grains rotate. Quantitatively, the VPSC model over

predicts the rate at which grains will rotate, particularly in multi-axial strain states such

as balanced biaxial.

Endnote

aCertain commercial equipment, instruments, or materials are identified in this paper

in order to specify the experimental procedure adequately. Such identification is not

intended to imply recommendation or endorsement by the National Institute of

Standards and Technology, nor is it intended to imply that the materials or equipment

identified are necessarily the best available for the purpose.

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

AC calculated the ODFs, compared the data and lead writing of the manuscript; LH conducted the VPSC simulations;

TGH measured and analyzed the experimental pole figures; ADR advised LH in VPSC, reviewed our work and assisted in

interpretation. All authors read and approved the final manuscript.

Acknowledgements

The authors thank Steve Banovic, Tim Foecke, Mark Iadicola and Mark Stoudt at the NIST Center for Metal Forming for

advice and assistance, and the Research Associate Program at the National Academies for supporting this work.

Author details

1National Institute of Standards and Technology, 100 Bureau Dr., 20899, Gaithersburg, MD, USA. 2Carnegie Mellon

University, 5000 Forbes Avenue, 15213, Pittsburgh, PA, USA. 3IBM Semiconductor Research and Development Center,

2070 Route 52, 12533, Hopewell Junction, NY, USA.

Received: 30 September 2013 Accepted: 20 November 2013

Published: 7 January 2014
Creuziger et al. Integrating Materials and Manufacturing Innovation 2014, 3:1 Page 18 of 19

http://www.immijournal.com/content/3/1/1

References

1. Bleck W, Grossterlinden R, Lotter U, Reip CP (1991) Textures in steel sheets. Steel Res 62(12): 580–586. ISSN 0177-4832

2. Bodin A, Sietsma J, van der Zwaag S (2002) Texture and microstructure development during intercritical rolling of

low-carbon steels. Metallurg Mater Trans A 33(6): 1589–1603. ISSN 1073-5623

3. Hutchinson B (1999) Deformation microstructures and textures in steels. Philos Trans R Soc Lond A

357(1756): 1471–1485. ISSN 1364-503X

4. Kestens L, Jonas JJ (1999) Deep drawing textures in low carbon steels. Metals And Mater-Korea 5(5): 419–427.

ISSN 1225-9438

5. Wronski S, Wrobel M, Baczmanski A, Wierzbanowski K (2013) Effects of cross-rolling on residual stress, texture and

plastic anisotropy in f.c.c. and b.c.c. metals. Mater Characterization 77: 116–126. ISSN 1044-5803

6. Yerra SK, Vankudre HV, Date PP, Samajdar I (2004) Effect of strain path and the magnitude of prestrain on the

formability of a low carbon steel: On the textural and microtextural developments. J Eng Mater Technol-Trans ASME

126(1): 53–61. ISSN 0094-4289

7. Creuziger A, Hu L, Gnäupel-Herold T, Rollett AD (2013) Data citation: crystallographic texture evolution in 1008 steel

sheet during multi-axial tensile strain paths. http://hdl.handle.net/11115/231 Accepted 11 Dec 2013

8. Oberg E, Jones FD, Horton HL, Ryffel H H (1996) Machinery’s handbook, 25th edition. Industrial Press Inc., 32

Haviland Street, Unit 2C | South Norwalk, CT 06854

9. (2013) MatWeb. AISI 1008 steel, CQ, DQ, and DQSK sheet. http://www.matweb.com/search/datasheet.aspx?

MatGUID=145867c159894de286d4803a0fc0fe0f

10. Ray RK, Jonas JJ (1990) Transformation textures in steels. Int Mater Rev 35(1): 1–36. ISSN 0950-6608

11. Kestens L, Jonas JJ (2006) ASM handbook, volume 14A:685–700. ASM, ASM International Materials Park, OH:

44073-0002

12. Inagaki H (1994) Fundamental aspect of texture formation in low-carbon steel. ISIJ Int 34: 313–321

13. Hutchinson B (1994) Practical aspects of texture control in low carbon steels. Mater Sci Forum 157–162: 1917–1928

14. Inagaki H, Kurihara K, Kozasu I (1977) Influence of crystallographic texture on strength and toughness of

control-rolled high-tensile strength steel. Trans Iron And Steel Institute of Japan 17(2): 75–81. ISSN 0021-1583

15. Luzin V, Banovic S, Gnaupel-Herold T, Prask HJ, Ricker RE (2005) Measurement and calculation of elastic properties in

low carbon steel sheet. In: Houtte PV, Kestens L (eds) ICOTOM 14: Textures Of Materials, volume 495-497 of Materials

Science Forum, pp 1591–1596

16. Mirshams RA, Crosby KE, Mohamadian HP, Burris CL (1993) The influence of biaxial stretching on the texture of

ultra-low carbon (ulc) steel sheets. Scripta Metallurgica Et Materialia 29(4): 433–438. ISSN 0956-716X

17. Baczynski J, Jonas JJ (1996) Texture development during the torsion testing of alpha-iron and two IF steels. Acta

Mater 44(11): 4273–4288. ISSN 1359-6454

18. Kang JY, Bacroix B, Regle H, Oh KH, Lee HC (2007) Effect of deformation mode and grain orientation on

misorientation development in a body-centered cubic steel. Acta Mater 55 15: 4935–4946. ISSN 1359-6454

19. Savoie J, Jonas JJ (1994) Simulation of the deformation textures induced by deep-drawing in extra low-carbon steel

sheets. Acta Metallurgica Et Materialia 42(12): 4101–4116. ISSN 0956-7151

20. Toth LS, Jonas JJ, Daniel D, Ray RK (1990) Development of ferrite rolling textures in low-carbon and extra low-carbon

steels. Metallurgical Trans A 21(11): 2985–3000. ISSN 0360-2133

21. Raabe D (1995) Simulation of rolling textures of b.c.c. matals considering grain interactions and crystalline slip on

{100}, {112} and {123} planes. Mater Sci Eng A 197: 31–37

22. ASTM International (2008a) Standard test method for atomic emission vacuum spectrometric analysis of carbon and

low-alloy steel. ASTM Standards E 415

23. ASTM International (2011) Standard test methods for determination of carbon, sulfur, nitrogen, and oxygen in steel,

iron, nickel, and cobalt alloys by various combustion and fusion techniques. ASTM Standards E 1019

24. ASTM International (2004) Standard test methods for determining average grain size E 112. ASTM Standards E 112

25. ASTM International (2008b) Standard test methods for tension testing of metallic materials E8/ E8M. ASTM

Standards E 8

26. Foecke T, Iadicola MA, Lin A, Banovic SW (2007) A method for direct measurement of multiaxial stress-strain curves

in sheet metal. Metallurgical and Materials Transactions A 38A: 306–313

27. Raghavan KS (1995) A simple technique to generate in-plane forming limit curves and selected applications.

Metallurgical and Mater Trans A 26A: 2075–2084

28. Brand PC, Prask HJ, Gneaupel-Herold T (1997) Residual stress measurements at the NIST reactor. Physica B

241–243: 1244–1245

29. Matthies S, Wenk HR (1992) Optimization of texture measurements by pole figure coverage with hexagonal grids.

Physica Status Solidi (a) 133: 253–257

30. Gnaeupel-Herold T (2012) Bt8 data analysis. http://www.ncnr.nist.gov/instruments/bt8/BT8DataAnalysis.htm.

Accessed 29 Jan 2013

31. Kallend UF, John S, Kocks, Rollett AD, Wenk HR (1991) popLA - an integrated software system for texture analysis.

Textures and Microstructures 14: 1203–1208

32. Matthies S, Vinel GW (1982) On the reproduction of the orientation distribution function of texturized samples from

reduced pole figures using the conception of a conditional ghost correction. Physica Status Solidi B

112(2): K111–K114. ISSN 0370–1972

33. Hielscher R, Schaeben H (2008) A novel pole figure inversion method: specification of the MTEX algorithm.

J Appl Crystallography 41: 1024–103

34. Creuziger A, Syed K, Gnaeupel-Herold T (2014) Measurement of uncertainty in orientation distribution function (odf)

calculations. Scripta Mat 72-73: 55–58

35. Bunge HJ, Haessner F (1968) 3-dimensional orientation distribution function of crystals in cold-rolled copper.

J Appl Phys 39(12): 5503–5515

36. Lebensohn RA, Tomé CN (1993) A self-consistent anisotropic approach for the simulation of plastic deformation and

texture development of polycrystals: application to zirconium alloys. Acta Metallurgica et Materialia 41: 2611–2624
Creuziger et al. Integrating Materials and Manufacturing Innovation 2014, 3:1 Page 19 of 19

http://www.immijournal.com/content/3/1/1

37. Hu L, Rollett AD, Iadicola M, Foecke T, Banovic S (2012) Constitutive relations for AA 5754 based on crystal plasticity.

Metallurgical and Mater Trans A 43A: 854–869

38. Dawson PR, MacEwen SR, Wu PD (2003) Advances in sheet metal forming analyses: dealing with mechanical

anisotropy from crystallographic texture. Int Mater Rev 48: 86–122

39. Pokharel R, Lind J, Kanjarala AK, Lebensohn RA, Li SF, Kenesei P, Suter RM, Rollett AD (2013) Polycrystal plasticity:

comparison between grain scale observations of deformation and simulations. submitted to Annual Reviews

Condensed Matter Physics

40. Li SF, Lind J, Hefferan CM, Pokharel R, Lienert U, Rollett AD, Suter RM (2012) Three-dimensional plastic response in

polycrystalline copper via near-field high-energy X-ray diffraction microscopy. J Appl Crystallography 45: 1098–1108

doi:10.1186/2193-9772-3-1

Cite this article as: Creuziger et al.: Crystallographic texture evolution in 1008 steel sheet during multi-axial tensile

strain paths. Integrating Materials and Manufacturing Innovation 2014 3:1.

Submit your manuscript to a 

journal and benefi t from:

7 Convenient online submission

7 Rigorous peer review

7 Immediate publication on acceptance

7 Open access: articles freely available online

7 High visibility within the fi eld

7 Retaining the copyright to your article

 Submit your next manuscript at 7 springeropen.com