Schmid Factor over Pole Figure in uniaxial tension

[maxencebuttard]

Hello everyone,

I am currently trying to plot Schmid Factor for a specific stress condition sigma = -stressTensor.uniaxial(vector3d.Y) over a pole figure. I know that it is possible to plot a Schmid Factor spherical function for different tension direction but I would prefer to plot a Schmid Factor for different orientations for a specific stress condition.

For now here is my code

% Define first slip system
n1 = Miller(1,1,1,ebsd('Cubic').CS,'hkl');
d1 = Miller(1,-1,0,ebsd('Cubic').CS,'uvw');

sS1 = slipSystem(d1,n1);
sSAll1 = sS1.symmetrise('antipodal');
sSLocal = ebsd('Cubic').orientations*sSAll1;

% Calculate Schmid factors for each orientation
SF= sSLocal .SchmidFactor(sigma);
[SFMax1, id1] = max(abs(SF), [], 2);

But I don't know how to plot it superimposed to a 100 pole figure.

[ralfHielscher]

How about

plotPDF(ebsd('Cubic').orientations, SFMax1, Miller(1,0,0,ebsd('Cubic').CS))

Ralf.