-
Notifications
You must be signed in to change notification settings - Fork 0
/
SS_setFCC_partial.m
43 lines (39 loc) · 1.85 KB
/
SS_setFCC_partial.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
CS = crystalSymmetry('-43m');
SlipSystem = struct('n',zeros(1,3),'b',zeros(1,3));
SlipSystem(1).n=[1,1,1]; SlipSystem(1).b =[0,-1,1];
SlipSystem(2).n=[1,1,1]; SlipSystem(2).b =[-1,0,1];
SlipSystem(3).n=[1,1,1]; SlipSystem(3).b =[-1,1,0];
SlipSystem(4).n=[-1,1,1]; SlipSystem(4).b =[1,1,0];
SlipSystem(5).n=[-1,1,1];SlipSystem(5).b =[1,0,1];
SlipSystem(6).n=[-1,1,1];SlipSystem(6).b =[0,-1,1];
SlipSystem(7).n=[1,-1,1]; SlipSystem(7).b =[1,1,0];
SlipSystem(8).n=[1,-1,1]; SlipSystem(8).b =[-1,0,1];
SlipSystem(9).n=[1,-1,1]; SlipSystem(9).b =[0,1,1];
SlipSystem(10).n=[1,1,-1];SlipSystem(10).b =[-1,1,0];
SlipSystem(11).n=[1,1,-1]; SlipSystem(11).b=[1,0,1];
SlipSystem(12).n=[1,1,-1]; SlipSystem(12).b =[0,1,1];
SlipSystem(13).n=[1,1,1]; SlipSystem(13).b =[0,1,-1];
SlipSystem(14).n=[1,1,1]; SlipSystem(14).b =[1,0,-1];
SlipSystem(15).n=[1,1,1]; SlipSystem(15).b =[1,-1,0];
SlipSystem(16).n=[-1,1,1]; SlipSystem(16).b =[-1,-1,0];
SlipSystem(17).n=[-1,1,1];SlipSystem(17).b =[-1,0,-1];
SlipSystem(18).n=[-1,1,1];SlipSystem(18).b =[0,1,-1];
SlipSystem(19).n=[1,-1,1]; SlipSystem(19).b =[-1,-1,0];
SlipSystem(20).n=[1,-1,1]; SlipSystem(20).b =[1,0,-1];
SlipSystem(21).n=[1,-1,1]; SlipSystem(21).b =[0,-1,-1];
SlipSystem(22).n=[1,1,-1];SlipSystem(22).b =[1,-1,0];
SlipSystem(23).n=[1,1,-1]; SlipSystem(23).b=[-1,0,-1];
SlipSystem(24).n=[1,1,-1]; SlipSystem(24).b =[0,-1,-1];
for si = 1:1:24;
n = SlipSystem(si).n;
b = SlipSystem(si).b;
hkl = Miller(n(1),n(2),n(3),CS);
uvw = Miller(b(1),b(2),b(3),CS);
beta_ij = SchmidTensor(uvw,hkl,'generalized');
m_ij = SchmidTensor(uvw,hkl,'symmetric');
q_ij = SchmidTensor(uvw,hkl,'antisymmetric');
SlipSystem(si).m = m_ij;
SlipSystem(si).beta = beta_ij;
SlipSystem(si).q = q_ij;
end
save('FCC_SS24_Set.mat','SlipSystem');