Plane matching across grain boundaries

[Ben Schuessler]

Hi All,

I was wondering if there is a way to determine which plane(s) are most closely "aligned" across grain boundaries? 

I have looked at the mPrime tutorial, however, if I understood it correctly, that requires some sort of loading on the sample to determine the active slip systems. I am looking to see if there is a way to do this without a loading scenario.

Could it be done with looking at the orientation and length of a given plane trace in the xy-plane? Looking for some advice. Long live MTEX!

Thanks,

Ben

[ruediger Kilian]

Hi Ben,

I'm not really sure what the exact goal is, but assuming you have two orientations (could be grain meanOrientation, orientations across boundaries ...) o1 and o2, then

h2 = inv(o2)*o1 * h1

where h1 is a crystal direction/plane with respect to o1 and h2 will be another crystal direction/plane with respect to o2 which is parallel to h1.

Is that what you are looking for? Or are you searching for the same crystal direction / plane which have the smallest angle to another?

Cheers,
Rüdiger

[Ben Schuessler]

Hey Rudiger!

Thanks for the equation! 

However, just like the second question you proposed, I'm looking to find the crystal direction/plane that has the smallest angle to each other. Determining which planes are most parallel with each other. 

Thanks,

Ben

[Rüdiger Kilian]

Hi Ben,
I think this sort of does what you are looking for, though it's looking for pairs of Millers on a grid, I think might be a better solution:

mtexdata titanium
grains = ebsd.calcGrains;
%%
o = ebsd(grains.boundary('t','t').ebsdId).orientations; % orientations along gb
g = Miller(equispacedS2Grid('resolution',3*degree),ebsd('T').CS); % grid of Millers
% minimum angle between all pairs of directions
[mang, id] = min(angle(o(:,1)*g,o(:,2)*g,'antipodal'),[],2);
g(id) %should be the Millers most closely aligned
%%
cK = ipfHSVKey(ebsd('T').CS);
plot(grains.boundary('t','t'),cK.Miller2Color(g(id)),'linewidth',3)
nextAxis

%% for grain meanOrientation
og = grains('id',grains.boundary('t','t').grainId).meanOrientation
[mang, idg] = min(angle(og(:,1)*g,og(:,2)*g,'antipodal'),[],2);
plot(grains.boundary('t','t'),cK.Miller2Color(g(idg)),'linewidth',3)
nextAxis
plot(cK)


So obviously, the result strongly depends on th choice of the grid, but maybe it's a starter.

Cheers,
Rüdiger



________________________________________
From: mtex...@googlegroups.com <mtex...@googlegroups.com> on behalf of Ben Schuessler <schuess...@gmail.com>
Sent: Friday, August 2, 2019 5:02:35 PM
To: MTEX
Subject: Re: {MTEX} Plane matching across grain boundaries

--
If you want to reduce the number of emails you get through this forum login to https://groups.google.com/forum/?fromgroups=#!forum/mtexmail, click "My membership" and select "Don't send me email updates". You can still get emails on selected topics by staring them.
---
You received this message because you are subscribed to the Google Groups "MTEX" group.
To unsubscribe from this group and stop receiving emails from it, send an email to mtexmail+u...@googlegroups.com<mailto:mtexmail+u...@googlegroups.com>.
To view this discussion on the web visit https://groups.google.com/d/msgid/mtexmail/f9205c5a-e0aa-4bd8-a2f7-df2f46d7ee5c%40googlegroups.com<https://groups.google.com/d/msgid/mtexmail/f9205c5a-e0aa-4bd8-a2f7-df2f46d7ee5c%40googlegroups.com?utm_medium=email&utm_source=footer>.

[Rüdiger Kilian]

Hi Ben,
actually it should be the plane normal to the misorientations axis, forget the grid:
pp = grains.boundary('t','t').misorientation.axis;
Cheers,
Rüdiger

[Ben Schuessler]

Hey Rudiger!

Sorry for just now getting back to you. Your suggestion was perfect! I'll reply again once I have time to play around with this topic a little more.

Cheers, mate. 

Ben