24小时热门版块排行榜    

查看: 9769  |  回复: 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, 给个红包,谢谢回帖
楼主, 您好。请问如果我知道两个相邻晶粒中所有原子的坐标,怎么计算取向差呢?还请您指点一二
3楼2019-06-21 21:25:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 13 个回答

凌云qq

铁杆木虫 (职业作家)


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

发自小木虫Android客户端
2楼2018-12-04 07:08:39
已阅   回复此楼   关注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的回帖

gogogo刘

新虫 (初入文坛)


小木虫: 金币+0.5, 给个红包,谢谢回帖
楼主您好,关于第三步线性变换不是很理解,这个参考文献2网上找不到,您可以分享给我吗。感谢楼主
5楼2021-04-07 09:31:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[论文投稿] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 G6APbkg8SA6w 2026-08-29 3/150 2026-08-30 02:33 by ZPa0EcMwuECS
[考博] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 G6APbkg8SA6w 2026-08-29 3/150 2026-08-30 02:22 by ZPa0EcMwuECS
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 G6APbkg8SA6w 2026-08-29 3/150 2026-08-30 02:18 by ZPa0EcMwuECS
[公派出国] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 gy1nBQXYQJqL 2026-08-29 4/200 2026-08-30 01:45 by ZPa0EcMwuECS
[找工作] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 gy1nBQXYQJqL 2026-08-29 6/300 2026-08-30 01:36 by ZPa0EcMwuECS
[博后之家] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +4 gy1nBQXYQJqL 2026-08-29 6/300 2026-08-30 01:36 by ZPa0EcMwuECS
[论文投稿] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 gy1nBQXYQJqL 2026-08-29 6/300 2026-08-30 01:25 by ZPa0EcMwuECS
[考博] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +6 ASdOkHsho7FD 2026-08-28 9/450 2026-08-30 01:04 by ZPa0EcMwuECS
[基金申请] 中青基了要发朋友圈吗? +4 349506619 2026-08-28 4/200 2026-08-29 22:41 by alongwaytogo
[基金申请] 国自然面上复盘~欢迎讨论 (金币+15) +15 晴天加油 2026-08-26 16/800 2026-08-29 18:28 by symmetry
[基金申请] 有没有仍没收到信息的 +5 德尚中行 2026-08-27 6/300 2026-08-29 18:10 by manplx
[基金申请] 为什么到现在没收到通知? +4 tannykie 2026-08-29 4/200 2026-08-29 17:14 by lmz0216
[基金申请] 我就是申请一个面上项目而已,这评审意见是按照杰青的条件评的吧? +6 gouxfjh 2026-08-28 9/450 2026-08-29 09:59 by jklily
[教师之家] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 ASdOkHsho7FD 2026-08-28 7/350 2026-08-29 09:26 by G6APbkg8SA6w
[基金申请] 为什么资助数各大高校都创新高,自己申请怎么就这么难 +12 Kittylucky 2026-08-27 13/650 2026-08-29 00:04 by superceng
[基金申请] 怎么查啊 +6 huang1991js 2026-08-26 6/300 2026-08-28 08:42 by winsaint
[基金申请] 基金未中,这种答复是模板吗? +5 zhaosm1982 2026-08-27 6/300 2026-08-27 16:00 by lfy8008
[基金申请] 国合现在查不到了吗? +10 chengyan1220 2026-08-24 20/1000 2026-08-26 08:57 by peasantsprig
[基金申请] 牛来!米来!面来! +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
信息提示
请填处理意见