How do I rotate my ebsd/odf data (Euler angles) to plot the pole figures of the loading direction in mtex

635 views
Skip to first unread message

Vishwanath Hegadekatte

unread,
Sep 16, 2014, 10:44:15 AM9/16/14
to mtex...@googlegroups.com
Hello Everyone,

I am new to mtex and also to texture analysis. At the moment I am trying to pickup crystal plasticity modeling. I am now faced with the plotting the pole figures using mtex. I use the following commands to plot the pole figures:

~~~
cd E:\mtex\mtex-3.5.0
startup_mtex
cs = symmetry('aluminum.cif')
ss = symmetry('triclinic')
fname = fullfile('E:UAT\','','texture.txt')
myebsd = loadEBSD_generic(fname,'cs',cs,'ss',ss,'bunge','degree','ColumnNames',{'Euler 1' 'Euler 2' 'Euler 3'},'header',2,'Columns', [1,2,3])
myodf = loadODF_generic(fname,'cs',cs,'ss',ss,'bunge','degree','ColumnNames',{'Euler 1' 'Euler 2' 'Euler 3'},'header',2,'Columns', [1,2,3])
plotpdf(myebsd,Miller(1,1,1))
figure
plotpdf(myodf,Miller(1,1,1))
~~~

I am using the aluminum.cif file from http://code.ohloh.net/file?fid=7pHHj8onBv8LniDjh960wl34gPo&cid=89dAYL7Sr-E&s=&fp=304052&mp&projSelected=true#L0
The file texture.txt is as follows: *----- Orientations -----* ang1 ang2 ang3 101.41 145.88 -112.30 135.67 88.12 -132.22 ... ... total of 500 Euler angles My crystal plasticity code creates the texture.txt file assuming that the loading axis is z-axis. But I have defined the geometry such that the loading axis is y-axis and therefore would like to rotate myebsd and myodf such that I am able to plot the pole figures of the loading direction. I have followed the instructions at http://merkel.zoneo.net/RDX/index.php?n=Texture.PlotIndividualOrientationsInMTex rot = rotation('Euler',0*degree,90*degree,0*degree) myodfRotated = rotate(myodf,rot) When I do the above I get the following error: Error using subsindex Function 'subsindex' is not defined for values of class 'EBSD'. Error in quaternion/subsref (line 10) q.a = subsref(q.a,s); What is it that I am doing wrong? Thank you all in advance. Vishwanath

Ralf Hielscher

unread,
Sep 22, 2014, 2:20:05 AM9/22/14
to mtex...@googlegroups.com
Hi Vishwanath,

as far as I can see everything looks fine. Could you please post which command exactly gives the error?

Thank you,

Ralf.

Vishwanath Hegadekatte

unread,
Sep 25, 2014, 2:16:06 AM9/25/14
to mtex...@googlegroups.com
Hi Ralf,

Thank you very much for your reply.

The following command: 

myodfRotated = rotate(myodf,rot)

results in the following error message:

Error using subsindex
Function 'subsindex' is not defined for values of class 'EBSD'.

Error in quaternion/subsref (line 10)
  q.a = subsref(q.a,s);

Vishwanath

PS: I am attaching the texture.txt file. Also I am using mtex version 3.5.0 on matlab version R2013a (8.1.0.604)
texture.txt

Ralf Hielscher

unread,
Oct 18, 2014, 4:22:31 PM10/18/14
to mtex...@googlegroups.com
Hi Vishwanath Hegadekatte,

I checked this on my computer and I can not reproduce the error.

Could you please post the output of the following commands

myodf

and

myebsd

Ralf.

On Tuesday, September 16, 2014 4:44:15 PM UTC+2, Vishwanath Hegadekatte wrote:
Reply all
Reply to author
Forward
0 new messages