In this section we discuss properties of grains that are related to the distribution of orientations within the grains, i.e.,
grain orientation spread
grain average misorientation
grain average misorientation axis
As usual, we start by importing some EBSD data and computing grains
Grain average orientation
As by construction grains consist of pixels with similar orientation. In order to access all the orientations that belong to a specific grain we make use of the property
ebsd.grainId where we have stored during grain reconstruction for every EBSD pixel to which grain it belongs. Hence, we may use logical indexing on our EBSD variable
ebsd to find all orientations that belong to a certain
We could now use the command
mean to compute the grain average orientation from these individual orientations. A more direct approach, however, is to access the property
grain.meanOrientation which has been filled with the mean orientations during grain reconstruction.
Misorientation to the grain mean orientation
Once we have a reference orientation for each grain - the grain meanorientation
ori_ref - we may want to analyze the deviation of the individual orientations within the grain from this reference orientation. The grain reference orientation deviation (GROD) is the misorientation between each pixel orientation to the grain mean orientation defined as
While the above command computes the misorientations to the grain mean orientation just for one grain the command
calcGrod computes it for all grains simultaneously
Grain orientation spread (GOS)
The full misorientation to the mean orientation is often not so simple to understand. A good starting point is to just look at the misorientation angles to the grain mean orientation. The average of the misorientation angles to the grain mean orientation is called grain orientation spread (GOS) and can be computed by the command
grainMean which averages arbitrary EBSD properties over grains. Here, we use it to average the misorientation angle for each grain separately.
It should be noted that the GOS is also directly available as the grain property
grainMean can also be used to compute the maximum misorientation angle to the mean orientation within each grain. To achieve this we have to pass the function
@max as an additional argument. Note, that you can replace this function also with any other statistics like the median, or some quantile.
Grain average misorientation (GAM)
A measure that is often confused with the grain orientation spread is the grain average misorientation (GAM). The GAM is defined as the kernel average misorientation (KAM) averaged over each grain. Hence, we can compute is by using the functions
While the GOS measures global orientation gradients within a grain the GAM reflect the average local gradient.
The misorientation axis (crystal dispersion axis)
Under certain conditions, deformation may result in the dispersion of orientations within a grain. This can usually be the case when deformation is accommodated by slip on one dominant slip system for each grain and conditions are such, that the resulting orientation gradients are preserved in the material (as it is the case in many geomaterials deforming at moderate temperatures). In such a case, we would expect the orientations inside a grain to be aligned along a line with a specific misorientation axis to the mean orientation. Such a line is called fibre and we can use the command
fibre.fit to find the best fitting fibre for a given list of orientations. Lets do this for a single grain.
fibre.fit has three output arguments. The first one
f is the fitted fibre. From this we can easily detect the prominent misorientation axis in specimen coordinates by
f.r and in crystal coordinates by
The second output argument
lambda are the eigenvalues of the orientation matrix. The largest eigenvalue indicates are localized the orientations are. The second largest eigenvalue is a measure how much the orientation distributed along the fitted fibre. The third and forth eigenvalue describe how much the orientations scatter off the fibre. The scatter off the fibre is more conveniently described in the last output argument
fit, which is the mean misorientation angle of the orientations to the fitted fibre.
Lets perform the above analysis for all large grains
And plot the fit, the third and the second largest eigenvalues. We clearly see how the fit is related to the third largest eigenvalue \(\lambda_2\).
The crystal dispersion axes in crystal coordinates
In order to visualize the crystal dispersion axes we first need to define an appropriate color key for crystal directions. This can be done with the command
HSVDirectionKey. Note, that we need to specify the option
'antipodal' since for the crystal dispersion axes we can not distinguish between antipodal directions.
Now we can use this colorkey to visualize the misorientation axes in the grain map
The crystal dispersion axes in specimen coordinates
Colorizing the crystal dispersion axes in specimen coordinates is unfortunately much more complicated. In fact, it is mathematically impossible to find a corresponding color key without color jumps. Instead MTEX visualizes axes in specimen coordinates by compass needles which are entirely gray if in the plane and get divided into black and white to indicate which end points out of the plane and which into the plane.
In many materials, a direct relation can be observed between the position of the crystal dispersion axis in specimen coordinates and the inferred type of flow. E.g. in many geomaterials which have undergone (close to) simple shear progressive deformation, the average of the crystal dispersion axes align parallel to the vorticity axis of flow; in pure shear progressive deformation, crystal dispersion axes form a girdle with a normal parallel to the shortening direction.
to get some idea about any preferred direction, we can add contours, weighted by the fit. grains with a large mean misorientation angle will also have a more well defined direction of the dispersion axis.
Here we do not see this clear of a picture (maybe because this is a piece of steel which might behave differently, maybe because we do not consider a large enough number of grains) Question: if this is processed steel, which sample directions is pointing to the east?
TODO: Testing on Bingham distribution for a single grain
Although the orientations of an individual grain are highly concentrated, they may vary in the shape. In particular, if the grain was deformed by some process, we are interested in quantifications.
Testing on the distribution shows a gentle prolatness, nevertheless we would reject the hypothesis for some level of significance, since the distribution is highly concentrated and the numerical results vague.