4. XRD定量与微结构分析

前几节我们完成了物相定性识别,知道样品里“有什么相”。但工程问题往往更进一步:各相占多少?晶粒多大?晶格有没有被拉压变形?晶粒有没有择优取向?这些才是决定材料性能的关键。本节围绕四条主线展开:Rietveld全谱拟合定量、晶粒尺寸与微观应变分离、织构判读、残余应力判读。它们共享同一套原始数据,却各自提取不同的微结构信息。

一、Rietveld精修:从全谱拟合到相含量

Rietveld方法的核心思想是:不再只用几个强峰,而是用整个衍射谱(含背景)做最小二乘拟合。模型谱由各相晶体结构(空间群、原子坐标、占位、温度因子)、仪器参数(峰形、零点、样品位移)和微结构参数(晶粒尺寸、应变、择优取向)共同计算得到,通过不断调整参数使计算谱逼近实测谱,最终用加权剩余因子Rwp、期望因子Rexp和拟合优度χ²= (Rwp/Rexp)²评价质量。工程上通常要求χ²接近1~2,Rwp在10%以内,且差值曲线无明显系统性起伏。

定量结果来自各相的标度因子Si,其质量分数近似为:

W_i = (S_i · M_i · V_i) / Σ(S_j · M_j · V_j)

S_i:第i相标度因子
M_i:第i相晶胞质量
V_i:第i相晶胞体积

注意:Rietveld给出的是“晶体学相含量”,若样品含非晶相,需另加内标(如刚玉α-Al₂O₃)才能得到绝对含量。实操中,精修顺序建议为:先修背景与零点→再修晶胞参数→再修峰形与微结构→最后修原子坐标与占位。切忌一次性放开所有参数,否则极易发散。

二、晶粒尺寸与微观应变:Williamson-Hall分离

衍射峰宽由仪器宽化、晶粒尺寸宽化和微观应变宽化共同贡献。先扣除仪器宽化(用标准LaB₆或Si标定),再对剩余宽化做分离。Scherrer公式给出尺寸贡献:

D = K·λ / (β_D·cosθ)      K≈0.9

而微观应变贡献:β_ε = 4·ε·tanθ

合并为 Williamson-Hall 形式:
β·cosθ = K·λ/D + 4·ε·sinθ

以4sinθ为横轴、βcosθ为纵轴作图,直线截距给出Kλ/D,斜率给出ε。若各点明显偏离直线,说明存在各向异性宽化,需改用修正W-H或Warren-Averbach方法。示例:某冷轧低碳钢,用Cu Kα(λ=0.15406 nm),测得(110)、(200)、(211)峰扣除仪器宽化后,拟合得截距0.0021 rad·nm,斜率0.0038,则D≈0.9×0.15406/0.0021≈66 nm,ε≈0.0038/4≈9.5×10⁻⁴。该结果说明晶粒已细化到纳米级且存在约千分之一量级的微应变,与冷轧加工硬化一致。

XRD微结构分析四条主线与数据流 原始衍射谱 2θ-强度 + 背景 ① Rietveld精修 全谱拟合 · 标度因子 Rwp / χ² 评价 → 各相质量分数 ② 尺寸/应变分离 扣仪器宽化 W-H 线性拟合 → D 与 ε ③ 织构判读 极图 / 反极图 取向分布函数ODF → 择优取向程度 ④ 残余应力判读 sin²ψ 法 峰位随ψ偏移 → σ 拉/压 工程输出:相含量 · 晶粒尺寸 · 微应变 · 织构强度 · 残余应力 用于工艺-性能关联、失效分析、质量判定 关键前提 仪器宽化标定 (LaB₆ / Si) 峰位准确 · 背景合理

三、织构判读:极图与取向分布

多晶材料中晶粒取向若随机分布,各(hkl)峰相对强度符合PDF卡片;若存在择优取向(织构),某些峰异常增强、某些减弱。判读织构最直接的工具是极图:固定某一(hkl),测量衍射强度随试样倾斜角χ和旋转角φ的变化,绘制等强度线图。若极图出现明显集中斑点或环,说明存在织构。更进一步可计算取向分布函数ODF,用欧拉角(φ₁, Φ, φ₂)空间描述取向密度,并给出织构体积分数。

工程判读要点:冷轧钢板常出现{111}或{100}织构,直接影响深冲性能;铝合金挤压棒材易出现<111>丝织构。若极图强度最大值与随机强度之比(即织构指数)超过2~3,就应视为有明显织构,定量相含量时需在Rietveld中引入择优取向修正(如March-Dollase模型),否则相含量会失真。

四、残余应力判读:sin²ψ法

残余应力通过晶面间距变化反映。当试样倾斜ψ角时,沿该方向测得的晶面间距dψ与无应力态d₀之差,与应力σ成正比:

ε_ψ = (d_ψ - d₀)/d₀ = [(1+ν)/E]·σ·sin²ψ - (ν/E)(σ₁+σ₂)

以 sin²ψ 为横轴、d_ψ 为纵轴作图:
斜率 m = [(1+ν)/E]·σ·d₀
则 σ = m·E / [(1+ν)·d₀]

若斜率为正,d随ψ增大而增大,对应拉应力;斜率为负则为压应力。工程上常用(211)峰测铁基材料、(311)峰测铝基材料。示例:某焊接接头热影响区,用Cr Kα测(211)峰,ψ取0°、15°、30°、45°,拟合得斜率m=+0.00042 nm,E=210 GPa,ν=0.30,d₀=0.1170 nm,则σ≈0.00042×210000/(1.30×0.1170)≈580 MPa,为拉应力,接近屈服强度,提示该区域存在开裂风险。若斜率非线性弯曲,说明存在应力梯度或织构干扰,需谨慎解释。

四条主线并非孤立:织构会改变峰强,影响Rietveld定量;残余应力会移动峰位,影响晶胞参数精修;晶粒尺寸与微应变又共同决定峰宽。因此实际分析应统一在一个精修模型中协同处理,而不是各做各的。
常见误区:直接用Scherrer公式算出的“尺寸”未扣仪器宽化,结果偏小;用单峰算应变,未做W-H分离,把尺寸贡献误算成应变;sin²ψ法未校正试样表面离焦和吸收,导致斜率失真。这些都会让微结构判读结论偏离真实。

综合来看,XRD定量与微结构分析的价值在于把“谱图”翻译成“材料状态”。Rietveld回答“各相多少”,W-H回答“晶粒与应变”,极图回答“取向”,sin²ψ回答“应力”。四者结合,才能对工艺—组织—性能关系给出有依据的判断。