4. 位错蠕变机制

当温度升高到约0.4Tm以上时,原子扩散能力显著增强,位错可以通过攀移绕过障碍物,蠕变行为从扩散主导转向位错运动主导。位错蠕变是高温合金、耐热钢和陶瓷材料在中高应力区间最主要的变形机制,理解其速率方程、控制环节和应力依赖关系,是建立寿命预测模型的基础。

一、Orowan方程:蠕变速率的基本框架

位错蠕变的宏观应变速率可以用Orowan方程描述:

ε̇ = ρ_m · b · v̄

其中ε̇为稳态蠕变速率,ρm为可动位错密度,b为柏氏矢量长度,v̄为位错平均运动速度。该方程表明,蠕变速率由两个因素共同决定:有多少位错在运动,以及它们运动得多快。位错在滑移面上运动时,会遇到三类主要障碍:林位错(其他滑移系上的位错)、析出相颗粒和晶界。位错不能直接切过这些障碍,必须通过攀移或交滑移绕过,而攀移需要空位扩散,因此高温下位错运动速度受扩散控制。

以Ni基单晶高温合金CMSX-4为例,在1050°C、150 MPa条件下,典型稳态蠕变速率约为1×10-8 s-1。若测得可动位错密度约为1×1012 m-2,柏氏矢量b≈0.25 nm,则可由Orowan方程反推位错平均速度约为4×10-8 m/s,即每秒仅移动约0.04微米。这一数量级说明,高温蠕变中位错运动极其缓慢,速率瓶颈不在位错滑移本身,而在绕过障碍物的等待时间。

二、回复控制:蠕变的速率决定环节

位错在滑移过程中不断增殖和堆积,若没有回复过程,加工硬化会使蠕变迅速停止。回复通过位错攀移和异号位错湮灭来降低位错密度、软化材料,使变形得以持续。因此,稳态蠕变速率实际上由回复速率控制,而回复速率又取决于空位扩散系数和层错能。

回复控制蠕变的一个核心推论是:蠕变激活能应接近自扩散激活能。对于纯铝,蠕变激活能约为142 kJ/mol,与铝的自扩散激活能126~142 kJ/mol吻合;对于α-Fe,蠕变激活能约为250 kJ/mol,与其自扩散激活能251 kJ/mol一致。这一规律成为判断蠕变是否由位错攀移控制的重要实验依据。若实测激活能明显低于自扩散激活能,则可能由管道扩散或晶界扩散主导;若明显偏高,则可能存在析出相粗化等附加过程。

层错能越低,位错交滑移越困难,回复速率越慢,蠕变速率也越低。因此,Co基合金中通过降低层错能可以提高蠕变抗力,这也是合金设计的重要思路。

三、门槛应力:蠕变是否发生的判据

位错蠕变并非在任何应力下都会发生。当外加应力低于某一临界值时,位错无法克服障碍物的阻力,蠕变速率趋近于零,该临界值称为门槛应力σth。门槛应力主要来源于析出相颗粒的Orowan阻力、固溶原子的拖曳效应和晶界滑移阻力。

对于Orowan绕过机制,门槛应力可估算为:

σ_th ≈ G·b / (2π·λ)

其中G为剪切模量,λ为颗粒间距。以镍基合金中γ′相为例,若G≈80 GPa,b≈0.25 nm,λ≈100 nm,则σth≈32 MPa。这意味着在低于约32 MPa的应力下,位错无法绕过γ′颗粒,蠕变速率极低。工程上常利用门槛应力来设定合金的长期服役应力上限,例如定向凝固合金在1000°C下的设计应力通常控制在80~120 MPa,远高于门槛应力但低于导致快速蠕变的临界值。

四、幂律蠕变与指数律蠕变

当外加应力高于门槛应力时,稳态蠕变速率与应力的关系通常用幂律方程描述:

ε̇ = A · σ^n · exp(-Q/RT)

其中A为材料常数,n为应力指数,Q为蠕变激活能,R为气体常数,T为绝对温度。应力指数n反映了蠕变机制:n≈3~5通常对应位错攀移控制;n≈5~7可能对应位错交滑移或析出相强化;n≈1~2则可能对应扩散蠕变或晶界滑移。

当应力继续升高到接近屈服强度时,幂律关系不再成立,蠕变速率随应力呈指数增长,进入指数律蠕变区:

ε̇ = B · exp(β·σ) · exp(-Q/RT)

其中β为指数律常数。指数律蠕变的物理机制是:高应力下位错无需等待攀移即可直接越过障碍,或发生热激活交滑移,变形速率对应力极为敏感。工程上,指数律区往往对应蠕变寿命急剧缩短的危险工况,应尽量避免。

位错蠕变机制与应力-速率关系 Orowan方程 ε̇ = ρ_m · b · v̄ ρ_m:可动位错密度 b:柏氏矢量 v̄:位错平均速度 速率瓶颈:绕过障碍的等待 回复控制 位错攀移 + 异号位错湮灭 Qcreep ≈ Qself-diffusion Al:142 kJ/mol ≈ 126~142 α-Fe:250 kJ/mol ≈ 251 层错能↓ → 回复慢 → 蠕变抗力↑ 门槛应力 σ_th σ_th ≈ G·b / (2π·λ) G ≈ 80 GPa,b ≈ 0.25 nm λ ≈ 100 nm → σ_th ≈ 32 MPa σ < σ_th:蠕变速率趋近于零 工程应力上限设定依据 稳态蠕变速率与应力的双对数关系 σ (对数) ε̇ (对数) σ_th 幂律区 n ≈ 3~5 位错攀移控制 指数律区 热激活交滑移 扩散区 指数律区:寿命急剧缩短 ε̇ = B·exp(βσ)·exp(-Q/RT) 工程上应避免进入此区 幂律区:设计工作区 ε̇ = A·σ^n·exp(-Q/RT) n值判断蠕变机制

五、机制判据与工程应用

在实际工程分析中,判断蠕变机制需要综合多个信息。首先看应力指数n:若n在3~5之间,且激活能接近自扩散激活能,可判定为位错攀移控制的幂律蠕变。其次看门槛应力:若服役应力低于门槛应力,则蠕变寿命可能极长,设计可偏保守;若接近或超过门槛应力,则需按幂律方程进行寿命外推。最后看微观组织:若析出相粗化速率快,λ增大,σth下降,蠕变速率会随时间上升,出现加速蠕变阶段。

以某火电机组用P91钢为例,在600°C、80 MPa条件下,实测应力指数n≈5,激活能Q≈550 kJ/mol,明显高于α-Fe自扩散激活能,说明存在析出相强化和位错亚结构演化的附加贡献。此时若简单用纯攀移模型外推寿命,会高估蠕变速率、低估寿命。工程上更可靠的做法是采用Larson-Miller参数法或Omega法,结合多温度多应力短时试验数据拟合,再外推到服役条件。

门槛应力并非固定值,它随析出相粗化、位错亚结构变化而降低。长期服役后,原本安全的应力可能超过新的门槛应力,导致蠕变速率突然上升。因此,寿命评估必须考虑组织退化对门槛应力的影响。

总结来看,位错蠕变的核心逻辑链是:Orowan方程给出速率框架,回复控制决定速率瓶颈,门槛应力划定蠕变发生的下限,幂律与指数律分别描述中应力和高应力区的速率-应力关系。掌握这四个环节,就能在给定温度、应力和组织状态下,判断蠕变机制、估算稳态速率并评估寿命风险。