Re: Orientation Relationship (between 2 phases)

[Ralf Hielscher]

As far as I can reproduce your problem

plotBoundary(grains,'property',rot,'delta',5*degree)

should also include misorientations between different phases. The only drawback I can see is that one can not restrict the output to specific phase transitions.

In a future release we plan to support the syntax where mori is a misorientation between two specific phases phase1 and phase2

mori = orientation('axis',xvector,'angle',60*degree,csPhase1,csPhase2)

in this case the command

plotBoundary(grains,'property',mori,'delta',5*degree)

should only plot grains boundaries between phase1 and phase2.

Can you please check this?

Ralf

[Robert Krakow]

Hey Ralf,

thank you for the response. I tried the following to extrapolate misorientation axis and angle of three misorientations:

cs01 = symmetry('hexagonal');
cs02 = symmetry('cubic');

%% class orientation

q_ex1 = euler2quat(57*degree,55*degree,48*degree,'Bunge');
q_ex2 = euler2quat(335*degree,15*degree,48*degree,'Bunge');
q_ex3 = euler2quat(139*degree,29*degree,37*degree,'Bunge');
or_ex1 = orientation(q_ex1,cs01);
or_ex2 = orientation(q_ex2,cs02);
or_ex3 = orientation(q_ex3,cs02);

angle(or_ex1,or_ex2)/degree
angle(or_ex3,or_ex1)/degree
angle(or_ex3,or_ex2)/degree

axis(or_ex1,or_ex2)
axis(or_ex3,or_ex1)
axis(or_ex3,or_ex2)

when using axis and angle for misorientation between orientation 1 and orientation 2:

figure('position',[200 100 1000 800])
plotBoundary(grains)
hold on
% abc are the entries for vector3d = axis(or1,or2)
rot4 =  rotation('axis',vector3d(0.466,0.133,0.874),'angle',56.0148*degree);
plotBoundary(grains,'property',rot4,'delta',5*degree)
legend('>2^\circ','56.01^\circ/[a b c]')
title('Misorientation: Angle-Axis-Pair 56.01°,(a,b,c)');

I only obtain the grain boundary map without the additional specified boundaries highlighted (see attachment).
However, there is also an error in orange:

Warning: Symmetry mismatch!
> In orientation.mtimes at 66
  In orientation.axis at 17
  In C_Class_Orientation at 29

Perhaps the misorientation axis has to be calculated differently?
E.g. when I follow the steps accordingly for crystallographic twins, I can obtain the right angle (i.e. 60 deg)
but the axis is 0,0,1 though I would expected 1,1,1.

 R.