Crystal Orientation as Coordinate Transformation edit page

In MTEX a crystal orientation is defined as the rotation that transforms crystal coordinates, i.e., a description of a vector or a tensor with respect to the crystal reference frame, into specimen coordinates, i.e., a description of the same object with respect to a specimen fixed reference frame.

In MTEX any orientation consists of two ingredients. A rotation

% lets take a random one
rot = rotation.rand
rot = rotation
 
  Bunge Euler angles in degree
     phi1     Phi    phi2
  79.3568 50.2577 294.557

and a description of the crystal lattice, which are represented in MTEX by variables of type crystalSymmetry

% lets take cubic crystal symmetry
cs = crystalSymmetry.load("Al-Aluminum.cif")
cs = crystalSymmetry
 
  mineral : Aluminum
  symmetry: m-3m    
  elements: 48      
  a, b, c : 4, 4, 4

Combining both ingredients allows us to define an orientation

ori = orientation(rot,cs)
ori = orientation (Aluminum → xyz)
 
  Bunge Euler angles in degree
     phi1     Phi    phi2
  79.3568 50.2577 294.557

As a consequence a variable of type orientation is at the same time of type rotation and hence allows for all operations that are available for rotations.

Crystal coordinates to specimen coordinates

Let us consider the following crystal direction

h = Miller(1,0,0,cs,'uvw')
h = Miller (Aluminum)
  u v w
  1 0 0

In a grain with orientation ori this direction h has the specimen coordinates

r = ori * h
r = vector3d
        x        y        z
  2.62519   1.2191 -2.83219

Similarly, orientations transform tensors given with respect to the crystal reference frame, e.g., the following single crystal stiffness tensor

C = stiffnessTensor(...
  [[2 1 1 0 0 0];...
  [1 2 1 0 0 0];...
  [1 1 2 0 0 0];...
  [0 0 0 1 0 0];...
  [0 0 0 0 1 0];...
  [0 0 0 0 0 1]],cs)
C = stiffnessTensor (Aluminum)
  unit: GPa              
  rank: 4 (3 x 3 x 3 x 3)
 
  tensor in Voigt matrix representation:
 2 1 1 0 0 0
 1 2 1 0 0 0
 1 1 2 0 0 0
 0 0 0 1 0 0
 0 0 0 0 1 0
 0 0 0 0 0 1

into a stiffness tensor with respect to the specimen reference frame

ori * C
ans = stiffnessTensor (xyz)
  unit: GPa              
  rank: 4 (3 x 3 x 3 x 3)
 
  tensor in Voigt matrix representation:
  2.4972  0.9427  0.5601  0.1377 -0.0851   -0.02
  0.9427  2.2007  0.8566  -0.247  0.0578  0.0626
  0.5601  0.8566  2.5832  0.1093  0.0273 -0.0426
  0.1377  -0.247  0.1093  0.8566 -0.0426  0.0578
 -0.0851  0.0578  0.0273 -0.0426  0.5601  0.1377
   -0.02  0.0626 -0.0426  0.0578  0.1377  0.9427

Objects that can be translated by orientations from crystal into specimen coordinates and vice versa include

Specimen coordinates into crystal coordinates

Conversely, we can go back from specimen coordinates to crystal coordinates by multiplying with the inverse orientation

inv(ori) * r
ans = Miller (Aluminum)
        h       k       l
  16.3991       0       0

Note, that in literature orientations are often defined to transform specimen coordinates into crystal coordinates, i.e., to coincide with the inverse orientations in MTEX. The consequences of this differences are exhaustively discussed in the topic orientation conventions.

Specimen Rotation

Rotations of the specimen ,i.e., changing the specimen reference frame, do also change the orientation. Assume the specimen is rotated about the X-axis about 60 degree. We may define this rotation by

rot = rotation.byAxisAngle(vector3d.X,60*degree);

Then an orientation ori is updated to the rotated reference frame by

ori_new = rot * ori
ori_new = orientation (Aluminum → xyz)
 
  Bunge Euler angles in degree
     phi1     Phi    phi2
  50.4216 78.6571 354.793

It should also be noted, that orientations are sensitive with respect to the alignment of the Euclidean reference frame \(\vec X\), \(\vec Y\), \(\vec Z\) with respect to the crystal axes \(\vec a\), \(\vec b\) and \(\vec c\). This issue is discussed in more detail in the topic Crystal Reference Frames.