I would like to plot Schmid factor maps for various slip systems of a single family. Removing "symmetrise" from the schmid factor statement would work?I would be happy if you could share the information on theSymmetrise function,what exactly it does.

[Bharat Reddy]

[tauMax,mActive,nActive,tau,ind] = calcShearStress(sigmaCS,n,b,'symmetrise');

For plotting Schmid factor maps of slip systems of a family, example 3 prismatic slip systems of hcp unit cell., removing of symmetrise function from the above statement would work?

Also I would like to know what exactly symmetrise function  does?

Thanks in advance

[grandr...@gmail.com]

No it wouldn't, the symmetrise function ensures the system calculates across all the crystallographically identical systems.  Of course if it were me I'd use the code I posted on this site a few days ago- change the name from magnesium to your HCP material and run the commands down to the section beginning 'to Plot' near the bottom (line 280 ish).  Then plot your schmid factor map using the below command.  Of course this code is overkill for you because it does all slip systems but you can just ignore that data. Check that the ones defined for Pyramidal are correct for you.


%to see the prismatic schmid factors
plot(ebsd('Magnesium'),PrismaticArray(:, 4))

Jessica

[Bharat Reddy]

Thanks for the reply.

The code was useful.

If I remove 'symmetrise' function from [tauMax,mActive,nActive,tau,ind] = calcShearStress(sigmaCS,n,b,'symmetrise'); and just write [tauMax,mActive,nActive,tau,ind] = calcShearStress(sigmaCS,n,b); , what values of schmid factor are plotted? Does it plot for individual slip system of a family?

Thanks in advance

[grandr...@gmail.com]

I suspect that what happens is that you only plot the individual slip system, in which case tauMax won't really be the maximum but just the value for that system.  To get the maximum you want it to examine all slip systems and compare to determine the max. 

Jessica

[ruediger Kilian]

Hi Bharat,
if you do not supply ‘symmetrise', it will only consider the plane m and direction n, no symmetry involved.
Have a look at:
>help calcShearStress
...
% Options
% symmetrise - consider also all symmetrically equivalent planes and directions
...
Try it out with/without and compare the results.

Cheers,
Rüdiger

> --
> 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.

[Bharat Reddy]

Thank you for the kind replies.

I actually want to plot the schmid factor maps for individual slip systems. Example:- (10-10)[1-210] slip system alone.

The problem statement is that I would like to check how active are individual slip systems rather than seeing at a maximum schmid factor for a family.

If I use the option 'symmetrise' for any set of available systems of a family example: {10-10}<1-210> , I am getting the same schmid factor map because symmetrise function symmetrises other independent slip systems of a single family which in turn gives the same result. 

 

On Tuesday, July 26, 2016 at 10:59:36 AM UTC+5:30, Bharat Reddy wrote:

[ruediger Kilian]

Hi,
so, if you just want to calculate the Schmid factor for a one slip system, just omit the ‘symmetrise’ and run it for each single plane/direction combination or have a look at the example below.
I’m not entirely sure what you mean by “how active individual slip systems are” but I guess you want to compare e,g, the Schmid factors for all members in certain slip systems.
Just have a look here:
https://mtex-toolbox.github.io/files/doc/PlasticDeformation.html#13
That’s from the examples:
CS = crystalSymmetry('432')
% one slip system
m = Miller(1,1,1,CS)
n = Miller(0,-1,1,CS)

% all slip systems
[mSym,l] = symmetrise(m,'antipodal');
[nSym,l] = symmetrise(n,'antipodal');
[row,col] = find(isnull(dot_outer(vector3d(mSym),vector3d(nSym))));
mSym = mSym(row)
nSym = nSym(col)

% some tensor
M = [[1 0 0];[0 0 0];[0 0 0]];
sigma = tensor(M,'name','stress')

% calc the Schmidtensor for each slip system
RSym = SchmidTensor(mSym,nSym)
% calculate the schmid factor for each one
tau = double(EinsteinSum(RSym,[-1,-2],sigma,[-1,-2],'name','Schmid factor'))'

The indices of tau correspond to those in nSym, mSym, so you have for each value also the corresponding slip system.

Cheers,
Rüdiger

[Bharat Reddy]

Thank you Kilian.

And yes I want to compare the Schmid factors of slip systems in a family.

I think this method should help me.

Thanks again for looking at this problem.

On Tuesday, July 26, 2016 at 10:59:36 AM UTC+5:30, Bharat Reddy wrote: