24小时热门版块排行榜

返回列表

yudonghua

铁虫 (小有名气)

应助: 0 (幼儿园)
金币: 3.9
散金: 10
帖子: 52
在线: 53.8小时
虫号: 3273570
注册: 2014-06-14
性别: GG
专业: 金属材料

[求助] 求助EBSD计算Schmid因子

求教，我做镁合金EBSD，实验软件，应该是牛津的，给出的基面滑移Schmid因子，和我自己参考论文和杨平的书上的计算结果完全不一样。下面是我写的matlab程序，请问问题出在哪里？
Eoule=xlsread('E:Eoule2.xls','sheet1'); %读取欧拉角
A1=sum(abs(Eoule'));
index=find(A1==0);
Eoule(index,

=[];%把标定率为0的去除
p1=Eoule(:,1);
p=Eoule(:,2);
p2=Eoule(:,3);
%%
h1=0;  k1=0;  i1=0;  l1=1;  %设置需要计算的滑移系
u1=1;  v1=1;  t1=-2;  w1=0;
dhkl=sqrt(2*(h1+k1)^2+3*k1^2+3*l1^2*(1/1.624)^2);  %滑移系正交归一化
duvw=sqrt(3*(u1+v1/2)^2+9*v1^2/4+w1^2*(1.624)^2);
hn=2*(h1+k1)/dhkl;
kn=sqrt(3)*k1/dhkl;
ln=sqrt(3)*l1*(1/1.624)/dhkl;
un=(2*u1+v1)*sqrt(3)/2/duvw;
vn=3*v1/2/duvw;
wn=w1*1.624/duvw;
%%
h=sin(p2).*sin(p); %将欧拉角转化为Mille指数
k=cos(p2).*sin(p);
l=cos(p);
for i=1:length(p1)
m1(i,1)=(un*h(i)+vn*k(i)+wn*l(i))/(sqrt(h(i)^2+k(i)^2+l(i)^2)*sqrt(un^2+vn^2+wn^2)); %计算schmid因子
m2(i,1)=(hn*h(i)+kn*k(i)+ln*l(i))/(sqrt(h(i)^2+k(i)^2+l(i)^2)*sqrt(hn^2+kn^2+ln^2));
m(i,1)=abs(-m1(i,1)*m2(i,1));
end
[n,xout]=hist(m,50);
bar(xout,n/length(m),'b');

回复此楼