EBSD 实操:空间统计
错取向统计再往前一步:MODF、ODDF 与 OCF 怎样分离织构和邻接
从随机配对的概率推导出织构基线,再用同一权重比较真实邻接,读懂小角度峰究竟来自哪里。
先记住这三点
- 随机取向的基线与给定织构下随机邻接的基线是两种分布。
- OCF 比较相同权重和相同分箱中的实际频率与期望频率。
- 相邻像素数量很大,并不意味着拥有同样多独立晶界。
先把统计对象固定下来
设 R 把晶体笛卡尔坐标变到试样坐标,两个晶粒的相对取向取 Δ=R1ᵀR2。若换成 R2ᵀR1,轴的正负和旋转方向随之改变;角度分布可能不变,带方向的相关分析却不能混用。每对取向均需按两侧晶体对称取等价代表。
MODF 或 MDF 描述实际相邻区域的错取向分布。ODDF 在本章专指忽略空间位置、按给定总体取向分布随机配对得到的分布。它回答“若只保留这些取向,各晶粒随机遇见,会出现什么错取向”,而不是把所有角度看成同样可能。
从独立配对推导 ODDF
令 f(R) 对归一化 Haar 测度 dμ(R) 的积分为 1。先取第一颗晶粒 R,再独立取第二颗 R2。固定 Δ 时 R2=RΔ,旋转不改变 Haar 测度,所以随机配对的错取向密度为 D(Δ)=∫f(R)f(RΔ)dμ(R)。这就是取向分布在旋转群上的自相关。
检验归一化:∫D(Δ)dμ(Δ)=∫f(R)[∫f(RΔ)dμ(Δ)]dμ(R)=1。做晶体对称约化后,需要将所有等价代表的概率质量归入同一箱,而不是只计算一个矩阵并丢掉等价质量。即使 f=1,约化后的最小错取向角直方图仍不平坦。
实际邻接的分布怎样得到
将真实界面分为小段 j,长度 lj,错取向 Δj。某分箱 B 的二维长度概率为 M(B)=Σj lj·1[Δj∈B]/Σj lj。三维数据改用面元面积 Aj。若问题是“一颗晶粒平均有何种邻居”,可以按邻接边条数统计,但得到的是另一种概率。
晶粒面积加权的 f 产生体积意义上的随机配对。若 M 按界面长度加权,与面积权重的 D 直接相除可能混入晶粒尺寸和形状的影响。较稳妥的空间检验是固定原有邻接图和界面权重,随机置换晶粒取向,再计算随机基线;若尺寸与取向强相关,置换还需在合理尺寸组内进行。这里讲的是统计设计,不要求读者为理解原理运行模拟。
两组分手算:峰高并不自动表示特殊邻接
设 A、B 两类取向各占一半;把一对的次序合并,随机配对给 P(AA)=1/4、P(BB)=1/4、P(AB)=1/2。若 A、B 都很尖锐,ODDF 必然在近零角和 A–B 特征错取向处出现峰,即使空间完全随机。
假设教学用的 100 条等权界面中,AA 有 10 条、BB 有 10 条、AB 有 80 条,那么 AB 的观测概率为 0.80,随机期望为 0.50,富集比为 1.60;同类总概率是 0.20/0.50=0.40。这里的数字只用于跟算,并非论文实验结果。
OCF 的分母、零值与误差
定义 OCF(B)=M(B)/D(B)。大于 1 表示比该织构的随机邻接更常见,小于 1 表示被抑制。它是相对富集,不是面积分数,也不是重新归一化的普通概率密度;它满足的是加权关系 ΣB OCF(B)D(B)=1。
当 D(B) 很小,少量观测就能产生极高比值;当 D(B)=0,该箱不可直接除。应标为缺少基线支持、合并有物理意义的邻近箱,或说明采用的平滑模型。不能把无穷大、缺失值和零相关都涂成同一种颜色。
为什么角度和轴要一起看
同样 60° 可以绕完全不同的轴旋转,只有角峰不能确认 Σ3。应先看完整轴角分布,再将特定轴附近的错取向投影到角度直方图。如果选择“最接近〈110〉的等价旋转”,随机基线也必须使用同样选择规则。
Rodrigues 向量 r=a tan(θ/2) 把定轴旋转放到一条直线上,适合识别特征关系;但直接在 r 空间做欧氏等体积分箱并不等于在旋转群上等概率分箱。绘图坐标方便,不会消除测度和对称的影响。
真实案例:高纯铝电容箔的连续再结晶
Engler 与 Huh(1999)将组织和织构共同用于研究高纯铝箔:早期退火可由扩展回复形成细小组织并保留轧制织构,后续立方取向通过不连续晶粒长大变得突出。第3版书印刷页383–384进一步展示了该问题中实测 MODF、织构生成的 ODDF 和二者之比。
读者应按三步读图:先识别总体织构使小角配对本来就有多常见,再观察真实相邻关系的额外富集,最后结合退火过程和组织变化讨论机制。单张 OCF 峰不能独立证明所有材料都发生连续再结晶;它是整条证据链的一部分。
拿到自己的结果时怎样写分析
先写相、取向约定、晶粒/亚晶粒定义,再写相邻关系和长度或面积权重;给 MODF 与匹配权重的随机基线,最后给 OCF。将峰对应的界面返回地图,确认它们来自许多独立区域,而非一条很长界面的重复像素。
自检题:若 A 占90%、B占10%,随机异类配对是多少?答案是2×0.9×0.1=0.18。若实际异类界面占30%,富集比约1.67,绝不能因为30%小于50%就称其被抑制。
来源与延伸阅读
以下链接用于核对概念、标准状态或原始文献。本页为中文教学整理,不复制受版权保护的标准正文。