分享：高钒合金型内熔化扩散温度场及钒含量分布模拟

余世浩1,李佳琪2,张琳琅1 

(1.武汉华夏理工学院智能制造学院,武汉 430223;2.武汉理工大学材料科学与工程学院,武汉 430070) 

摘 要:利用 ANSYS 软 件 建 立 V9Cr4 高 钒 合 金 型 内 熔 化 扩 散 有 限 元 模 型,模 拟 了 浇 铸 5CrNiMo合金液后凝固过程中高钒合金棒不同位置处的温度;基于合金棒轴向温度变化不大的模 拟结果,将合金元素的三维扩散简化成二维扩散,建立径向钒元素含量分布数学模型,计算了径向 钒元素含量并进行了试验验证。结果表明:浇铸时合金棒的温度高于其固相线温度,说明高钒合金 棒与5CrNiMo合金能实现冶金结合;由温度分布曲线确定钒元素的扩散时间为810s,将其代入钒 元素含量分布数学模型,计算得到的不同位置处钒元素含量与测试结果的相对误差小于1%,说明 建立的型内熔化扩散有限元模型较为准确,可以利用其模拟结果来计算钒元素含量分布。 

关键词:合金型内熔化扩散;温度场;元素扩散 中图分类号:TG27 文献标志码:A 文章编号:1000-3738(2021)05-0071-05

0 引 言 

在装备制造、冶金机械、矿山开采等工业领域, 许多设备零部件因使用工况非常恶劣,需兼具高强 韧性和高耐磨性以提高其使用安全性和使用寿命; 但是采用单一均质材料制造的零部件往往难以兼顾 这两种性能。合金型内熔化扩散方法[1]是一种固液 复合铸造方法,根据钒与碳反应形成高耐磨碳化钒 颗粒的原理,在模具型腔底部预置高钒合金棒,浇入 高温液态金属进行凝固成型,通过高温使高钒合金熔化并与金属液之间发生元素互扩散和钒碳反应, 从而在铸件表层生成呈梯度分布的碳化钒硬质相, 实现表层耐磨性和深层韧性的良好匹配。 目前,大多采用有限元方法对金属铸造过程进 行预测进而实现有效控制。已有研究大多集中在通 过有限元软件模拟常规金属铸造时的凝固温度场, 结合缩孔、疏松等缺陷的热力学条件,对铸件中的缩 孔、疏松以及热裂纹等缺陷进行预测[2-4]。合金型内 熔化工艺由于型腔中内置的合金棒与液态金属的理 化特性不同,凝固时的传质、传热过程也与常规浇铸 过程不同。合金元素的扩散是影响铸件综合性能的 关键因素,然而目前关于温度场对该工艺下合金元 素扩散规律的影响研究尚不多见。余世浩等[5-6]研 究了型内熔化扩散过程中合金元素的浓度分布,建 立了合金元素扩散方程。作者采用有限元方法模拟 了高钒合金型内熔化扩散温度场,根据温度模拟结 果建立了钒元素含量分布的数学模型,模拟了钒元 素沿铸件径向的含量分布曲线并进行了验证。

 1 温度场计算模型

 1.1 基本假设 

合金型内熔化扩散铸造方法如图1所示。在型 腔内放置高钒合金棒会导致金属液的凝固过程更加 复杂,若在求解时考虑所有影响因素,则会造成计算 过于复杂。因此,对凝固过程作如下假设[7]:液态金 属在充满型腔的瞬间开始凝固;凝固时液态金属没 有相对流动,即只考虑热传导作用;砂型温度变化较 小,其热物性参数为与温度无关的常数;不考虑过冷 度,即低于液相线温度后,液态金属开始凝固,冷却 到固相线温度时,凝固过程结束。 图1 合金型内熔化扩散铸造方法示意 Fig.1 Diagramofalloyin-moldmeltingdiffusionmethod

 1.2 初始及边界条件

 初始条件 指 时 间 为 0 时 型 腔 内 的 温 度 分 布, 包括高温金属液 在 充 满 型 腔 时 的 初 始 温 度(试 验 温度为1550℃),合金棒与砂型的初始温度(取室 温25 ℃)。 基于整个铸造系统进行温度场计算,故边界条 件主要指温度场的表面换热条件。金属液充满型腔 后,其凝固过程实际上也是砂型向外逐渐散热的过 程。在模拟时,假设铸件(扩散件)和砂型之间不存 在间隙,则整个铸造系统的传热过程主要包括高温 金属液、合金棒和砂型之间的传热,以及砂型与外部 环境之间的对流和辐射传热。金属液浇注进型腔 后,砂型表面升温不明显,热辐射较小,并且辐射传 热具有显著的非线性特征。因此,为减少计算量、提 高分析效率,将辐射传热造成的影响换算成对流散 热,从而简化边界条件[8]。

 1.3 热物性参数 

在型内熔化扩散过程中,合金棒的表面温度应 高于其固相线温度,从而使界面处发生冶金反应实 现冶金结合。因此在金属液浇铸凝固过程中,合金 棒材料会发生相变,需要考虑相变潜热对温度场的 影响。 ANSYS软件可以通过定义材料参数来计算相 变潜热,计算公式为 H =∫ρc(T)dT (1) 式中:H 为相变潜热;ρ 为密度;c(T)为比热容;T 为热力学温度。 在进行相变分析时,将相变潜热作为材料特性 输入模型[9]。砂型的热物性参数与温度无关[10],其 导热系数为0.52 W·m -1·K -1,密度为 1630kg· m -3,比热容为1220J·kg -1·K -1。V9Cr4高钒合 金和5CrNiMo合 金 的 热 物 性 参 数 分 别 见 表 1 和 表2。 表1 V9Cr4高钒合金的热物性参数 Table1 ThermophysicalparametersofV9Cr4highvanadiumalloy 温度/ ℃ 导热系数/ (W·m -1·K -1) 比热容/ (J·kg -1·K -1) 密度/ (kg·m -3) 焓/ (J·m -3) 0 18.00 137.00 0 1100 29.00 475.00 3.06×10 9 1147 29.47 869.73 7800 3.30×10 9 1164 29.64 547.15 3.40×10 9 1200 30.00 595.11 3.56×10 9 

2 凝固过程温度场有限元模拟及结果

 2.1 有限元模型

 整 个铸造体系结构为对称结构,取其1/4进行 72 余世浩,等:高钒合金型内熔化扩散温度场及钒含量分布模拟 表2 5CrNiMo合金的热物性参数 Table2 Thermophysicalparametersof5CrNiMoalloy 温度/ ℃ 导热系数/ (W·m -1·K -1) 比热容/ (J·kg -1·K -1) 焓/ (J·m -3) 0 62.5 468 0 1452 65.2 632 8.24×10 9 1502 58.2 619 1.05×10 10 1532 58.2 539 1.40×10 10 研究,如图2所示。利用 ANSYS软件实体建模工 具构建三维实体模型,选用 Solid70单元类型进行 网格划分。由于合金棒和砂型的尺寸相差很大,若 采用自由网格划分,会导致网格畸变,影响分析精 度。因此,采用映射网格划分,结果如图3所示。 图2 1/4铸造体系结构的几何图 Fig.2 Geometryof1/4castingsystem 图3 铸件有限元模型及网格划分 Fig.3 Finiteelementmodelofcastingandmeshgeneration

 2.2 加载求解 

分别设定砂型、合金棒、液态金属的初始温度, 并在砂型外表面施加对流载荷。设定对流换热系数 为65 W·m -2·K -1,中心截面处为绝热边界。首先 进行时间步长为0.01s的稳态分析,将分析得到的 温度场认定为整个瞬态分析的初始温度,然后指定 分析类型为瞬态求解,载荷为阶跃施加,将自动时间 步长设置为开,终止时间设置为1000s,研究该段 时间内的温度场变化规律。

 2.3 温度场模拟结果分析 

为分析合金型内熔化扩散过程的温度场变化, 取图4中 A,B,C 三点进行分析。其中,A,C 为合 金棒上的点,A 点为合金棒端面和金属液接触的中 心点,C点为与砂型端面平行的合金棒的中心点,B 点为砂型端面与合金棒接触的金属液上的点。 图4 扩散件结构和尺寸示意 Fig.4 Diagramofdiffuserstructureandsize 图5 模拟得到图4中不同点的温度随时间的变化曲线 Fig.5 Simulatedcurvesoftemperaturevstimeatdifferentpoints showninFig 4 a pointsAandCand b pointsBandC 如图5(a)所示,金属液浇入砂型后,通过热传 导使合金棒温度快速上升,合金棒轴向 A 点和 C点 先后达到最高温度,而后随铸件一起冷却。A 点直 接接触高温金属液,其达到的峰值温度高于 C 点。 然而由于合金棒的尺寸不大,热传导很快,冷却过程 中 A 点和 C点的温差相对较小(约30 ℃)。因此可 认为合金棒周围的温度沿轴向变化不大,将合金棒 径向和轴向的三维扩散简化成径向的二维扩散进行 分析。由图5(b)可以发现,C 点的温度低于 B 点, 这是由于 C点不与高温金属液直接接触,且合金棒 的导热系数比砂型大造成的,这与实际情况相符。 液-固双金属要实现稳固的冶金结合,界面温度 要达到固态金属固相线以上20~40 ℃。V9Cr4高 73 余世浩,等:高钒合金型内熔化扩散温度场及钒含量分布模拟 钒合金的化学成分(质量分数)为2.5%C,0.5%Si, 0.8%Mn,8.9%V,4.2%Cr,根据合金 熔 点 计 算 公 式[11-12]求得其熔点约为1340℃。从C点的温度变 化曲线可知,合金棒发生熔化,这说明该型内熔化扩 散工艺可以实现5CrNiMo合金和 V9Cr4高钒合金 的冶金结合。

 3 钒元素含量分布计算及测试 

3.1 钒元素含量计算方法 

高温金属液浇入型腔后,高钒合金棒表面发生 熔化,钒元素向金属液中扩散迁移。整个铸造体系 中各个位置点的钒元素含量随时间不断变化,为不 稳定扩散,可以用 Fick第二定律来描述,即: ∂C ∂t = ∂ ∂x D ∂C ∂x (2) 式中:C 为钒元素质量浓度;D 为扩散系数;t 为扩 散时间;x 为钒元素扩散深度。 由图5中的温度模拟结果可知,合金棒轴向温 度变化不大,即可认为合金棒周围温度场分布呈圆 柱状。因此,可将钒元素在径向和轴向的三维扩散 简化为径向的二维平面扩散。将合金棒截面圆心设 置为直角坐标系的原点,则扩散过程中钒元素沿径 向的含量与时间及空间位置的关系为 C(x,t)= 1 2 πDt∫ r -r C0exp - (x-ξ)2 4Dt ?? ?? ???? ?? ?? ???? dξ (3) 式中:C0 为基体金属毛细管内的钒元素质量浓度;r 为合金棒端面半径;ξ为扩散区中任一点的位置。

 3.2 钒元素含量分布计算结果及测试验证

 对于型内熔化过程中合金元素的径向扩散行 为,主要影响因素是扩散系数和时间。在高温状态 下,钒元素在5CrNiMo合金毛细管内的扩散系数为 3.580×10 -5~3.733×10 -5cm 2·s -1[6]。 一般来说,合金元素在液态金属中的扩散系数 比在固体金属中高2到3个数量级。因此,在型内 熔化过程中,可忽略固态金属中合金元素的扩散行 为,只考虑液态金属或者固液混合时合金元素的扩 散,即可将液态金属/固液混合时的扩散分布看作该 合金元素的最终含量分布,而扩散时间则为金属液 充满砂型后到凝固结束需要的时间。5CrNiMo合 金的固液转变温度为1399 ℃,根据 B 点的温度模 拟结果,如图6所示,可推知钒元素的扩散时间为 810s。 将 扩 散 系 数 和 扩 散 时 间 代 入 式 (2),利 用 MATLAB软件计算得到钒元素在型内熔化扩散过 图6 B点温度模拟结果 Fig.6 TemperaturesimulationatpointB 程中的含量分布。以合金棒横截面的圆心为原点, 利用JXA-8230型电子探针沿径向对铸件进行钒元 素含量检测,得到钒元素含量沿径向的分布曲线。 由图7可以看出,钒元素含量分布的计算结果 与测试结果基本吻合,不同位置钒元素含量的相对 误差小于1%,说明利用型内熔化温度模拟结果来 计算钒元素含量分布是可行的。建立的型内熔化扩 散有限元模型和钒元素含量分布数学模型,可以用 于型内熔化工艺优化和钒元素扩散控制,即选择合 适的型内熔化工艺参数,通过调整温度场来实现合 理的钒元素含量分布。 图7 铸件径向钒含量分布曲线的实测结果与计算结果 Fig.7 Measuredandcalculatedresultsofvanadiumcontent radialdistributioncurvesofcastings 

4 结 论 

(1)建立了高钒合金型内熔化有限元模型,模 拟结果表明合金棒温度沿轴向变化不大,可将其径 向和轴向合金元素的三维扩散简化为径向的二维扩 散;在浇铸金属液后,合金棒端面及内部温度均高于 其固相线温度,说明合金棒熔化并与5CrNiMo合金 实现冶金结合。 (2)由有限元模拟得到的温度分布曲线,可知 钒元素的扩散时间为810s;将由型内熔化温度模拟 结果确定的钒元素扩散时间代入钒元素含量分布数 学模型,计算得到的径向不同位置处钒元素含量与 测试结果的相对误差小于1%,说明建立的型内熔 74 余世浩,等:高钒合金型内熔化扩散温度场及钒含量分布模拟 化扩散有限元模型较为准确,可以利用其模拟结果 来计算钒元素含量分布。

来源：材料与测试网