| 查看: 1741 | 回复: 2 | ||
yudonghua铁虫 (小有名气)
|
[求助]
求助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'); |
回复此楼» 猜你喜欢
FileCode能看出啥？
已经有15人回复
-
filecode
已经有8人回复
-
确定了，国自然21号放榜
已经有7人回复
-
帮忙看看fileCode
已经有5人回复
-
综述论文作为代表作会不会影响评审专家的印象分？
已经有11人回复
-
关于代码变化问题，想知道的进来
已经有24人回复
-
为什么网上很多人说本周 12号出结果
已经有3人回复
-
静等基金结果
已经有16人回复
-
我的国基提前知道中了，可是同事的操作让我实在接受不了，怎么会有这样的人
已经有11人回复
-
关于Filecode分析方法
已经有3人回复
zzuzdd
铁杆木虫 (著名写手)
- 应助: 2 (幼儿园)
- 金币: 13235.7
- 散金: 211
- 红花: 1
- 帖子: 1745
- 在线: 213.6小时
- 虫号: 1478332
- 注册: 2011-11-05
- 专业: 金属结构材料
2楼2018-06-09 23:42:13
yudonghua
铁虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 3.9
- 散金: 10
- 帖子: 52
- 在线: 53.8小时
- 虫号: 3273570
- 注册: 2014-06-14
- 性别: GG
- 专业: 金属材料
3楼2018-06-10 11:35:13
