the identification of variant id

[SEEHAIHAI]

As Niessen mentioned in his famous paper 'Parent grain reconstruction from partially or fully transformed microstructures in
MTEX', the variant ID in Mtex was consistent with Morito's paper 'The morphology and crystallography of lath martensite in
Fe-C alloys'. When i checked this with ortools, the misorientation between variants is diffrent from Table 2 in Morito's paper.
The code is as follows:
csFCC = crystalSymmetry('432',[3 3 3],'mineral','FCC')
csBCT = crystalSymmetry('432',[3.3 3.3 3.3],'mineral','BCT')
oriFCC = orientation.cube(csFCC)
R = orientation.KurdjumovSachs(csBCT,csFCC)%BCT-->FCC
oriFCC_variant = variants(R,oriFCC)
plotPDF(oriFCC_variant,Miller(0,0,1,csBCT,'hkl'))
for i = 1:24
misori(i) = inv(oriFCC_variant(1))oriFCC_variant(i);
misaxis = [misori(i).axis.h,misori(i).axis.k,misori(i).axis.l]/norm([misori(i).axis.h,misori(i).axis.k,misori(i).axis.l]);
misangle(i) = mori(i).angle180/pi;
disp(['v1 to v', num2str(i), ' | axis: ', num2str(misaxis), ' | angle: ', num2str(misangle(i)), '°'])
end

[AzdiarGazder]

Hi @SEEHAIHAI,

Apologies for the delay in replying.
I think there's a misunderstanding of the commands that enable you to go about parent grain reconstruction and analysis.
Please have a detailed read of the sub-section labelled "Phase transitions" in the MTEX documentation.

In the meantime, here is the correct script to generate what you are after.

clear all; close all; clc

% define the parent symmetry
cs_fcc = crystalSymmetry('m-3m', [3.6599 3.6599 3.6599], 'mineral', 'fcc');
% define the child symmetry
cs_bcc = crystalSymmetry('m-3m', [2.866 2.866 2.866], 'mineral', 'bcc');

% define the Kurdjumov-Sachs parent fcc to child bcc orientation relationship
p2c = orientation.KurdjumovSachs(cs_fcc,cs_bcc);

% calculate the variants
fcc2bcc = p2c.variants;
bcc2bcc = p2c.variants.*inv(p2c.variants(1));
variantAngle = angle(bcc2bcc);
variantAxis = axis(bcc2bcc,cs_bcc);

for ii = 1:24   
    disp(['V1 to V', num2str(ii), ' | axis: ', num2str(round(variantAxis(ii).uvw,4)), ' | angle: ', num2str(variantAngle(ii)./degree), '°'])
end

% define a symmetric fcc parent orientation
ori_fcc = orientation.cube(cs_fcc);
ori_fcc = ori_fcc.symmetrise;
% compute all child bcc orientations related to the symmetrised fcc parent orientation
ori_bcc = ori_fcc * inv(p2c);

% plot the pole figure
cMap = jet(24);
for ii = 1:24
  plotPDF(ori_bcc(ii),Miller(0,0,1,cs_bcc),...
      'equal','antipodal',...
  'markerSize',10,...
  'markerFaceColor',cMap(ii,:),...
  'markerEdgeColor','k');
  hold all;
end
colormap(jet(24));
clim([1 24]);
colorbar('location','eastOutSide','lineWidth',1.25,'tickLength', 0.01,...
    'YTick', 1:24,...
    'YTickLabel',num2str((1:24)'), 'YLim', [1 24],...
    'TickLabelInterpreter','latex','FontSize',14,'FontWeight','bold');

**Lastly, and importantly, if an explanation from the maintainers resolves your issue, please press the "Mark as answer" button on the appropriate response. Your assistance in this regard enables us to close discussions as complete. **

Hope this helps.

Warm regards,
Azdi