Define HCP Slip Systems

[Karthik Palaniappan]

I tried to define Slip system for calculating the Taylor and Schmid Factor map (EBSD data)

r = vector3d(1,-1,0);

r=yvector;

BasalFactor=3;

PrismaticFactor=1;

%type pyramidal

Pyramidal_a_Factor=3;

%type I pyramidal

Pyramidal_ab_I_Factor=7;

sSBasal = slipSystem.basal(CS{2},3);

sSBasal = sSBasal.symmetrise('antipodal');

sSPrismatic = slipSystem.prismatic2A(CS{2},1);

sSPrismatic = sSPrismatic.symmetrise('antipodal');

% type II pyramidal

sSPyramidal = slipSystem.pyramidal2CA(CS{2},7);

sSPyramidal = sSPyramidal.symmetrise('antipodal');

rCS = ebsd_ci('Titanium(Alpha)').orientations \ r;

% Basal slip calculations

sFBasal = abs(sSBasal.SchmidFactor(rCS));

plot(ebsd('Titanium (Alpha)'),max(sFBasal,[],2))

But I see the following error:

Error using phaseList/checkSinglePhase (line 287)

There are no indexed data in this variable!

Error in phaseList/get.CS (line 154)

      id = checkSinglePhase(pL);

Error in EBSD/get.orientations (line 148)

      ori = orientation(ebsd.rotations,ebsd.CS);

Error in EBSD/subsref (line 40)

  [varargout{1:nargout}] = builtin('subsref',ebsd,s);

I couldn't fix it. Can anyone help. 

[ruediger Kilian]

What is the output of

ebsd_ci('Titanium(Alpha)').orientations

?

Cheers,
Rüdiger

[Karthik Palaniappan]

I get the following error for  " ebsd_ci('Titanium(Alpha)').orientations"

Error using phaseList/checkSinglePhase (line 287)

There are no indexed data in this variable!

Error in phaseList/get.CS (line 154)

      id = checkSinglePhase(pL);

Error in EBSD/get.orientations (line 148)

      ori = orientation(ebsd.rotations,ebsd.CS);

Error in EBSD/subsref (line 40)

  [varargout{1:nargout}] = builtin('subsref',ebsd,s);

[ruediger Kilian]

So this tells you that there's nothing in ebsd_ci('Titanium(Alpha)').orientations.

If you were still using ebsd_ci = ebsd(ebsd.ci>0.1) - taken from your previous post, maybe it should be ebsd.prop.ci - or maybe 0.1 is too large? What's does ebsd return? What is sum( ebsd.prop.ci>0.1)?

Cheers,
Rüdiger

[Karthik Palaniappan]

When I replaced  ebsd_ci = ebsd(ebsd.ci>0.1) to  ebsd_ci = ebsd(ebsd.prop.ci>0.1) I see the same error.  And 

sum( ebsd.prop.ci>0.1)  is 2935998

With regards,

Karthik 

[ruediger Kilian]

Ok, and what is the output of:

ebsd(ebsd.prop.ci>0.1) ?

[Karthik Palaniappan]

ans = EBSD (show methods, plot)

 

 Phase   Orientations           Mineral        Color  Symmetry  Crystal reference frame

     1  2652743 (90%)  Titanium (Alpha)   light blue       622       X||a, Y||b*, Z||c*

     2  283255 (9.6%)   Titanium (Beta)  light green       432                         

 

 Properties: ci, fit, iq, sem_signal, unknown1, unknown2, unknown3, unknown4, x, y

 Scan unit : um

With regards,

Karthik

[ruediger Kilian]

Ah, ok. Try to give the second titanium phase in your CS definition when you read in your data a different name, such as ... 'beta_titanium' .

> --
> 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.
> Visit this group at https://groups.google.com/group/mtexmail.
> For more options, visit https://groups.google.com/d/optout.

[Karthik Palaniappan]

I replaced the CS 'Titanium(Beta) ' with beta_titanium but the error exists. So I replaced Titanium(Alpha) with '1'. I didn't see the error. 

But, I see another error for the command  "plot(ebsd('1'),max(sFBasal,[],2))" 

Error using EBSD/plot (line 77)

The number of values should match the number of ebsd data!

With regards,

[ruediger Kilian]

Maybe posting your entire script might help.
Cheers,
Rüdiger

[Karthik Palaniappan]

I have attached the .ang file and the matlab code. Actually, I also need to find the Schmid factor map for Simple shear load. In that case how to include this in the code.  

 ald.ang

Karthik Palaniappan,

Sender notified by Mailtrack 09/21/18, 9:59:12 PM

[ruediger Kilian]

Replacing ebsd('1') with ebsd_ci('1') when plotting the Schmid factors should actually work.

Cheers,
Rüdiger

[Karthik Palaniappan]

Thank you very much for your help. The code worked successfully.

But I need to replace the r vector with a simple shear tensor to calculate the Schmid factor map. How can I make the necessary changes?

With regards, 

Karthik