织构分析:PF · IPF · ODF

ODF 取向分布函数:从体积分数定义到欧拉截面,完整推导一次

从微小取向范围占多少材料出发,推出 ODF 归一化、sinΦ 体积元、PF 与 ODF 的积分关系,以及 EBSD 核密度估计。

阅读设置

先记住这三点

  1. ODF 的自变量是完整取向,不是试样上的 x、y、z。
  2. 随机 ODF 为 1,但欧拉角空间中的等尺寸小盒子不具有相同旋转体积。
  3. 峰值、织构指数与组分体积分数回答不同问题。
01

先不谈复杂函数,只数材料体积

取一个小取向范围 Ω,例如与某个 Cube 理想取向相差不超过一个很小角度。把落在 Ω 中的晶粒体积加起来,除以总体积 V,得到该范围的体积分数。取向范围越来越小时,这个份额与范围大小之比,就成为取向密度。

写作 dV/V=f(g)dμ(g)。g 是一个完整取向,可由三个欧拉角表示;dμ 是旋转空间的归一化体积元;f(g) 就是 ODF。它非负,且对所有取向积分等于 1。这里 f 本身不是百分比,f 乘上取向空间体积才成为份额。

02

为什么体积元里会出现 sinΦ

完整取向可以分两步描述:先决定一根晶体轴朝向哪里,再决定绕该轴自转多少。轴的方向在球面上,球面面积元是 sinΦ dΦ dφ₁;绕轴自转提供 dφ₂。两者相乘,得到旋转体积元与 sinΦ dφ₁dΦdφ₂ 成正比。

在完整 SO(3) 表示中,φ₁、φ₂ 各取 0–2π,Φ 取 0–π。总积分为 (2π)×2×(2π)=8π²,所以 dμ=(1/8π²)sinΦ dφ₁dΦdφ₂。采用晶体或试样对称约化后,积分域和归一化必须相应处理;不能把完整空间的常数随意搬到一个缩小的角盒子里。

03

随机取向为什么不是三个角都均匀

随机取向要求对旋转空间均匀,即 f(g)=1。此时 Φ 的概率密度与 sinΦ 成正比,靠近 90° 的等角宽区间拥有更大旋转体积。若把 φ₁、Φ、φ₂ 分别均匀取值,就会在两极附近过度采样。

因此,判断一个欧拉散点云是否“均匀”,不能只看每个小方格里的点是否一样多。正确的比较要除以该小格的旋转体积。这个 sinΦ 因子同样进入 ODF 积分与纹理组分的定量计算。

跟着推导 · 教学几何

ODF:密度为什么不等于比例

完整 SO(3) 的 Bunge 参数范围中,固定 Δφ₁ = Δφ₂ = 60°,让一个宽 20° 的 Φ 小区间上下移动。Euler 图上同样大小的方块,并不代表同样多的随机取向。

0°45°90°135°180°随机取向的 Φ 边缘密度 ∝ sinΦΦ

dμ = sinΦ dφ₁ dΦ dφ₂ / (8π²)

μ(小区间) = (60/360) × [cos(35°) − cos(55°)]/2 × (60/360)

随机取向基线 = 0.341%;若局部 f ≈ 5,则该区材料比例 ≈ 1.705%

为什么 f = 5 不能直接读成 5%?

f 是单位取向空间体积上的相对密度。实际比例还必须乘该区域的取向空间体积,或者对变化的 f 做积分。这里仅计算一个小区域的局部积分,没有假设整个 ODF 都等于 5;完整 ODF 仍需满足积分为 1。

采用完整旋转空间以解释 Haar 测度;实际按晶体/试样对称约化后,须使用对应归一化,不能直接照搬这里的区间常数。示意曲线不是材料实测 ODF。

04

为什么 PF 是 ODF 的一类积分

固定晶体方向 h 和试样方向 y,所有满足 Rh=y 的取向仍可绕 y 自转,记自转角为 γ。PF 在 y 处把这些取向的贡献加在一起。在相容的单极、归一化约定下,可写 Pₕ(y)=(1/2π)∫₀²π f(g(γ))dγ;面族还要处理等价极的平均。

从三维取向分布沿一个自由度积分,剩下二维方向分布,这正是 PF 信息少于 ODF 的原因。反过来从 PF 求 ODF 就是一个反问题:许多不同的三维分布可能给出相近投影,需要多面族数据、对称性及约束共同限制。

05

EBSD 为什么可以走另一条路

EBSD 已经给出每个测点的完整取向 gᵢ,可以直接在旋转空间中进行统计。最直观的做法是分箱计数;更常用的是核密度估计 f̂(g)=Σᵢ wᵢKₕ(gᵢ⁻¹g)/Σᵢwᵢ,其中核 Kₕ 要在旋转空间归一,并按晶体对称处理等价取向。

wᵢ 决定统计对象:等面积网格逐点等权,反映截面面积权重;每晶粒平均取向等权,反映晶粒数权重;每粒均值配晶粒面积,则保留面积权重而舍弃晶内散布。三种做法都有用途,但必须说明,尤其在大晶粒与特定取向相关时。

06

核宽改变了什么,怎样避免只追求尖峰

核宽小,把每个取向分得更细,尖峰高、细节多,但对样本量与角噪声敏感;核宽大,把近邻取向平滑在一起,可能降低峰值、合并组分。它控制的是角分辨尺度,不是扫描步长。

做对比时先统一核类型、核宽定义与统计权重,再检查较窄和较宽设置下主要结论是否稳定。不能只挑使某个试样峰值最大的设置。软件中的 halfwidth、bandwidth、harmonic degree 也不是可以互换的同一个数字。

07

欧拉截面究竟怎样读

常见 ODF 图把 φ₂ 固定为 0°、5°、10°……,每张小图的横轴为 φ₁、纵轴为 Φ,颜色为 f。它们像一本三维分布的切片册,但切的是取向空间,不是材料厚度。相邻小图不是相邻空间位置。

读取一个峰,先写完整三元组 (φ₁,Φ,φ₂),再考虑对称等价和已知组分。只看到“φ₂=45° 截面很亮”不能命名为某种纤维,因为同一截面包含许多取向。应同时看峰随各截面的延伸,并用理想组分或纤维轨迹作参照。

08

把峰高与份额分开

一个非常窄的峰可以有很高的 f,但只占很小体积;较宽的峰虽然不高,却可能包含更多材料。组分体积分数应算 VΩ/V=∫Ω f(g)dμ,Ω 可以定义为到某理想取向的最小错配角不超过 α 的集合。

举一个纯教学例子:在旋转体积约为 0.01 的小区域内,若平均 f≈6,则贡献约为 0.06,即 6%。这里用的是区域平均值和区域体积,而不是把峰顶的 6 直接加上百分号。实际区域必须包含对称处理,并避免重复计数。

09

织构指数 J 从哪里来

定义 J=∫f(g)²dμ。在归一化体积为 1 的空间上,由 ∫(f−1)²dμ≥0 展开,得到 J−2∫f dμ+∫1dμ=J−1≥0,所以随机分布 J=1,非均匀分布通常 J>1。

J 衡量整体分布集中程度,却不知道峰朝向哪个方向。同样的 J 可以对应完全不同的各向异性。比较 J 还要固定平滑尺度与归一化;理想离散尖峰需要核平滑后才谈有限的数值。不能用“J 更大”替代性能分析。

10

一份完整 ODF 报告应该留下什么

说明相、晶体对称、试样对称、取样区域、取向数量与独立晶粒数量、权重、核或展开阶数、截面设置、色标、组分区域的角容差。若来自 XRD,还应提供测量的面族、覆盖角范围和反演回算残差;若来自 EBSD,说明未索引点和清理。

到这里,读者可以把三个对象重新串起来:取向矩阵描述一个晶粒;ODF 描述许多晶粒的完整取向分布;PF 和 IPF 是对指定方向的统计视图。它们之间可以互相计算某些量,但不会无损地彼此替代。

来源与延伸阅读

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