2. 晶粒组织表征与工程需求

相场模拟晶粒长大，最终目的不是得到一张好看的彩色图，而是回答工程问题：这个材料的晶粒组织能否满足使用要求。因此，在建立模型之前，必须先明确"表征什么"和"要求什么"。本节围绕晶粒尺寸分布、取向差、拓扑特征三类核心表征量，以及实验获取手段和工程性能要求展开，为后续相场参数标定提供依据。

一、晶粒尺寸分布：平均值之外的工程信息

晶粒尺寸最常用的描述量是平均晶粒直径 d̄，但仅凭平均值无法区分"均匀细小"与"粗细混杂"两种组织。工程上更关注尺寸分布，通常用等效圆直径 d = 2√(A/π) 表示单个晶粒尺寸，其中 A 为截面面积。对三维组织，截面法测得的尺寸需乘以修正因子（如 1.5 左右）才能近似真实三维直径。

分布形态常用对数正态分布拟合，其概率密度为：

f(d) = 1/(d·σ√(2π)) · exp[-(ln d - μ)²/(2σ²)]

μ：ln d 的均值    σ：ln d 的标准差
几何平均 d_g = exp(μ)   几何标准差 σ_g = exp(σ)

例如某低碳钢实测 d_g = 18 μm、σ_g = 1.35，说明约 68% 的晶粒落在 13.3～24.3 μm 区间。若 σ_g 超过 1.6，则表明存在明显的异常长大晶粒，这类组织往往对应局部力学性能薄弱区。相场模拟中，初始晶粒结构应尽量复现实验测得的 d_g 与 σ_g，而不是只匹配平均值。

二、取向差与晶界特征

晶粒之间的取向差 θ 决定晶界能与迁移率。对立方晶系，常用最小取向差表示：

θ = arccos[(tr(R) - 1)/2]

R：两晶粒取向矩阵的相对旋转矩阵
θ 取使取向差最小的对称等价解

按 θ 可将晶界分为小角度晶界（θ < 15°）和大角度晶界（θ ≥ 15°）。小角度晶界能较低，迁移率也低；大角度晶界中，某些特殊取向差（如 Σ3 孪晶界，θ ≈ 60° 绕 <111>）具有极低界面能，常被称为"特殊晶界"。工程上通过晶界工程提高特殊晶界比例，可改善抗晶间腐蚀和抗蠕变性能。

相场模型中，取向差通常通过多序参量或取向场间接体现。若模型只追踪晶粒形状而不区分取向，则无法正确再现特殊晶界钉扎效应，这是建模时必须提前判断的简化边界。

三、拓扑特征：配位数与晶粒形状

拓扑特征描述晶粒之间的连接关系，最核心的是配位数（每个晶粒的相邻晶粒数）。二维理想正六边形组织的配位数为 6；实际组织中，配位数分布偏离 6 的晶粒往往处于长大或消失的竞争状态。根据 von Neumann 定律，二维单个晶粒面积变化率与其边数 n 相关：

dA/dt = -M·γ·(π/3)·(n - 6)

n > 6：晶粒长大    n < 6：晶粒收缩消失
n = 6：面积近似不变

该关系说明：拓扑边数直接决定晶粒的演化命运。相场模拟中可通过统计每个晶粒的邻居数验证模型是否合理——若模拟得到的配位数分布过于集中，往往意味着初始结构过于规则，与真实组织不符。

四、实验表征手段与数据获取

金相法是获取晶粒尺寸分布最直接的手段：抛光腐蚀后在光学显微镜下成像，用截线法或面积法统计。截线法在 n 条随机测试线上数晶界截点数 P，平均截距 L̄ = L_total / P，再换算为平均晶粒直径。该方法效率高，但难以区分取向差。

EBSD（电子背散射衍射）可同时给出晶粒形貌和取向信息，是获取取向差分布和特殊晶界比例的标准手段。典型输出包括：取向成像图（IPF）、取向差角分布直方图、晶界特征分布（GBCD）。对 Σ3 晶界，EBSD 可依据 Brandon 判据 Δθ ≤ 15°/√Σ 进行判定。三维组织则需 3D-EBSD 或 FIB 断层扫描，成本高但能避免截面统计偏差。

五、工程材料对晶粒组织的性能要求

不同工程场景对晶粒组织的要求差异显著，下表给出典型示例（参数为教学示例，实际以具体标准为准）：

这些要求最终转化为相场模型的输入约束：初始晶粒尺寸分布、取向差分布、晶界能各向异性参数。若模拟目标为"细化晶粒"，则需在模型中引入第二相粒子钉扎或溶质拖曳；若目标为"提高 Σ3 比例"，则需在晶界能函数中体现特殊晶界的低能特征。表征数据越完整，模型标定越可靠。