钢中奥氏体晶粒长大规律

钢中奥氏体晶粒长大规律
刘文月;任毅;王爽;张帅;高红
【摘要】评述了奥氏体正常长大过程中影响晶粒尺寸的因素、晶粒尺寸计算常用
公式与晶粒尺寸确定方法等,归纳了钢中奥氏体长大的定性与定量规律.目前计算奥
氏体晶粒尺寸的公式多为经验公式,尚不具备预测的功能.
【期刊名称】《上海金属》
【年(卷),期】2019(041)004
【总页数】6页(P88-93)
【关键词】奥氏体;晶粒长大;Zener-Hollomon参数;Zener方程
【作者】刘文月;任毅;王爽;张帅;高红
【作者单位】海洋装备用金属材料及其应用国家重点实验室,辽宁鞍山114009;鞍
钢集团钢铁研究院,辽宁鞍山114009;海洋装备用金属材料及其应用国家重点实验室,辽宁鞍山114009;鞍钢集团钢铁研究院,辽宁鞍山114009;海洋装备用金属材料
及其应用国家重点实验室,辽宁鞍山114009;鞍钢集团钢铁研究院,辽宁鞍山114009;海洋装备用金属材料及其应用国家重点实验室,辽宁鞍山114009;鞍钢集
团钢铁研究院,辽宁鞍山114009;海洋装备用金属材料及其应用国家重点实验室,辽
宁鞍山114009;鞍钢集团钢铁研究院,辽宁鞍山114009
【正文语种】中文
Fe- C平衡相图中，温度超过727 ℃，在整个钢的成分范围内都可能存在奥氏体。

向Fe- C二元合金中，添加不同种类的合金元素，可以改变奥氏体存在的温度区间与成分范围[1]。

奥氏体晶粒长大规律符合Avrami法则，即晶粒尺寸是温度的幂
指数函数，是时间的抛物线函数[2- 3]。

奥氏体晶粒的长大不是没有限制的，具有一个极限尺寸，Zener认为该尺寸由钉
扎粒子的尺寸和数量决定[4- 6]。

奥氏体晶粒长大过程中，会出现正常长大与异常长大两种现象。

正常长大过程中，最大晶粒尺寸约为平均尺寸的2.5倍[7]；异常
长大过程中，仅有少数晶粒尺寸急剧增加，最大尺寸一般为平均尺寸的6倍以上，有时可达十几倍[8]。

奥氏体晶粒的异常长大，也被称为二次再结晶[9]。

细化奥氏体晶粒，可以显著改善钢材的力学性能。

在奥氏体单相区，对奥氏体施加一定量的塑性变形使之发生再结晶从而细化晶粒，是现代钢铁产品获得优异性能的常规工艺。

无论是原奥氏体晶粒还是完成再结晶后的奥氏体晶粒，都具有相同的长大规律。

了解奥氏体晶粒的长大规律，对控制钢材的最终性能具有重要意义。

现有预测晶粒尺寸的数学公式皆为经验公式，存在需要通过试验或模拟计算才能确定的未知常数，因此，需要进一步寻找更理想的模型来定量计算奥氏体晶粒尺寸。

本文仅评述奥氏体正常长大过程中影响晶粒尺寸的因素、晶粒尺寸计算常用公式与晶粒尺寸确定方法等，尝试获得钢中奥氏体长大的定性与定量规律。

1 晶粒长大基本规律
奥氏体晶粒长大受Fe、C原子的扩散速率控制，其长大规律有如下定性认识[9]：1)晶粒长大是通过晶界移动而不是相邻晶粒合并的方式；2)晶界移动是不连续的，移动速率与方向都会发生改变；3)同一个晶粒的晶界可以向不同的方向移动，即部分晶界增加晶粒尺寸的同时，也允许另一部分晶界减小该晶粒的尺寸；4)晶界移动方向通常指向其曲率中心；5)二维晶界夹角的稳定数值为120°，非此角度晶界通
过晶界弯曲与移动达到此稳定值；6)溶质原子与第二相粒子具有改变晶界移动速度的作用。

晶粒长大是分阶段的，不同阶段有不同的长大机制。

Simpson等[10]在纯Pb、
Cd及Pb基、Cd基低合金中，观察到了晶粒的四阶段生长模式。

在奥氏体晶粒的等温生长过程中[11]，晶粒生长一般都以具有最低激活能的方式首先生长；温度越低，激活能越小的生长机制持续时间越长；然后随着温度升高依次过渡到激活能更高的生长机制；高温时激活能低的生长机制很快完成而失去作用。

高振英等[12]研究了中碳钢奥氏体晶粒的4种生长机制：1)Fe原子位错扩散控制机制；2)Fe原子位错、界面扩散控制的混合机制；3)界面扩散控制机制；4)Fe原子体扩散控制机制。

2 影响晶粒长大的因素
影响奥氏体晶粒长大的因素众多，主要有材料类型、温度、时间、长大机制、溶质原子数量与尺寸、钉扎粒子数量与尺寸及温度梯度等。

晶粒的等温长大过程可以不考虑温度梯度的影响。

而在非等温长大过程中，温度梯度对晶粒生长的影响较大。

定量描述奥氏体晶粒长大规律，特别是区分各个影响因素对晶粒长大的作用，比较困难。

1948年，Beck等首次给出了计算晶粒尺寸的经验公式[7,13]：
(1)
式中：D为平均晶粒尺寸，μm；k为与材料和温度有关的系数；t为保温时间，s；k0为材料常数；Q为激活能，J/mol；R为气体常数，8.314 J/(mol·K)；T为绝对温度，K。

从Beck方程中可以看出，晶粒尺寸除了与材料有关，还与温度、时间和激活能有关。

温度、时间对晶粒尺寸的影响满足Avrmi方程。

激活能对晶粒尺寸的影响是
一个综合性概念[14]，它包含了溶质原子、钉扎粒子、长大机制、温度梯度等诸多因素的影响。

不同溶质原子对奥氏体晶粒长大的影响不同[1,8,15- 17]，如表1所示。

Moon等
[2]、Lee等[16]的研究结果证实了Mo、Cr、Ni对奥氏体晶粒长大的抑制作用，
提供了C、Mn抑制奥氏体晶粒长大的例证。

表1 溶质原子对奥氏体晶粒长大的作用[8,15]Table 1 Effect of solute atoms on austenite grain growth [8,15]作用元素促进长大C、Mn、P微弱抑制Si、Ni、Cu[8]中等抑制Mo、Cr、W强烈抑制Nb、V、Ti、Al、Zr[8]、N、O
除溶质拖曳外，粒子钉扎是影响晶粒长大的最重要因素，可以用Zener方程[5,18]来描述：
(2)
(3)
式中：Rc为临界晶粒尺寸；r为粒子尺寸；f为粒子体积分数；K为Zener模型试验解析解系数；δ为指数常数；FZ为粒子钉扎力；KZ为Zener系数，γ为界面能。

由Zener方程可见，材料的临界晶粒尺寸Rc由材料中第二相粒子(包括夹杂物)的尺寸r和体积分数f决定[5]。

因此，第二相粒子尺寸越小，数量越多，粒子钉扎力FZ越大，则材料的平均晶粒尺寸越小[19]。

Zener方程描述的粒子钉轧力，会随着温度、时间的变化因粒子尺寸与数量的变
化而变化[1,19- 21]。

随着温度的升高，粒子会发生析出、长大和溶解等一系列变化。

钉扎粒子数量与尺寸的变化，导致在不同温度区间和时间段内的粒子钉扎力大小不同，这是Zener方程中系数K与指数δ取值多变(表2)的一个原因[1,5]。

粒子的溶解温度可以用固溶度积公式计算[26]：
(4)
式中：KS为平衡常数；X为微合金元素的质量分数，%；Y为碳或氮的质量分
数，%；C1、C2为常数。

表2 不同指数δ取值时Zener模型试验解析解系数K的取值[1,5]Table 2 Analytical solutions of factor K in Zener model test with different exponential δ values[1,5]δ值K值10.05～0.67; 0.73～0.91; 4/3 (1.33); 0.88～2.83∗0.932/9(0.22)0.921.540.876.10.550.60.51.03; 1.18; 1.41; 1.45; 1.7; 30.5(1.73); 1.8; 20.331.80.31±0.024.5±0.8
*当钉扎粒子为TiN共格粒子时，K取值为2.10
由于试验和推导的基础条件不同，对于同一种粒子所得到的固溶度积公式也不同[22] (见表3)。

在Nb的质量分数为0.021%的合金中[23]，因Nb/TiN之间的交互作用，不同温度下Nb固溶量的计算值是试验值的十几倍，甚至数十倍以上(见表4)。

除了公式本身的缺陷外，取样或制样不当也会给结果带来误差。

1956年，Beck 首先报道了晶粒长大过程中的“板厚效应”[7]，即当晶粒长大到一定尺寸时，其长大速度会因试样尺寸过小而减小。

若出现这种情况，会造成激活能的计算值严重偏大。

为了避免这种情况，建议试样的最小截面尺寸至少为试验后平均晶粒尺寸的十倍。

另外，为了消除观测误差，可通过冷加工的方式将样品表面去掉一部分，以消除试验过程对试样表面的影响，如脱碳、氧化等。

如果样品尺寸允许，建议至少去掉2 mm。

表3 Nb溶解温度计算公式[22]Table 3 Formula for calculating Nb dissolution temperature[22]公式名称[Y]C1C2T/℃Meyer[C]3.047 2901
112Mori[C]0.873.187 7001 177Irvine[C+12/14N]2.166 7701
267Siciliano[C+12/14N]2.266 440-838[Mn]0.246+1 730[Si]0.5941
245Dong[C+12/14N]3.14+0.35[Si]-0.91[Mn]8 049-1 371[Mn]+923[Si]1 304Mori[N]0.65[C]0.244.0910 4001 297Meyer[N+C]1.545 8601 326
表4 Nb的质量分数为0.021%的合金中Nb固溶量的计算值与试验值[23]Table 4 Calculated and experimental contents of solid solution Nb in alloys with
0.021% Nb by mass[23]奥氏体化温度/℃计算值/%试验值/%1 2500.021 00.003 21 2000.013 30.000 99000.000 00.000 5
在奥氏体晶粒的非等温长大过程中，除了以上影响因素，还有其特有的影响因素。

钢中奥氏体晶粒的非等温长大主要发生在铸造、热轧与焊接等过程中。

与晶粒等温长大相比，非等温长大由于温度场的存在而具有不同的特征(见表5)[4]。

表5 非等温过程晶粒长大特征[4]Table 5 Characteristics of grain growth during non- isothermal processes[4]铸造热轧焊接(焊接热影响区)温度变化缓
慢变化缓慢或快速变化高温度梯度与强热循环典型特征完全柱状或等轴状或两者混合晶粒结构;浇注温度低与温度梯度小,促进柱状晶向等轴晶转变加热速度慢,晶粒尺寸大晶粒尺寸空间差异明显;晶粒尺寸分布与拓扑类型分布不受空间影响;热钉扎源
于高温度梯度
3 常用的描述奥氏体晶粒长大的经验公式
为了描述方便，本文将与奥氏体晶粒长大相关的数学公式中所采用的符号做了统一。

Sellars等[24]引入了初始晶粒尺寸的影响，将Beck方程(式(1))做了等价改写(式(4))，改写后的形式已被广泛接受[1,5,7]。

近年来，研究人员[25]联立Beck模型
与Sellars模型，获得了第3种形式的Beck- Sellars模型(式(5))。

(5)
式中：Dt为t时刻平均晶粒尺寸；D0为起始时刻平均晶粒尺寸；n、n′为与材料、温度有关的指数；k为与材料、温度有关的系数；t为保温时间。

Beck、Sellars与Beck- Sellars 3种模型都存在与材料和温度有关的未知常量，为
求得这些常量，需获得材料在不同温度和保温时间条件下的平均晶粒尺寸，然后对试验结果回归分析，求得这些未知常量的试验解，最终得出计算平均晶粒尺寸的经验公式。

任何满足Avrami方程的晶粒长大模式，都可用这些模型来描述[1,13]。

在实际应用中，Beck、Sellars与Beck- Sellars 3种模型公式都需要做一定的变换才能使用。

为了区分材料与试验温度的影响，将参数k、M改写为Arrhenius方
程的形式(式(6)、式(7) )。

对式(1)求导数，得到晶粒长大速度式(8)；式(9)是晶粒
长大速度的另一种表现形式，联立这两个方程，可得m与n之间的关系式(式10)。

考虑Zener钉扎力的影响，可得式(11)。

式(1)、式(5)、式(8) 和式(11)都是较为典型的晶粒尺寸表达形式，其典型取值如表6所示。

(6)
(7)
(8)
(9)
(10)
(11)
式中：k、M为与材料、温度有关的系数；k0、M0为与材料有关的系数；n、m
为与材料、温度有关的指数；Q为激活能，J/mol；R为气体常数，8.314
J/(mol·K)；T为绝对温度，K；D为平均晶粒尺寸；ν为晶界迁移速率；ΔF为晶界
驱动力；γ为晶界表面能；r为粒子尺寸；f为粒子体积分数；KZ为Zener系数。

式(6)～式(11)没有考虑应变、应变速率对奥氏体晶粒尺寸的影响，因此，再结晶
晶粒尺寸Ddrx的计算，通常使用Zener- Hollomon参数[18,26]来构造(式(12)、式(13))，它们分别是Sellars幂指数模型[18,24,27]与Sekine半对数模型，前者
获得了广泛的试验支持。

Ddrx=cZ-nd
(12)
Ddrx=lgaZ-b
(13)
式中：Ddrx为再结晶平均晶粒尺寸；nd、a、b、c为常数；Z 为Zener- Hollomon参数，s-1。

获得Ddrx首先要确定Z参数，Z参数的确定需测得材料在不同温度和应变速率下的应力- 应变曲线，然后对曲线进行拟合得到Z参数中未知的试验常数，再对Z参数与晶粒尺寸Ddrx进行拟合，确定Sellars幂指数模型与Sekine半对数模型中
的未知试验常数，即可获得晶粒尺寸Ddrx与Z参数之间的函数关系。

Sellars幂
指数模型是用于计算动态再结晶晶粒尺寸的使用最广泛的经验公式，表6给出了
模型中系数c与指数-nd的一些取值实例。

Z参数有3种形式[24]，分别适用于不同的应力范围，式(14)适用于整个应力范围，式(15) 适用于低应力范围(ασ<0.8)，式(16) 适用于高应力范围(ασ>1.2)。

式(14)～式(16)不能直接用于计算再结晶晶粒尺寸，需要辅助一定的公式形式变化与试验。

首先取应力峰值σp代入式(14)，获得应变速率的表达式(式(17))。

对式(17)求偏导，得到式(18)，再结合不同温度和应变速率下的应力峰值，通过线性拟合获得q、α、Qd的试验解，将这些试验解代入式(17)可得A的试验解。

至此，
式(14)和式(17)获得了可计算的表现形式。

同理，对式(15)和式(16)两边取对数，
通过线性拟合试验数据，可以获得A1、A2、q1、β的试验解，且β与α、q1存在如式(19)所示的限定关系。

表6 描述奥氏体晶粒长大的经验模型试验解Table 6 Experimental results on empirical models describing austenite grain growth等式形态nk0Q/(kJ·mol-1)备注D1nt-D1n0=k0tn′exp-
QRT()0.371.23×105161.8n'=0.404[28]1.235.548×1014385.2n'=0.897[29]D1 nt-D1n0=k0texp-
QRT()0.100.100.500.140.500.100.130.140.143.87×10325.02×10321.44×1012 1.45×10270.72×10122.60×10284.20×10271.97×10221.30×1023400 914 267 400 267.1437 400 275 308 T>1 273 K [20]T≤1 273 K[20]文献[20]文献[30]Dt=k0tnexp-QRT()0.217.67×10589+3.58C+1.21Ni+1.44Cr+4.03Mo文献[16]0.501.43×105128.3文献[11]dDdt=nk0exp-QRT()1D-FZ[]1n-
10.241.72×1021352.1+21.83C+20.0Mn+7.19Cr+7.38Ni文献[1,17]
表7 Sellars幂指数模型中D- Z关系式试验解Table 7 Experimental relationship between D and Z in the Sellars power exponent model钢种cndAISI 304[27]5.2×1030.17S45CVMn[18]1.204×1020.056C-
Mn[31]1.6×1040.23
获得Z参数的可计算形式后，可将材料的应变与峰值应力描述为式(20)、(21)，通过热变形试验获得试验拟合解后，即可建立该材料高温变形的本构方程，用于描述材料在再结晶阶段的应力- 应变关系。

结合式(12)、式(13)、式(20)与式(21)即可获得奥氏体再结晶晶粒尺寸。

(14)
Z=A1σq1
(15)
Z=A2exp(βσ)
(16)
(17)
(18)
β=α·q1
(19)
ε=A′Zq′
(20)
σp=A″Zq″
(21)
式中：为应变速率；Qd为激活能，J/mol；R为气体常数，8.314 J/(mol·K)；T 为绝对温度，K；A、A1、A2、A′、A″、α、β、q、q1、q′、q"为材料常数；σ为流变应力；σp为应力峰值；ε为应变。

4 结语
奥氏体晶粒长大遵循扩散控制机制，其晶粒尺寸是温度的幂指数函数，是时间的抛物线函数。

钢中奥氏体晶粒长大是分阶段的，不同阶段有不同的控制机制，表现为Fe、C原子扩散方式的改变。

晶粒尺寸除了与溶质原子有关，还与温度、时间和激活能有关，另外第二相粒子对晶粒的长大速度也有一定影响。

为了定量描述这些影响因素，建立了一系列的经验公式。

公式中含有与温度和时间相关的试验常数，需要进行一系列的试验来确定。

虽然理论计算或计算机模拟可以在一定程度上确定钢中奥氏体的晶粒尺寸，但是具有很大的局限性。

而应用最广泛的经验公式具有不
可预测的缺点。

因此，开发出具有可预测的计算公式或者使用广泛的模拟计算方法，是今后奥氏体晶粒定量化研究的方向。

参考文献
【相关文献】
[1] 崔忠圻. 金属学与热处理[M]. 北京: 机械工业出版社, 2000.
[2] MOON J, LEE J B, LEE C H. Prediction for the austenite grain size in the presence of growing particles in the weld HAZ of Ti- microalloyed steel[J]. Materials Science and Engineering A, 2007, 459(1/2):40- 46.
[3] AVRAMI M. Kinetics of phase change. I General theory[J]. Journal of Chemical Physics, 1939, 7(12):1103- 1112.
[4] MISHRA S, DEBROY T. Non- isothermal grain growth in metals and alloys [J]. Materials Science and Technology, 2006, 22(3):253- 278.
[5] MANOHAR P A, FERRY M, CHANDRA T. Five decades of the Zener equation [J]. ISIJ International, 1998, 38(9): 913- 924.
[6] SMITH C S. Grain, phases and interface: An interpretation of microstructure [J]. Transactions of the Metallurgical Society of AIME, 1948, 175:15- 51.
[7] FELTHAM P. Grain growth in metals [J]. Acta Metallurgica, 1957, 5(2): 97- 105.
[8] 王达. 镁/锆处理对船板钢奥氏体晶粒长大行为的影响[D]. 沈阳:东北大学, 2014.
[9] BURKE J E, TURNBULL D. Recrystallization and grain growth[J]. Progress in Metal Physics, 1952, 3: 220- 292.
[10] SIMPSON C J, AUST K T, WINEGARD W C. The four stages of grain growth[J]. Metallurgical Transactions, 1971, 2(4): 987- 991.
[11] POUS- ROMERO H, LONARDELLI I, COGSWELL D, et al. Austenite grain growth in a nuclear pressure vessel steel[J]. Materials Science and Engineering A, 2013, 567:72- 79.
[12] 高振英, 孙福玉. 中碳钢奥氏体晶粒生长动力学的研究[J]. 钢铁研究学报, 1989, 1(4): 59- 63.
[13] HU H, RATH B B. On the time exponent in isothermal grain growth [J]. Metallurgical Transactions, 1970, 1(11): 3181- 3184.
[14] MEDINA S F, HERNANDEZ C A. General expression of the Zener- Hollomon parameter as a function of the chemical composition of low alloy and microalloyed steels [J]. Acta Materialia, 1996, 44(1):137- 148.
[15] 黄英, 张珉. 不同因素对C- Mn钢晶粒长大趋势的影响[C]//中国金属学会轧钢学会钢管学术委
员会五届三次年会. 济南:中国金属学会, 2007:391- 396.
[16] LEE S J, LEE Y K. Prediction of austenite grain growth during austenitization of low alloy steels [J]. Materials and Design, 2008, 29(9):1840- 1844.
[17] UHM S, MOON J, LEE C H, et al. Prediction model for the austenite grain size in the coarse grained heat affected zone of Fe- C- Mn steels: Considering the effect of initial grain size on isothermal growth behavior [J]. ISIJ International, 2004, 44(7):1230- 1237.
[18] DOHERTY R D, HUGHES D A, HUMPHREYS F J, et al. Current issues in recrystallization:
A review [J]. Materials Science and Engineering A, 1997, 238(2): 219- 274.
[19] MAALEKIAN M, RADIS R, MILITZER M, et al. In situ measurement and modelling of austenite grain growth in a Ti/Nb microalloyed steel [J]. Acta Materialia, 2012, 60(3): 1015- 1026.
[20] MA H, LIAO S L, WANG S F. Effect of Ti on austenite grain growth behavior in high carbon steels [J]. Journal of Iron and Steel Research, International, 2014, 21(7):702- 709. [21] QUISPE A, MEDINA S F, GóMEZ M, e t al. Influence of austenite grain size on recrystallisation- precipitation interaction in a V- microalloyed steel [J]. Materials Science and Engineering A, 2007, 447(1/2):11- 18.
[22] 胡晓宏, 刘雅政, 张永青, 等. Nb微合金化HRB400钢筋加热温度的研究[J]. 首钢科技,
2008(2):10- 13,17.
[23] KENNETT S C, KRAUSS G, FINDLEY K O. Prior austenite grain size and tempering effects on the dislocation density of low- C Nb- Ti microalloyed lath martensite [J]. Scripta Materialia, 2015, 107:123- 126.
[24] SELLARS C M, WHITEMAN J A. Recrystallization and grain growth in hot rolling [J]. Metal Science, 1979, 13(3/4):187- 194.
[25] ANELLI E. Application of mathematical modeling to hot rolling and controlled cooling of wire rods and bars [J]. ISIJ International, 1992, 32(3):440- 449.
[26] 何宜柱, 陈大宏, 雷廷权, 等. 形变Z因子与动态再结晶晶粒尺寸间的理论模型[J]. 钢铁研究学报, 2000, 12(1):26- 30.
[27] DEHGHAN- MANSHADI A, BARNETT M R, HODGSON P D. Hot deformation and recrystallization of austenitic stainless steel: Part I. Dynamic recrystallization [J]. Metallurgical and Materials Transactions A, 2008, 39(6):1359- 1371.
[28] 党淑娥, 宿展宁, 刘志龙, 等. 30Cr2Ni4MoV钢铸态加热过程中奥氏体晶粒的长大行为[J]. 材料热处理学报, 2014, 28(9): 675- 683.
[29] 何智勇, 贺晓军, 刘洋, 等. GH3690合金冷轧热处理过程晶粒长大模型[C]//中国热处理技术路线图高层论坛暨第八届中国热处理活动周论文集. 厦门:中国热处理学会,2012:188- 192.
[30] JUNG K H, LEE H W, IM Y T. Numerical prediction of austenite grain size in a bar rolling process using an evolution model based on a hot compression test [J]. Materials Science and Engineering A, 2009, 519 (1/2):94- 104.
[31] HODGSO P D, GIBBS D R K. A mathematical model to predict the mechanical properties of hot rolled C- Mn and microalloyed steels [J]. ISIJ International, 1992, 32(12):1329- 1338.。