EBSD 实操:取向梯度

KAM 完整推导:邻域、阈值、步长与局部变形

手算一个邻域的 KAM,再推导连续取向梯度与噪声如何进入结果,说明为什么红色不能直接换成塑性应变。

阅读设置

先记住这三点

  1. KAM 是中心点与规定邻点之间错取向角的平均。
  2. 邻域物理大小、阈值、晶界排除与角噪声共同决定数值。
  3. KAM 可用于比较局部取向变化,但没有材料通用的应变或位错换算。
01

先定义两点之间的角

每点取向 Ri 将晶体坐标转到试样。两个同相像素的相对旋转可由 RjᵀRi 表示,再在晶体对称等价中取最小旋转角。对一个有效旋转 Δ,cosω=(trΔ−1)/2;数值实现需将余弦限制在 [−1,1],避免浮点误差使反余弦失效。

这里算的是完整三维旋转差,不是两个 IPF 颜色的色差,也不是 Euler 三个数逐项相减。不同相之间一般不能直接用同一最小错取向定义进入 KAM。

02

写清核与有效邻点

设中心 i 的几何邻域为 Ni。筛选后有效集合 Ei={j∈Ni:同相、有效索引、满足晶粒限制且 ωij<θcut}。则 KAMi=(Σj∈Ei wijωij)/(Σj∈Ei wij)。等权时 wij=1,分母就是有效邻点数。没有有效邻点时结果应记为缺失,而不是自动填 0。

有的软件按邻域阶数选点,有的按物理半径;方格第一圈与六角网格第一圈点数不同;是否包括对角点也会改变距离分布。报告“一级 KAM”而不说网格和算法,别人未必能复做。

03

跟着手算一个核

课堂示例的中心点与四个近邻错取向分别为 0.2°、0.4°、0.6°、7.0°。不筛选时平均为 2.05°。若 θcut=5°,最后一个点被排除,则 KAM=(0.2+0.4+0.6)/3=0.4°。分母是 3,不能仍除以原来的 4。

若 7.0° 来自另一晶粒,排除它可以避免把晶界跳变当粒内梯度;若它来自真实强烈局部转动,阈值也可能删去研究对象。阈值不是“去掉坏数据”的万能开关,必须结合晶界与原始质量判断。交互图采用这一组课堂数值,可改变阈值观察分母与结果。

跟着推导 · 教学几何

KAM:让有效邻点参与平均

四个邻点已用完整取向和晶体对称算好错取向角。这里练习 KAM 最容易算错的一步:排除邻点后,分母也要改变。

0.2°计入0.4°计入0.6°计入7.0°排除中心点
有效邻点数3 / 4

分子 = 0.2 + 0.4 + 0.6

KAM = 1.20° ÷ 3

0.400°
结果变小,是材料真的变形少了吗?

这次只是改变分析阈值或课堂邻点,不能说明材料状态变化。实际分析还要固定步长、核范围、相和晶粒限制,并检查原始质量。没有有效邻点时记为缺失,才能区分“没法计算”与“测得零错取向”。

四近邻教学核,数字用于手算,不是实验数据。真实二维 KAM 需按所用网格与算法构建邻域。

04

从连续取向场推到步长效应

在小范围、小旋转条件下,用旋转向量 θ(x) 表示局部取向偏离。邻点距离为 rj,有 θ(x+rj)−θ(x)≈κrj,其中 κ=∂θ/∂x 是旋转梯度。于是 ωij≈||κrj||,KAM≈邻域中 ||κrj|| 的平均。

若真实梯度不变,把所有邻点距离放大两倍,理想信号 KAM 也约放大两倍。对于一维线性旋转 θ(x)=kx,间距 h 的左右邻点各给 |k|h;若一阶核还含与梯度正交的上下邻点,平均还会受几何因子影响。这解释了不同邻域形状不能直接互换。

05

为什么同样一级核不能跨步长比较

步长 0.1 μm 的一级核与步长 1 μm 的一级核看的是不同物理尺度。后者可以跨过更多真实取向变化,也更可能跨过亚晶结构。保持相同像素阶数没有保持相同测量尺度。

比较两组图时,可按相同物理邻域规划分析,并讨论有效空间分辨率;或者在明确模型下考察角变化随距离的趋势。简单除以步长也不能解决所有问题,因为噪声、边界筛选、有限分辨率和多方向梯度仍在。

06

角噪声为什么使零梯度也有正 KAM

测量取向可理解为真实取向叠加小误差。两点的误差相减后即使均值为零,其长度或绝对值的平均仍为正。KAM 平均的是非负角度,正负误差不会像带符号数那样抵消。

若真实角变化接近噪声水平,KAM 主要反映测量能力。提升花样质量、可靠索引和适当降噪有帮助,但降噪也可能削弱真实梯度。Kamaya 的研究专门讨论测量准确度与局部梯度定义,适合在比较不同设置前阅读。

07

KAM、塑性应变与 GND 之间差哪一步

塑性应变是累计形变,KAM 是当前取向场在指定邻域内的不均匀性。位错经过并离开可造成塑性形变,却不一定留下同样的梯度;正负位错相互抵消也可使梯度小。因此相同塑性应变可能得到不同 KAM。

量纲上,在非常简化的晶格弯曲模型里可写 ρGND≈C·KAM/(bL),KAM 必须用弧度,b 为 Burgers 矢量长度,L 为物理距离,C 取决于几何假设。这个估算丢掉旋转轴和梯度分量信息,不等同于完整 Nye 张量反演;C 更不能作为所有材料共用的常数。

08

教程:正确比较一组 KAM 图

先固定相、坐标、晶粒规则和未索引处理;记录步长与核的物理尺度;给出 θcut、是否限制同晶粒、是否使用权重;保留原始质量图;用同一色标展示地图与分布;说明平均是按像素、面积还是晶粒。对晶界附近结果,检查有效邻点数量是否系统变化。

若退火后 KAM 降低,可以先报告在同一分析条件下粒内局部错取向减小。若要进一步说回复、再结晶或位错减少,需结合组织、边界、GOS、热处理与其他证据。不能把色条上 1° 读成 1% 应变。

09

为什么邻域大小有时比色图更重要

一个宽缓转动区和一条很窄的强梯度带,可能在某个核尺度下得到相近平均值。改变核尺度后两者趋势可能不同。KAM 是一个压缩后的局部统计量,它没有保存完整取向场。

读论文时把邻域、步长、阈值和清理方式写在同一行,再看作者的物理解释。若这些条件缺失,跨论文比较绝对 KAM 值通常缺少基础。

来源与延伸阅读

以下链接用于核对概念、标准状态或原始文献。本页为中文教学整理,不复制受版权保护的标准正文。

Olaf Engler; Stefan Zaefferer; Valerie Randle:《Introduction to Texture Analysis: Macrotexture, Microtexture, and Orientation Mapping》(2024,第 3 版)第 9.4.4 节,印刷页 315–319,KAM 和变形图。MTEX:Kernel Average Misorientation官方 KAM 定义、邻域与阈值。 官方定义、邻域阶数、阈值和晶粒 ID 的作用。示例参数属于示例材料。Masayuki Kamaya (2011). Assessment of local deformation using EBSD: quantification of accuracy of measurement and definition of local gradient局部梯度及测量准确度原始研究。Stuart I. Wright; Matthew M. Nowell; David P. Field (2011). A Review of Strain Analysis Using Electron Backscatter DiffractionEBSD 变形分析综述。J. Jiang; T.B. Britton; A.J. Wilkinson (2013). Measurement of geometrically necessary dislocation density with high resolution electron backscatter diffraction: Effects of detector binning and step size步长与噪声对位错估计的影响。