-
Notifications
You must be signed in to change notification settings - Fork 5
/
Mtex_commands.txt
38 lines (24 loc) · 1.46 KB
/
Mtex_commands.txt
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%% To Calculate Volume of fiber along certain axis %%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
ori=calcOrientations(odf, 10000);
h = Miller(0,0,1,CS,'uvw');
r = zvector;
100 * fibreVolume (ori, h, r, 15 * degree)
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%% To Calculate ODF sections plot %%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
plot(odf,'phi2',[45]*degree,'contourf')
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%% To Calculate ODF from Pole Figures (XRD data) %%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
pf(pf.intensities<0) = 0;
odf = calcODF(pf,'silent', 'resolution',5*degree);
%%%%%%%%%%%%%%%%%%%%%%%%%%%% ODF calculation from VPSC output %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
texfile1= 'TEX_PH1.OUT';
CS = crystalSymmetry('m-3m');
SS = specimenSymmetry('mmm');
data= loadOrientation_generic(texfile1,'CS',CS,'SS',SS,'ColumnNames',{'Euler1' 'Euler2' 'Euler3'}, 'Bunge');
odf = calcDensity(data);
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%