24小时热门版块排行榜    

查看: 9765  |  回复: 12
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

fanarsenal

木虫 (正式写手)

[交流] 从欧拉角(Euler angle)计算取向差(misorientation) 已有7人参与

最近要根据ebsd计算晶粒取向差,搜索论坛之后发现不少虫友也有类似问题,但大家的解答的都不那么直接,特地发帖讨论一下。注意:1)此处只是讨论已知两组欧拉角之后,如何计算晶粒的取向差,关于如何从ebsd处理软件得到对应晶粒的欧拉角,请自行查找相关资料;2)此处仅讨论立方系。

第一步:分别将两组欧拉角转化成取向矩阵g1 和g2
立方系欧拉角(φ1,Φ,φ2)和取向矩阵g 转换关系如下[1]:
[ cos φ1 cos φ2   - sin φ1 sin φ2  cos  Φ                       sin  φ1 cos φ2  +cos φ1 sin φ2 cos Φ                            sin φ2 sin Φ  
-cos φ1 sin φ2  - sin φ1 cos φ2 cos Φ                         -sin φ1 sin φ2  + cos φ1 cos φ2 cos Φ                           cos φ2 sin Φ
           sin φ1  sin Φ                                                                    -cos φ1  sin Φ                                                            cos Φ        ]

第二步,计算旋转矩阵 r12, 即g2乘以g1的逆矩阵。
R12 = g2 inv(g1) .

第三步,对矩阵进行线性变换,使得其对角元素为其所在行和列的最大值 [2],记矩阵元素为r_ij.
立方晶系对称操作太多,对应同一个旋转矩阵,有1000多种旋转可能。计算取向差一般都是角度最小的那个旋转角度及其对应的旋转轴。另外,如果将对角元素从11到33按照从大到小排列,则旋转轴的三个指数也是从大到小排列。

第四步 根据第三步得到的旋转矩阵,得出旋转角和旋转轴
theta = arccos(0.5*(r11+r22+r33-1));
c = [c1, c2, c3] = 1/(2 sin theta ) [ r32-r23, r13-r31, r21-r12];
此处c的分量都是小数。

[1] 杨平,电子背散射衍射技术及其应用.
[2] w. bollmann, crystal lattices, interfaces, matrices: an extension of crystallography.

下面给个例子
两个晶粒 欧拉角分别是[281.7 33.3 2.7] 和 [138.7 45.6 38.6]
g1 = [    0.2411   -0.9702    0.0259
    0.8080    0.2154    0.5484
   -0.5376   -0.1113    0.8358]
g2 = [   -0.8752    0.1879    0.4457
    0.1078   -0.8226    0.5584
0.4716    0.5368    0.6997]

r12 =[ -0.3818   -0.4222    0.8222
    0.8384    0.2161    0.5003
   -0.3889    0.8803    0.2715]
第一行和第三行交换,第一列和第二列交换
r12 = [  0.8803 -0.3889     0.2715
      0.2161  0.8384    0.5003
-0.4222  -0.3818   0.8222]
所以 θ =39.6057
c = [ -0.6919    0.5441    0.4745] 转化成整数大概是[-4 3 3]

Edit:最开始用图来表示,发现大部分公式都挂了,还是蹩脚的打上去吧。
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

with you

木虫 (著名写手)


小木虫: 金币+0.5, 给个红包,谢谢回帖
请问楼主,线性变换使其对角线元素为其所在行和列最大值,是将所有元素取绝对值进行比较大小吗,如果上述r12矩阵的第一行变成了 [-0.8803 -0.3889     0.2715], 那么最大值是 第1个还是第3个?
9楼2023-02-06 10:57:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 13 个回答

凌云qq

铁杆木虫 (职业作家)


小木虫: 金币+0.5, 给个红包,谢谢回帖
谢谢楼主,虽然现在不用,后面可能用得到

发自小木虫Android客户端
2楼2018-12-04 07:08:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

with you

木虫 (著名写手)


小木虫: 金币+0.5, 给个红包,谢谢回帖
楼主, 您好。请问如果我知道两个相邻晶粒中所有原子的坐标,怎么计算取向差呢?还请您指点一二
3楼2019-06-21 21:25:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fanarsenal

木虫 (正式写手)

引用回帖:
3楼: Originally posted by with you at 2019-06-21 21:25:07
楼主, 您好。请问如果我知道两个相邻晶粒中所有原子的坐标,怎么计算取向差呢?还请您指点一二

你知道所有原子坐标,直接用立体几何不就好了?
4楼2019-12-07 06:50:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[博后之家] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +4 ASdOkHsho7FD 2026-08-28 5/250 2026-08-29 14:24 by jCd0dEvKHShX
[考博] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 ASdOkHsho7FD 2026-08-28 7/350 2026-08-29 14:15 by jCd0dEvKHShX
[论文投稿] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 ASdOkHsho7FD 2026-08-28 4/200 2026-08-29 14:13 by jCd0dEvKHShX
[基金申请] 麻烦专家们看看评委们的意见(F口面上) +5 gdd2018 2026-08-28 10/500 2026-08-29 14:10 by Jacob678
[找工作] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +4 ASdOkHsho7FD 2026-08-28 6/300 2026-08-29 11:29 by jCd0dEvKHShX
[考研] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +6 ASdOkHsho7FD 2026-08-28 9/450 2026-08-29 11:18 by jCd0dEvKHShX
[基金申请] 中青基了要发朋友圈吗? +3 349506619 2026-08-28 3/150 2026-08-29 10:49 by 木风9012
[基金申请] 国自然评审意见 +13 wangmingqi 2026-08-28 19/950 2026-08-29 10:22 by Poppy1104
[基金申请] 2026年叶企孙基金 +4 bud_bud 2026-08-27 7/350 2026-08-29 07:23 by foolishmani
[基金申请] 系统查不到 +11 董八千 2026-08-26 11/550 2026-08-28 18:06 by Leogzhya
[基金申请] 怎么看青基中了没有啊 +5 叶九微 2026-08-26 5/250 2026-08-27 10:35 by l_zh2008
[基金申请] 我不理解! +15 Edward_pc 2026-08-26 23/1150 2026-08-26 20:34 by zzuzxg
[基金申请] 2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战 +7 process2012 2026-08-23 10/500 2026-08-26 19:23 by hmhminy
[基金申请] 2026年8月25日国自然放榜前突然收到列入评审专家邮件,有关系吗? +25 木水思豆 2026-08-25 28/1400 2026-08-26 14:53 by draco1987
[基金申请] 出来了 +9 trojank 2026-08-26 9/450 2026-08-26 14:25 by 宝贝虫子
[基金申请] 国际合作可查了,中了面上 (EPI+1)(金币+50) +18 Ldrop2023 2026-08-26 18/900 2026-08-26 11:15 by cmrandy
[基金申请] 系统进不去 +4 yanglien 2026-08-26 5/250 2026-08-26 11:10 by wenfengw83
[基金申请] 今天务委会开完了,明天出结果吗 +19 angus9576 2026-08-25 23/1150 2026-08-26 10:03 by zp519
[基金申请] 牛来!米来!面来! +8 beefly 2026-08-26 8/400 2026-08-26 08:37 by xuzhipiao
[基金申请] 某些机构,以效率低为荣,以效率低作为存在感 +9 yuleib84 2026-08-25 10/500 2026-08-25 17:14 by alexon
信息提示
请填处理意见