combining data from 2 EBSD data sets of the same sample

[Élise Cossette]

Dear Ralf and MTEX users,

I am still using MTEX 3.3.1

Due to some bad luck, I have two different ebsd data sets from the same sample - one contains quartz orientations and the other contains muscovite orientations. I wonder if it is possible to combine the orientations in order to calculate the average material tensor of these two phases. For now I can only calculate the tensor for single phases, since the calcTensor command only allows one entry of type 'ebsd'.

Many thanks,

Élise

[Ralf Hielscher]

Hi Elise,

you may want to combine those ebsd data sets by

ebsd = [ebsd1,ebsd2]

I hope this helps,

Ralf.

[Élise Cossette]

Thank you Ralf,

Is it possible that the combination doesn't work since ebsd1 is from an OIM-TSL program and ebsd2 is from an Oxford Channel program? I get this error when trying to combine the two:

Reference to non-existent field 'ci'.

Error in EBSD/horzcat (line 16)

    ebsd.options.(cfn) = vertcat(ebsd.options.(cfn),varargin{i}.options.(cfn));

Élise

[Ralf Hielscher]

Hi Élise,

if the EBSD data are from different vendors then you need to make them compatible, i.e., you have to remove every property that is not present in both data sets.

As an example you can remove the field 'ci' MTEX is complaining about by

ebsd1 = set(ebsd1,'ci',[])

You should apply this syntax to any property that is not present in both data sets. Hope this helps,

Ralf.

On Thursday, December 4, 2014 9:27:16 PM UTC+1, Élise Cossette wrote:

[David Mainprice]

I would suggest the following.

A) Verify by plotting pole figures that the orientation of the data from ebsd1 (OIM-TSL ) quartz (?) and orientation of the data from ebsd2 (Oxford Channel ) muscovite (?)

are correctly aligned in specimen coordinates within MTEX

B) Define elasttic stiffness tensor of alpha-quartz and muscovite

You may need to adjust the tensor CS a,b,c, alpha,beta,gamma to agree with you CS of your imported data

BUT do not change the 'x||a','z||c’ or 'x||a','z||c*’ of the tensor CS as they are tensor reference axes.

%************************************************

% Define elastic stiffness tensor (GPa)

%************************************************

%

% Reference Elastic constants

% alpha-Quartz McSkimin et al 1965 J.Appl.Phys.v.36p.1624-1632

%

% Reference Crystal Structure

% alpha quartz Trigonal SG P3221 25°C d=2.648 g/cm3

%

% Define density (g/cm3)

 rho_quartz= 2.6480;

%

% Define Cartesian Tensor crystal symmetry & frame

cs_Tensor = symmetry('-3m',[  4.9134  4.9134  5.4052],...

 [  90.0000  90.0000 120.0000]*degree,'x||a','z||c',...

 'mineral','Alpha-Quartz 1965');

%

% elastic Cij stiffness tensor (GPa) as matrix M

 M =.... 

 [[   86.80    7.04   11.91  -18.04    0.00    0.00];...

 [     7.04   86.80   11.91   18.04    0.00    0.00];...

 [    11.91   11.91  105.75    0.00    0.00    0.00];...

 [   -18.04   18.04    0.00   58.20    0.00    0.00];...

 [     0.00    0.00    0.00    0.00   58.20  -18.04];...

 [     0.00    0.00    0.00    0.00  -18.04   39.88]];

%

% M as stiffness tensor C with MTEX tensor command

C_quartz = tensor(M,cs_Tensor);

%

%************************************************

% Define elastic stiffness tensor (GPa)

%************************************************

%

% Reference Elastic constants

% Muscovite Vaughan and Guggenheim JGR 1986

%

% Reference Crystal Structure

% Rothbauer, R., (1971) Untersuchung eines 2M1-Muskovits

% mit Neutronenstrahlen, Neues J. Mineral. Monal., 143-154.

% Muscovite C2/c x||a z||c*

% KAl2AlSi3O10(OH)2  2M1 

% Define density 2.834 (g/cm3)

rho_muscovite = 2.834;

%

% Define Cartesian Tensor crystal symmetry & frame

cs_Tensor = symmetry('2/m',[  5.1918 9.0153 20.0457],...

 [  90.00  95.74  90.00]*degree,'x||a','z||c*',...

 'mineral','Muscovite');

%

% elastic Cij stiffness tensor (GPa) as matrix M

 M =.... 

 [[  181.00   48.80   25.60    0.00  -14.20    0.00];...

 [    48.80  178.40   21.20    0.00    1.10    0.00];...

 [    25.60   21.20   58.60    0.00    1.00    0.00];...

 [     0.00    0.00    0.00   16.50    0.00   -5.20];...

 [   -14.20    1.10    1.00    0.00   19.50    0.00];...

 [     0.00    0.00    0.00   -5.20    0.00   72.00]];

%

% M as stiffness tensor C with MTEX tensor command

C_muscovite = tensor(M,cs_Tensor);

C)  Calculate averages for each mineral

[C_Voigt_quartz,C_Reuss_quartz,C_Hill_quartz] = ...

    calcTensor(ebsd1(‘Quartz'),C_quartz)

[C_Voigt_muscovite,C_Reuss_muscovite,C_Hill_muscovite] = ...

    calcTensor(ebsd2(‘Muscovite '),C_muscovite )

D) Define volume fractions of Quartz (vf_quartz) and Muscovite (vf_muscovite) from optical microscope or SEM and calculated weighted average

For example Hill average

C_Hill_Aggregate= vf_Quartz*C_Hill_quartz + vf_muscovite*C_Hill_muscovite

density_Aggregate = vf_Quartz*rho_quartz + vf_muscovite*rho_muscovite

all the best David

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

David Mainprice

15 Les Romarins

34270 Les Matelles

Tel. +33 4 67 67 03 30

(In France 04 67 67 03 30)

E-mail: davidma...@icloud.com