forked from tori-miller/Ti64-VPSC-CoRotation-Analysis
-
Notifications
You must be signed in to change notification settings - Fork 1
/
UnusedCodeBlocks.m
133 lines (97 loc) · 2.42 KB
/
UnusedCodeBlocks.m
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
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
From Corotationdataprocessing:
%% Data processing
% syntax: countourf(r,value)
% stress-strain data first:
% at yield
% figure
% contourf(r, stress(:,1))
% mtexColorMap white2black
% colorbar
%
% figure
% contourf(r,stress(:,end))
% mtexColorMap white2black
% colorbar
% figure(1)
% contourf(r,stress(:,1))
% mtexFig = gcm;
% contourf(r,stress(:,end),mtexFig.nextAxis)
% mtexColorMap white2black
% setColorRange('equal')
% mtexColorbar('multiple')
%
%%
% calc difference from oringinal angle.
adiff = abs(a - 1.79); %0.79 is the value at strain==0
adiff = adiff * (180/pi);
for j = 1:4:segments
% contour plots
figure
contourf(r,adiff(:,j))
colorbar
mtexColorMap white2black
end
%
%% Data processing: slip mode activity
% initial
figure
for k = 1:length(r(:))
hold on
plot(r(k),'MarkerFaceColor',(1/sum([act(k,1,1) act(k,1,2) act(k,1,3)]))*[act(k,1,1) act(k,1,2) act(k,1,3)],'MarkerSize',5);
end
%%
% final
figure
for k = 1:length(r(:))
hold on
plot(r(k),'complete','MarkerFaceColor',(1/sum([act(k,41,1) act(k,41,2) act(k,41,3)]))*[act(k,41,1) act(k,41,2) act(k,41,3)],'MarkerSize',5);
end
%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
From Final_Corotation_Processing
%% Plot initial orientations
% figure
% plotPDF(ori_a0,h_a)
%
% figure
% plotPDF(ori_b0,h_b,'MarkerFaceColor','r')
%%
for j = 2:4:segments
% contour plots
figure
contourf(r,adiv(:,j))
colorbar
mtexColorMap cool
figure
contourf(r,bdiv(:,j))
colorbar
mtexColorMap hot
end
%% Calculate and plot the deviation from the ideal alpha beta separation.
% calc difference from oringinal angle.
adiff = abs(a - 0.7900); %0.79 is the value at strain==0
adiff = adiff * (180/pi);
for j = 2:4:segments
% contour plots
figure
contourf(r,adiff(:,j))
colorbar
mtexColorMap white2black
end
%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
From MisorientationProcessing
%% Plot Alpha Difference From Original
figure
contourf(r,adiv(:,22))
colorbar
mtexColorMap blue2red
%% Plot Beta Difference From Original
figure
contourf(r,bdiv(:,22))
colorbar
mtexColorMap blue2red
%% Plot Alpha Misorientation From Original
figure
contourf(r,mangle(:,22))
colorbar
mtexColorMap blue2red