摘要
建立了扩散过程控制下的奥氏体连续冷却转变模型, 用于描述先共析铁素体界面位置随温度的变化规律. 模型考虑了界面移动过程中可能出现的软碰撞情况以及冷速对界面奥氏体侧C浓度的影响. 运用该模型对不同初始C浓度、奥氏体晶粒尺寸和冷速下的奥氏体连续冷却转变过程进行模拟, 得到不同冷却条件下的先共析铁素体界面位置、界面奥氏体侧的C扩散长度以及C浓度分布随温度的变化规律. 对晶粒尺寸为17 μm的Fe-0.17C合金进行不同冷速下的转变过程模拟, 所得结果与文献吻合较好.
关键词:
奥氏体向先共析铁素体转变, 既包括溶质C从界面处向母相奥氏体长程扩散的过程, 也包括界面处由fcc结构的母相奥氏体向bcc结构的铁素体转变的界面反应过程[1]. 研究[2-10]表明, 铁素体界面位置随时间的移动规律可以反映奥氏体转变量的变化规律. 对于奥氏体等温转变过程, 目前有2种描述其界面移动的经典模型: 扩散控制模型和界面控制模型[2,11-14]. 扩散控制模型认为, 由于界面迁移引起的自由能损耗可以忽略不计, 界面移动速率仅由扩散过程决定[11,15-17]; 界面控制模型认为, 相变驱动力完全被界面迁移过程所消耗, 界面移动速率与驱动力成正比[12,18]. 与奥氏体等温转变不同, 在连续冷却过程中, 奥氏体和铁素体中的平衡C浓度随温度降低不断变化, 且需要考虑冷速对界面移动的影响. 能否借鉴奥氏体等温转变的经典模型, 研究连续冷却过程中的界面位置移动规律, 目前还未见相关文献报道.
图 1扩散过程控制的奥氏体连续冷却转变界面移动和C浓度分布示意图
Fig.1Schematic of the diffusion fields evolution during austenite transforming in continuous cooling process for the diffusion-controlled growth model (x0is the position of austenite/ferrite interface, L is the length of diffusion field, Cgand Caare the equilibrium carbon concentrations in austenite and ferrite, C0is the bulk carbon concentration, Cmis the carbon concentration at the center of austenite, X is the centre of austenite, T1~T7are temperatures)
在实际生产中, 奥氏体连续冷却过程更为普遍. 目前描述连续冷却条件下先共析铁素体转变的模型主要是结合等温JKMA方程和Scheil叠加原理建立的[19-22]. 由于JKMA方程中参数的选取与相变条件有关, 需要基于大量实验数据得到, 并且该方程不能揭示奥氏体转变的微观机理. 因此, 本工作根据奥氏体向先共析铁素体转变的C扩散理论, 考虑了界面移动过程中可能出现的软碰撞情况以及冷速对界面奥氏体侧C浓度的影响, 建立了奥氏体连续冷却条件下转变动力学模型, 并分析了初始C浓度、奥氏体晶粒尺寸和冷却速率对转变过程的影响.
1 模型建立
根据扩散控制理论, 连续冷却奥氏体转变过程中, 界面迁移和奥氏体侧C浓度分布随温度的变化如图 1所示. 图中, x0表示奥氏体/铁素体界面位置, L表示C扩散长度, Cg和Ca分别表示界面奥氏体侧和铁素体侧的C浓度, C0表示奥氏体初始C浓度, Cm表示软碰撞发生后奥氏体晶粒中心处的C浓度, X表示奥氏体晶粒中心的位置, T1~T7表示温度的变化. 奥氏体连续冷却转变过程中, 奥氏体和先共析铁素体相界面从奥氏体晶界位置向奥氏体晶粒中心移动, 转变初始阶段界面奥氏体侧C扩散长度没有接触到奥氏体晶粒中心, 奥氏体晶粒中心C浓度保持为初始C浓度. 随着转变过程的进行, 奥氏体侧C扩散长度接触到奥氏体晶粒中心, 进入软碰撞阶段, 此时的奥氏体晶粒中心C浓度发生改变.
1.1 界面奥氏体侧C浓度分布
界面奥氏体侧的C浓度分布C(x)可近似用抛物线形式描述[5], 即C浓度扩散曲线随界面位置x的变化关系满足下式:
式中, A1, A2和A3由边界条件确定.
当等温转变或冷速非常小时, 界面奥氏体侧C浓度可达到平衡浓度
式中,
1.2 未发生软碰撞阶段
在该阶段, 奥氏体晶粒中心处的C浓度保持为C0, C扩散长度终点处C浓度变化率为0, 表达式分别如下:
联立式(1)~(4), 可得到未发生软碰撞阶段的C浓度分布为:
由于铁素体平衡C浓度在冷却过程中的变化非常小, 因此认为界面铁素体侧的C浓度始终保持为平衡C浓度, 即
联立式(5)和(7)可得L由下式确定:
式中,
此外, 由于C在界面处的扩散通量为0, 因此有:
式中,
对于奥氏体连续冷却转变过程, 由式(9)得:
联立式(5), (8)和(10)可得界面位置随温度的变化关系式:
式中, T为奥氏体转变过程中的温度, Ts为奥氏体开始转变温度.
1.3 发生软碰撞阶段
在软碰撞阶段, 奥氏体晶粒中心处的C浓度不再保持为C0, 而是随温度变化而变化. 因此有:
同样, 联立式 (1), (2), (4)和(12)可以得到C浓度分布为:
C在重新分配过程中质量守恒:
联立式(13)和(14)可得:
结合式(9), 可得:
其中,
奥氏体连续冷却转变过程中, 奥氏体和铁素体的平衡C浓度随温度的降低不断变化. 为了研究连续冷却过程中的界面位置变化以及界面奥氏体侧C浓度分布规律, 首先需要得到不同温度下的平衡C浓度, 本工作通过相变热力学公式计算各温度下的平衡C浓度. C在奥氏体中的扩散系数采用?gren[23]给出的公式计算.
未发生软碰撞阶段的界面位置可以由式(11)得到; 软碰撞发生后的转变过程, 可以看作是温度为T, 时间为Dt的等温过程的叠加, 结合式(13), (15)和(16)可得到界面位置、界面奥氏体侧C浓度分布以及奥氏体晶粒中心C浓度.
2 模型验证
利用本工作提出的奥氏体连续冷却转变模型, 计算C浓度为0.17% (质量分数, 下同), 初始奥氏体晶粒尺寸为17 mm, 不同冷速下的铁素体体积分数(界面位置与奥氏体晶粒半径的比)随温度的变化规律, 并与Kumar等[24]的计算结果进行对比, 如图 2所示. 式中的系数
此外, 由界面奥氏体侧C浓度修正系数A的表达式可以看出, 随着冷速的增加, A不断减小, 说明冷速大时, 界面奥氏体侧C浓度达到平衡C浓度变得困难.
3 模拟结果及讨论
本工作仅研究先共析铁素体转变, 因此将研究的温度下限选为600 ℃[25]. 为了研究不同奥氏体晶粒尺寸、冷却速率和C含量对铁素体转变过程的影响, 分别模拟奥氏体晶粒尺寸为15和30 mm, 冷却速率为1和5 ℃/s, 初始C含量为0.1%, 0.3%和0.5%时的铁素体转变过程.
图 2不同冷速下的铁素体体积分数计算值和文献值[24]对比
Fig.2Comparison between calculation results and literature results[24]at different cooling rates
图 3为不同奥氏体晶粒尺寸的Fe-0.1C, Fe-0.3C和Fe-0.5C合金, 在冷速为1和5 ℃/s时的相界面位置随温度的变化曲线. 其中实线表示考虑了软碰撞对相变过程的影响, 虚线表示没有考虑软碰撞对相变过程的影响, 圆圈表示软碰撞发生. 可以看出, 连续冷却过程中的界面位置与初始C浓度、奥氏体晶粒尺寸以及冷却速率有关. 随着初始C浓度、奥氏体晶粒尺寸和冷速增大, 界面位置随温度的变化速率减小. 若不考虑软碰撞的影响, 界面位置随温度的变化速率始终很快, 相界面会很快达到奥氏体晶粒中心位置, 表明铁素体相变完成, 这与实际转变过程不符, 也会给计算结果带来较大误差.
图3相界面位置随温度变化曲线
Fig.3Interface position as a function of temperature (d—austenite grain size,j—cooling rate)(a) d=15 mm,j=1 ℃/s (b) d=15 mm,j=5 ℃/s (c) d=30 mm,j=1 ℃/s (d) d=30 mm,j=5 ℃/s
软碰撞发生后, 奥氏体晶粒中心C浓度和C扩散长度是描述该阶段铁素体转变的2个重要参数. 奥氏体晶粒尺寸为30 mm的Fe-0.3C合金, 在冷速为1 ℃/s时界面奥氏体侧C扩散长度和奥氏体晶粒中心C浓度随温度变化曲线如图 4所示. 软碰撞前, 扩散长度随着界面迁移不断增加, 奥氏体晶粒中心的C浓度不受扩散的影响, 保持为初始C浓度; 软碰撞后, 扩散长度开始减小, 奥氏体晶粒中心的C浓度不断升高. 因此, 可以通过扩散长度或奥氏体晶粒中心C浓度变化来判断软碰撞发生的时刻.
晶粒尺寸为15 mm的Fe-0.3C合金在不同冷却速率下的C扩散分布曲线如图 5所示. 界面铁素体侧的C浓度变化很小, 界面奥氏体侧的C浓度不断升高.
图 4奥氏体晶粒尺寸为30 mm的Fe-0.3C合金在冷速为1 ℃/s时的界面奥氏体侧C扩散长度和奥氏体晶粒中心C浓度随温度变化曲线
Fig.4Diffusion length and solute concentration at the center of the austenite phase as a function of temperature in Fe-0.3C alloy with austenite grain size of 30 mm and cooling rate of 1 ℃/s
3.1 初始C浓度的影响
相同的初始奥氏体晶粒尺寸和冷速下, 不同初始C浓度下界面位置和C扩散长度随时间的变化规律如图 6所示. 可见, 随着C浓度增加, 软碰撞发生的时间推迟, 但由于界面移动速率减小, 因此相变结束时, 界面向奥氏体中心推进的距离减小, 表明铁素体体积分数和铁素体晶粒尺寸都减小. 奥氏体晶粒尺寸小或冷速小时, C浓度对软碰撞前的扩散长度的影响很小, 如图 6a~c所示.
图 5晶粒尺寸为15 mm的Fe-0.3C合金在不同冷却速率下的C扩散分布曲线
Fig.5Carbon profile evolution in the austenite phase in Fe-0.3C alloy with austenite grain size of 15 mm underj=1 ℃/s (a) andj=5 ℃/s (b)
图 6界面位置和C扩散长度随时间变化曲线
Fig.6Interface position and diffusion length as a function of time(a) d=15 mm,j=1 ℃/s (b) d=15 mm,j=5 ℃/s (c) d=30 mm,j=1 ℃/s (d) d=30 mm,j=5 ℃/s
3.2 冷速的影响
奥氏体晶粒尺寸为15 mm的Fe-0.3C合金在冷速为1和5 ℃/s时的界面位置随温度的变化曲线如图 7所示. 可见, 冷速大时, 铁素体相变开始温度低, C扩散系数小, 界面位置随温度变化缓慢.
3.3 奥氏体晶粒尺寸的影响
相同的初始C浓度和冷速下, 不同初始奥氏体晶粒尺寸下的界面位置和C扩散长度随时间的变化规律如图 8所示. 可见, 初始奥氏体晶粒尺寸增加, 软碰撞发生的时间推迟, 且由于界面位置移动速率快, 所以相变结束时, 界面向奥氏体中心推进的距离增加.
3.4 C扩散长度与界面位置的关系
奥氏体等温转变时, 式(8)中的修正系数A为1, 且平衡C浓度保持不变, 因此, S保持不变, 软碰撞前的C扩散长度随界面位置线性增加. 奥氏体连续冷却转变时, 由于平衡C浓度不断变化, S不再保持为定值. 图 9为晶粒尺寸为30 mm 的Fe-0.1C, Fe-0.3C和Fe-0.5C合金在不同冷速下的C扩散长度随界面位置的变化曲线. 冷速小时, 可以近似认为软碰撞前的C扩散长度随界面位置线性增加, 如图 9a所示. 冷速大时, 在奥氏体相变初期, C扩散长度随界面位置近似呈线性变化, 随着转变过程的进行, C扩散长度随界面位置的变化逐渐偏离线性关系, 如图 9b所示. 尽管随着C浓度增加, 软碰撞发生的时间会推迟, 但发生软碰撞时界面移动的距离减小.
图 7奥氏体晶粒尺寸为15 mm的Fe-0.3C合金在冷速为1和5 ℃/s时的界面位置随温度的变化曲线
Fig.7Interface position as a function of temperature in Fe-0.3C alloy with austenite grain size of 15 mm and cooling rates of 1 and 5 ℃/s
图 8Fe-0.3C合金在冷速为1和5 ℃/s时的界面位置和C扩散长度随时间变化曲线
Fig.8Interface position and diffusion length as a function of time in Fe-0.3C alloy underj=1 ℃/s (a) andj=5 ℃/s (b)
图 9晶粒尺寸为30 μm 的Fe-0.1C, Fe-0.3C和Fe-0.5C合金在不同冷速下的C扩散长度随界面位置变化曲线
Fig.9Diffusion length as a function of the interface position in Fe-0.1C, Fe-0.3C and Fe-0.5C alloys with austenite grain size of 30 μm underj=1 ℃/s (a) andj=5 ℃/s (b)
3.5 界面位置随时间的变化规律
奥氏体等温转变过程中, 未发生软碰撞阶段的界面位置随时间t的关系可以用下式表示[5]:
其中,
由式(17)可知, 该阶段的界面位置与t1/2呈线性关系. 软碰撞发生后, 由于奥氏体晶粒中心处的C浓度不断升高, 界面位置与t1/2将不再呈线性关系.
奥氏体连续冷却转变过程中, 界面位置随t1/2的变化曲线如图 10所示, 圆圈表示软碰撞发生. 当奥氏体晶粒尺寸小或者冷速慢时, 可以近似认为软碰撞前的界面位置与t1/2呈线性关系, 如图 10a~c所示. 与等温转变不同, 即使软碰撞发生前, 界面位置与t1/2也并不是始终保持为线性关系, 如图 10d所示. 这是由于随着温度降低, C扩散系数和平衡C浓度不断变化导致的.
图 10Fe-0.1C, Fe-0.3C和Fe-0.5C合金在不同晶粒尺寸和冷速下的界面位置随t1/2变化曲线
Fig.10Interface position as a function of t1/2in Fe-0.1C, Fe-0.3C and Fe-0.5C alloys with different austenite grain sizes and cooling rates(a) d=15 mm,j=1 ℃/s (b) d=15 mm,j=5 ℃/s (c) d=30 mm,j=1 ℃/s (d) d=30 mm,j=5 ℃/s
4 结论
根据奥氏体向先共析铁素体转变的C扩散理论, 建立了奥氏体连续冷却条件下转变动力学模型, 并分析了初始C浓度、奥氏体晶粒尺寸和冷却速率对转变过程的影响. 初始C浓度、奥氏体晶粒尺寸和冷速是影响奥氏体连续冷却转变的3个重要因素, 只有在奥氏体晶粒尺寸小或冷速慢时, 才可以近似认为软碰撞前的界面位置与t1/2呈线性关系. 冷速快时, 奥氏体转变过程中界面奥氏体侧C浓度不能用平衡C浓度近似, 其修正系数与冷速相关. 该模型的计算结果可以较好地反映奥氏体连续冷却转变的动力学规律.