作者：孙千 本文转载自公众号：老千和他的朋友们。原文地址：https://mp.weixin.qq.com/s/E1rc2JumOOuDbuNm4RgX1Q
在金属、陶瓷、功能薄膜的微观表征日常工作里，电子背散射衍射（EBSD）是各大高校测试平台、课题组统一使用的数据分析工具。很多初次接触EBSD的学生完成SEM扫描、完成花样标定后，仅能看懂软件自带IPF取向衬度图，面对论文高频出现的极图、反极图、ODF 取向分布函数时常无从下手。
晶体织构本质是多晶内部全部晶粒空间朝向的统计分布规律，晶粒择优排布会直接改变材料弹性、导电、冲压成型、疲劳开裂等力学与物理性能的各向异性，也是连接加工工艺与材料性能的核心分析纽带。
多晶体由大量独立单晶颗粒拼接而成，每一颗晶粒都拥有专属的空间晶体朝向，EBSD 探测器采集到的扫描数据会存储截面上每一个点位对应的离散取向信息，我们可以对完整三维微观信息做逐层简化，剥离晶粒尺寸、形貌、空间位置这类次要信息，只保留决定材料性能的取向数据。
这种由繁至简的信息抽象过程可以结合商用奥氏体不锈钢重构结果直观理解，完整三维重构图像完整保留晶粒几何与取向全部细节，剔除空间形貌后只留存全部晶粒离散轴–角取向集合，最后再降维生成单一晶面投影的 {001} 极图。数据降维能够大幅压缩数十万量级的EBSD原始数据，降低读图难度，但过度简化会丢失关键取向特征。
图1 EBSD微观结构信息的逐级抽象。(a) 商用奥氏体不锈钢的三维微观结构重构结果，包含晶粒尺寸、形貌与空间取向信息；(b) 晶体材料织构的轴–角表示，完整表征所有存在的晶体取向，丢弃空间特征；(c) 具有铜型织构材料的 {001} 晶面赤平投影极图
一、旋转与取向的基础定义
在EBSD软件中设置样品轧制坐标系时，旋转约定是极易被新手忽略。学界存在主动旋转、被动旋转两套互不相同的物理定义，二者虽可数学换算，但空间逻辑完全区分，主流 EBSD 软件（Oxford AZtec、ESPRIT QUBE 、TSL OIM、MTEX）底层计算统一采用被动旋转规则，全混用两套逻辑会让所有取向数据系统性出错。
我们可以用生活化场景区分二者，主动旋转是坐标系固定不动，仅转动代表晶体的矢量，好比桌面保持原位，手动旋转放置其上的笔，物体自身发生空间偏转；被动旋转则固定晶体与矢量，只转动观测参考坐标系，如同笔静止在桌面，观测者转动座椅更换观察角度，晶体本身没有任何朝向变化。二者实现等效相对取向的条件为：被动旋转的顺序与主动旋转相反，且旋转角度大小相等、符号相反。
两种旋转等效变换的几何逻辑可以通过矢量分步变换示意图直观区分，分步主动旋转与分步坐标系被动旋转能达成完全一致的相对空间效果，从随动坐标系视角观察，被动旋转呈现的视觉效果等同于原旋转的逆变换。
图 2 主动约定与被动约定的旋转对比。(a) 主动约定：矢量初始沿 x 轴，先绕 z 轴旋转 π/2，再绕 x 轴旋转 π/2，最终与 z 轴重合；(b) 被动约定：坐标系先绕 x 轴旋转 π/2，再绕新生成的 z’ 轴旋转 π/2，与 (a) 实现等效的空间相对转动；(c) 从随动坐标系观察者的视角看，(b) 中被动旋转的表观效果为原旋转的逆操作
理清旋转规则后，我们就能读懂晶体取向的完整定义，预先设定一套标准参考晶体，令其晶轴与板材的RD（轧制方向）、ND（轧面法向）、TD（横向）完全重合，将这套标准晶体整体旋转至与实测晶粒朝向完全匹配，这一次完整旋转对应的正交变换矩阵，就是该晶粒对应的晶体取向。
另外，新手最容易困惑的取向对称冗余问题可以借助立方体骰子来理解，立方晶体具备完整立方点群对称性，同一套实际晶粒朝向能够对应24 种等价的旋转矩阵，骰子六面均可朝上，每一面朝上又能完成四次周向旋转，二十四种摆放方式在晶体学层面完全等价，对应数学表达OCj，Cj代表晶体点群对称操作。
如果材料经过轧制、拉伸塑性加工，样品会形成正交统计对称性，等价旋转矩阵数量会进一步增多，整体满足SiOCj变换形式。无需手动处理繁杂的对称等价矩阵，EBSD软件在生成极图、导出欧拉数据集时会自动完成对称归一化，筛选旋转角度最小、落在标准赤平三角内的唯一旋转作为对应取向代表，常规测试分析无需手动参与对称计算。
如果后续需要开展批量数据建模、织构数值模拟，才需要深入理解对称冗余背后的群论逻辑，立方、六方等高对称晶系取向重复现象明显，单斜、三斜低对称晶系几乎不存在等价取向，这部分仅作为拓展内容，仅依靠EBSD软件完成基础绘图的研究者无需深挖推导。日常软件操作里最常见的失误是切换织构模块、多次导出数据时随意改动旋转设置，建议每一组样品分析初期就锁定软件默认被动旋转坐标系，全程不改动底层旋转参数，避免因旋转约定不一致导致取向图谱整体偏移、失真。
二、极图与反极图：晶面取向的二维投影
极图是软件织构模块里使用频率最高、期刊论文配图最通用的二维织构表征手段，上手门槛最低，也是所有 EBSD 分析工作必须熟练掌握的核心内容。极图的底层逻辑依托晶面法向量与单位球面极点模型建立，任意晶面的空间朝向都能由其法向量唯一确定，法向量延伸至单位球上半球形成交点即为极点，对于立方晶系材料，{hkl} 晶面对应极点与<hkl>晶向几何等价。
为了把三维球面的极点分布转化为可直接打印、观察的二维平面图，行业统一采用赤平极射投影（又称吴尔夫投影）作为标准映射方式，成像逻辑可以类比透明地球仪，以单位球面南极作为光源，将北半球全部极点投射至北极平面，完成三维到二维的降维转换，整套球面投影几何模型能够清晰展示极点到平面点位的完整变换路径，简易投影变换公式写作
x,y,z为球面原始三维坐标，x’,y’是投影后平面坐标。
图 3 晶面取向的球面表示与赤平投影原理。晶面法向量与单位球面交于极点 p；通过位于南极的光源 Q，将上半球面上的极点 p 投影至 z=1 的投影平面上，得到投影点 p’，实现三维球面到二维平面的映射
赤平投影具备保角不保面积的固有特征，球面上任意两条直线的夹角经过投影后不会改变，方便我们直观判断晶粒取向之间的角度差异，但图谱边缘点位会出现视觉放大变形，类似世界地图高纬度陆地面积失真，仅凭肉眼观察点云疏密无法精准判定织构强弱，如果需要定量统计极点密度，可在EBSD软件极图设置面板一键切换等面积投影（又称施密特投影），代价是图谱内部角度关系会出现轻微畸变，定量分析类论文优先采用等面积投影的极图。
针对轧制板材样品，软件通常会默认以 ND 轧面法向作为投影中心轴、RD 轧制方向作为水平参考，图谱内部蓝色离散点位对应单颗晶粒极点，点云聚集程度越高，代表该晶面择优排布的晶粒占比越大，织构强度越强。
面心立方（FCC）金属的织构按加工状态分为形变织构与再结晶织构两类：冷轧变形过程中会形成铜型 {112}<111>、黄铜型 {110}<112>、S {123}<634 > 三类典型形变织构；经过再结晶退火后，会形成以立方织构 {001}<100 > 为代表的再结晶织构。
高层错能金属（如纯铝、多数铝合金）冷轧以铜型和S织构为主，低层错能金属（如奥氏体不锈钢、黄铜）冷轧后黄铜型织构占比显著升高。三类典型形变织构与立方再结晶织构的 {111} 极图拥有辨识度极高的分布特征，是我们出图后对照判断织构类型的标准参考样式。
图 4 三种典型织构的 {111} 极图。(a) 立方织构；(b) 铜型织构；(c) 黄铜织构。ND 与 RD 分别代表样品的法向与轧制方向
结合EBSD软件完整实操流程，导入完成标定、降噪处理的 EBSD 面扫数据集后，打开顶部Texture织构工具栏，调用 Pole Figure 绘图功能，自主勾选论文常用 {001}、{111} 特征晶面，按需切换赤平投影或等面积投影模式，一键生成矢量格式极图，可直接导出 tif、eps 格式用于论文排版。实操中建议统计点数不低于 1000 个有效晶粒，避免统计误差过大。
实操中需要留意极图本身无法规避的固有局限，即便同时生成多张不同晶面的极图，也无法建立图谱极点与单个晶粒完整三维取向的一一对应关系，会丢失大量三维取向细节，仅能完成定性或半定量观察。想要定量计算材料各向异性、统计各类织构组分占比，必须依托软件 ODF 取向分布函数模块完成分析。
和极图形成互补的表征工具是反极图，二者观测视角完全相反，新手极易混淆解读逻辑，我们可以借助教室场景通俗区分，极图站在样品整体视角，统计全部晶粒同一晶面朝向样品坐标系内的各个方位，如同统计教室所有学生面部朝向；反极图切换至晶体自身视角，统计样品特征轴向在晶体坐标系内对应的晶向，好比固定教室讲台，观察讲台对应每一位学生的正面、侧面还是其他方位。
反极图更适配线材、纤维、单向拉伸薄膜这类单轴向受力材料，板材样品一般生成法向反极图，用来判断轧面法向主要集中于 <100> 或是 < 111 > 晶向。
依托立方晶体二十四套等价对称空间，完整反极图可以简化为单一标准取向三角形（赤平三角）展示，大幅精简图谱排版，立方、铜型、黄铜织构的完整反极图与简化三角视图能够形成清晰对照，也是软件反极图模块默认输出样式。
图 5 三种织构的法向反极图。每列上方为完整反极图（含 24 个等价赤平三角，细线为三角边界），下方为标准赤平三角区域。(a)(b) 立方织构；(c)(d) 铜型织构；(e)(f) 黄铜织构
三、离散取向的参数化方法
单次完整 EBSD 面扫会产生数十万组独立离散取向数据，三维空间任意定点旋转仅需要三个独立参数就能完整描述，这也是学界衍生出多套取向参数化体系的底层逻辑，不同表达方案如同描述地理位置的经纬度、地标方位，不存在绝对最优形式，仅适配不同使用场景。其中Bunge 约定的 Z-X’-Z” 欧拉角是EBSD软件原生输出、金属织构论文通用标准，属于必须掌握的基础内容；轴–角、Rodrigues向量、四元数仅能将软件导出欧拉文本导入MTEX二次计算实现可视化，属于拓展内容，常规论文绘图无需深入学习。
EBSD软件全程采用 Z-X’-Z” 被动旋转序列定义 Bunge 欧拉角，整套取向由三次连续坐标系旋转构成，先绕样品原始 z 轴旋转φ1，再绕旋转生成的全新 x’ 轴旋转Φ，最后绕更新后的 z” 轴完成φ2旋转，三次旋转的几何变换逻辑可以通过欧拉角旋转示意图完整展示。
Bunge 欧拉角的原生完整取值范围为φ₁: 0°~360°，Φ: 0°~180°，φ₂: 0°~360°；针对同时具备立方晶体对称、板材正交加工对称的轧制样品，软件会自动压缩至 φ₁: 0°~90°，Φ: 0°~90°，φ₂: 0°~90° 的简化取向基本区，降低切片绘图的复杂程度。
欧拉角三维参数空间能够沿 φ₂轴均匀切片，生成成套的二维平面分布图，这也是金属织构领域沿用数十年的论文标准配图形式。我们可以在软件内导出纯文本欧拉数据集，导入软件后自定义切片间隔，批量输出成套 ODF 截面图谱。
图6 欧拉角的取向定义。取向由三次连续旋转确定，依次绕初始 z 轴转 φ₁、绕中间 x’ 轴转 Φ、绕最终 z” 轴转 φ₂
欧拉角体系存在无法规避的万向锁奇点问题，当中间旋转角度Φ=0时，第一次、第三次旋转轴重合，φ1与φ2失去独立参数意义，数值微分、拟合建模过程会出现计算不稳定，仅做定性看图不会产生明显影响，开展织构定量模拟时需要提前规避奇点对应的参数区间。
轴–角参数表达依托经典欧拉旋转定理建立，理论上任意空间旋转都等价于绕固定转轴 n 旋转固定角度 ω，整套表达整合为矢量形式a=nω，矢量方向对应旋转轴，矢量长度对应旋转角度，为避免重复表达，旋转角度仅取0至π区间，绕同一根轴正向、反向旋转 180° 会得到完全一致的晶体取向，对应轴–角空间为半径 π 的实心球体，球体直径两端点位代表同一套旋转操作。
这套表达几何直观性最强，课堂理论讲解更容易理解，但旋转角度接近 180° 时矢量参数会发生边界跳变，微小取向变化会造成数值跨球体突变，数值计算稳定性较差，同时两套连续旋转的换算公式繁琐复杂，EBSD软件无原生输出功能，仅用于课堂基础科普。
图7 轴 – 角表示与新欧拉映射的对应关系。(a) 任意系列旋转的最终效果，都等价于绕转轴 n 旋转角度 ω；转轴 n 可由球坐标 θ 和 φ 表示；(b) 新欧拉映射中，参数向量 s=nf(ω) 沿转轴方向，长度随旋转角单调变化
Rodrigues 向量在轴 – 角模型基础上完成数学变换，定义式为r=ntan(ω/2)，转轴方向保持不变，矢量长度替换为半旋转角的正切值，这套变换能够大幅简化旋转叠加的运算公式，依托晶体对称划分出边界全部为平面多面体的取向基本区，可沿坐标轴均匀切片实现二维可视化，立方晶系铜型织构在 Rodrigues 截角立方体空间内的分布与切片视图是这套表征的典型展示形式。
该方法核心短板是旋转角趋近 180° 时，tan(ω/2)数值趋向无穷大，无法用有限数值完整表达 180° 旋转，因此不作为EBSD软件原生输出格式，仅用于 MTEX 织构可视化二次开发。
图8 立方晶系铜型织构的 Rodrigues 空间基本区表示。(a) 含原点的立方对称基本区为截角立方体；(b) 沿 r₃轴等距切片的二维分布
四元数是四维归一化旋转表达工具，归一化标准形式写作q=[cos(ω/2),nsin(ω/2)]，标量与三维矢量满足归一约束q02+q12+q22+q32=1，几何上对应四维单位超球面，仅标量为正的半区即可覆盖全部晶体取向，+q与−q对应完全等同的物理旋转。
四元数体系不存在全局奇点，全程数值运算平滑稳定，是 3D 仿真、MTEX 脚本开发的核心工具，仅依靠简单换算r=q/q0就能和 Rodrigues 向量互相转换，常规依托EBSD软件完成论文绘图的新手无需掌握四维代数运算，只有自主开发织构模拟程序时才需要拓展学习。
图 8 立方晶系铜型织构的四元数表示。(a) 四元数向量部分 q 构成的三维空间取向分布；(b) 将三维空间按旋转角划分为同心球壳，分别对上半球做赤平投影后的二维呈现
四、取向分布函数：谐函数展开体系
EBSD软件输出的离散取向点是孤立分散的数据集合，无法直接用来定量计算材料弹性、塑性各向异性参数，取向分布函数ODF能够将离散取向点拟合为连续三维密度场，直观展示任意取向在材料内部出现的体积概率，是织构定量分析的核心工具，整套拟合逻辑和傅里叶曲线平滑、图像模糊降噪思路高度相似，依靠多组基础谐函数叠加还原完整连续的取向密度分布。
一维圆周分布依靠圆谐函数完成傅里叶展开，球面极点对应的二维极图分布依托球谐函数拟合，球谐函数由连带勒让德函数与圆谐函数组合构建，依靠阶数l、阶次m区分不同基函数，相当于球面上的傅里叶级数，专门用来拟合极图内部极点密度。
将软件导出的极点数据集导入MTEX后，可自定义最高展开阶数lmax控制图谱平滑程度：低阶展开会过度抹平织构细节，高阶展开容易生成虚假的密度负值与伪峰。实际分析中一般开启软件内置的正定性约束消除不合理的负密度值，同时常规分析推荐展开阶数取 lmax=16~22，兼顾细节保留与降噪效果。立方织构、铜型织构、黄铜型织构经过球谐展开生成的连续平滑极图，是定量分析类论文的标准配图样式。
图9 三种织构的连续 {111} 极图（球谐展开，l_max=15）。(a) 立方织构；(b) 铜型织构；(c) 黄铜织构。有限阶展开会引入少量伪密度，图中已施加正定性约束
目前EBSD软件内置 ODF 计算模块底层全部采用 Bunge 建立的广义球谐函数体系，依托欧拉角三维参数空间完成函数展开，广义球谐函数基础形式为Tlmn=eimφ2Plmn(cosΦ)einφ1，由两组指数函数与广义连带勒让德函数组合而成，展开系数为复数形式，铜型等轧制形变织构在欧拉角空间内会呈现标志性的管状连续分布，对应的三维等值面与沿φ2切片二维视图，也是软件 ODF 模块默认输出的经典形式。
ODF 图谱内 MRD（Multiple of Random Distribution，随机取向倍数）数值是定量分析核心指标，完全随机取向多晶材料 MRD 数值等于 1，数值越高代表对应取向晶粒体积占比越高、织构强度越强，依托 ODF 密度场可以直接计算各类特征织构体积分数，代入本构模型完成塑性变形、弹性模量各向异性数值模拟。
图10 铜型织构的连续欧拉角分布（广义球谐展开，l_max=12）。(a) 分布函数的单等值面，呈现清晰的管状形态；(b) 沿 φ₂轴切片的二维呈现
总之，结合EBSD 数据分析流程，我们可以梳理一套适配标准化的织构分析逻辑，完成样品 SEM 扫描、EBSD 花样标定与基础降噪处理后，首先通过软件原生 IPF 取向图观察晶粒尺寸、晶界、形貌等微观基础特征；如果仅需要定性识别晶粒择优取向类型，直接打开Texture织构工具栏调用极图、反极图绘图模块，对照 FCC 经典轧制织构图谱完成初步判断；需要开展定量性能分析、统计织构组分体积占比时，再启用 ODF 取向分布函数模块完成密度统计。
不同织构表征手段对应各自适配的使用场景，期刊论文定性结果配图优先选用可视化直观的极图与反极图；涉及力学、弹性性能定量计算、织构组分统计的内容，统一采用 Bunge 欧拉角空间下的 ODF函数作为核心分析依据。
对于初学者而言，无需吃透全部取向数学表达，优先熟练掌握EBSD软件内极图、反极图、欧拉ODF三大核心织构模块的操作与基础读图方法，就足以支撑绝大多数实验测试，待软件操作、图谱解读形成完整认知后，再根据自身科研方向按需学习进阶拓展内容，循序渐进降低晶体织构理论的整体学习门槛。
参考资料
Representations of Texture
Introduction to texture analysis
