FCC single crystal overshoot

[bharat]

Hi all,

I am trying to plot the inverse pole figure for fcc single crystal under tension. I am using the below code. My problem is what ever the euler angle it is going into the same triangle. For overshoot i need two slip systems to be active. i.e., under tension the texture will overshoot from one slip system to another and go towards 111. Could anyone please guide me how to plot this.

I tried changing crystal symmetry to 432. It is just giving me the mirror image of other slip system.


CS = crystalSymmetry('m-3m', [4.0496 4.0496 4.0496], 'mineral', 'Aluminum', 'color', 'k');
SS = specimenSymmetry('-1');

path = '/home/mtex-5.0.3';
texture = {'/ivp.txt'};

[tt,nfiles] = size(texture);
disp(nfiles);

plota2east;
for i=1:nfiles
    fname = fullfile(path,texture{i});

    test = load_generic(fname,'header',0,'Columns',[1,2,3]);
    [m,n] = size(test);
    clear o;
    for j=1:m
        o(j) = orientation('Euler',test(j,1)*degree,test(j,2)*degree,test(j,3)*degree,CS,SS);
    end
  
   plotIPDF(o,[[xvector],[yvector],[zvector]],'MarkerSize',6,'MarkerEdgeColor','K','MarkerFaceColor','R')
  hold off
end

[ruediger Kilian]

Hi,
plotIPDF projects crystal directions into the fundamental sector by default. You can try to use the option 'complete' which will show the full hemispheres with all symmetrically equivalent directions or you can try if plot(inv(o)*xvector) does what you expect.

Hope this helps.
Cheers,
Rüdiger

[bharat]

Hi Ruediger,

Thank you so much for your reply. But this is not i am expecting. Please see the attached figure for overshoot. with complete option (figure attached) I am able to see all slip systems but they are mere reflections of one another. I need to plot my texture orientation which moves from one slip system to another slip system. But with the above code(in my previous mail) all my euler angles are going into only one triangle (as you mentioned which is the fundamental region). So, as i understand from your reply i need to define my fundamental region as two slip systems (which primary and Conjugate). Can you give me any suggestion regarding this.

[Rüdiger Kilian]

Hi Bharat,

did you try to plot the crystal directions (to my understanding, what we are plotting in your case are crystal directions; a slip system is something different) with plot(inv(o)*reference_direction) ?

Cheers,

Rüdiger

[bharat]

Hi Rudiger,

Please see the attached file for plot(inv(o)*zvector). I don't see this is also useful. And you are right i am plotting the motion of a single crystal. Basically in that overshoot picture it shows the motion of a single crystal when loaded in 125 direction. First it will move towards the 111-110 line and it overshoots because of hardening of the conjugate slip system and moves towards the stable orientation 111. That is why i need to plot the primary and its conjugate slip systems and show how my crystal is oriented according to the loading direction. May i get any suggestion for this. Thanks in advance.

Bharat

[Rüdiger Kilian]

Hi Bharat,
your plot isn't labeled, so it's hard to tell if it shows what you like to see.
e.g.
cs=o.CS
hs =[Miller(0,0,1,cs) Miller(1,1,0,cs) Miller(1,1,1,cs)]
annotate(hs,'symmetrised','labeled')

If this shows that crystal directions are actually where they should be, I do not see why

o.CS = o.CS.properGroup
plotIPDF(o,somerefdir)

wouldn't show the same, just maybe in another fundamental sector.

Cheers,
Rüdiger

[bharat]

Hi Rudiger,

Sorry for the late reply. Your are right.


 o.CS = o.CS.properGroup
 plotIPDF(o,[zvector],'MarkerSize',3)

the above code is showing what i need but i want to treat both triangles separately i.e., no reflection (please see the attachment). and the i tried the plot(inv(o)*zvector,'MarkerSize',3)  also. It is showing me upper and lower hemispheres separately.As per my understanding o.CS.properGroup is showing crystal symmery related to 432 and  I am not able to understand what plot(inv(o)) it is actually showing (sorry, i am not able to annotate on this plot). Please correct me if i am wrong.


Thanks in advance.

Regards,
Bharat

[Rüdiger Kilian]

Hi Bharat,
for a crystal with orientation o, inv(o)*vector3d(x,y,z) gives the crystal direction parallel to the specimen direction (x,y,z). That is what is under the hood of plotIPDF.

I assume you're having something still not correct in your script, since the following does not give orientations with inversion:

cs=crystalSymmetry('m-3m');
o=orientation.rand(1,cs)
plotIPDF(o,zvector)
figure
o.CS = o.CS.properGroup
plotIPDF(o,zvector)

Maybe share the revlevant parts of your script.

Cheers,
Rüdiger

[bharat]

Hi Rudiger,

Thank you so much. Please see the script. and attached file. Can you please tell me where it is going wrong. One more request. Can you please tell me how to connect 111 and 100 with a line.

CS = crystalSymmetry('m-3m', [4.0496 4.0496 4.0496], 'mineral', 'Aluminum', 'color', 'k');

SS = specimenSymmetry('triclinic');
path = '/home/mtex-5.0.3';
    fname = fullfile(path,'GAMMA.TXT');

    test = load_generic(fname,'header',0,'Columns',[1,2,3]);
    [m,n] = size(test);
    clear o;

    for j=1:m,

        o(j) = orientation('Euler',test(j,1)*degree,test(j,2)*degree,test(j,3)*degree,CS,SS);
    end
     

    figure(1)
    plotIPDF(o,[[xvector],[yvector],[zvector]],'MarkerSize',3,'MarkerEdgeColor','k','MarkerFaceColor','r')
%     hold off
%    
    figure(2)
    o.CS = o.CS.properGroup
    plotIPDF(o,zvector,'MarkerSize',3)


Regards,
Bharat

[Rüdiger Kilian]

Hi Bharat,
mind sharing one of those input files?

For the line this would be an option:

o.CS = o.CS.properGroup
plotIPDF(o,zvector,'MarkerSize',3)

hold on
line(Miller(0,0,1,cs),Miller(-1,1,1,cs),'Linewidth',3)
hold off

Cheers,
Rüdiger

[bharat]

Hi Rudiger,

I am sorry. I should have shared the input file in the previous mail. Please see the attached input file. Thank you,

Regards,
Bharat

[Rüdiger Kilian]

Hi,

the issue was that you SS defined as 'triclinic' atually gives '-1' instead of  '1', completeyl omitting the definition of SS works as well as it defaults to '1':

This:

SS = specimenSymmetry('1');
CS = crystalSymmetry('432');
  
ori = loadOrientation_generic('GAMMA.txt','CS',CS,'SS',SS, 'ColumnNames', ...
      {'Euler1' 'Euler2' 'Euler3'},'Columns',[1,2,3],'Bunge')
plotIPDF(ori,zvector,'MarkerSize',3)
hold on
line(h(1),h(2),'Linewidth',3)
hold off

does the attached figure.

Cheers,

Rüdiger

[bharat]

Thank you so much Rudiger. Your help is incredible. 

Regards,
Bharat