| 查看: 1749 | 回复: 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'); |
回复此楼» 猜你喜欢
澳大利亚 Murdoch University 全奖博士招生（3个名额）地质化工冶金领域
已经有66人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有30人回复
-
我的国基提前知道中了，可是同事的操作让我实在接受不了，怎么会有这样的人
已经有15人回复
-
关于Filecode分析方法
已经有12人回复
-
咱们一起用铁证分析2026国家社科基金中标与否
已经有11人回复
-
小木虫上这么多卖论文的，真有人买论文么？感觉没必要啊
已经有10人回复
-
阿姨
已经有4人回复
-
一作与独作在应聘高校教师时区别大吗
已经有4人回复
-
静等基金结果
已经有21人回复
-
FileCode能看出啥？
已经有28人回复
yudonghua
铁虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 3.9
- 散金: 10
- 帖子: 52
- 在线: 53.8小时
- 虫号: 3273570
- 注册: 2014-06-14
- 性别: GG
- 专业: 金属材料
3楼2018-06-10 11:35:13
zzuzdd
铁杆木虫 (著名写手)
- 应助: 2 (幼儿园)
- 金币: 13235.7
- 散金: 211
- 红花: 1
- 帖子: 1745
- 在线: 213.6小时
- 虫号: 1478332
- 注册: 2011-11-05
- 专业: 金属结构材料
2楼2018-06-09 23:42:13
