粉末高温合金VIDP熔炼过程中碳脱氧动力学的计算方法技术

本发明专利技术涉及金属加工技术领域，尤其是涉及一种粉末高温合金VIDP熔炼过程中碳脱氧动力学的计算方法。包括如下步骤：S1：建立待求解的VIDP炉的几何模型；S2：赋予材料物性参数；S3：添加物理场并设置初始值和边界条件；包括磁场、湍流、化学和稀物质传递物理场；通过洛伦兹力将磁场与湍流耦合；S4：对几何模型进行网格划分；S5：求解计算得到VIDP熔炼的脱氧过程；S3中，在稀物质传递物理场中，将坩埚壁面及底面设置为通量边界，熔池上表面为反应边界，并设置C和O的扩散系数。本发明专利技术的计算方法得到的熔池平均氧含量随时间的变化曲线与实验测定值较为接近，能较好地预测VIDP熔炼过程中碳脱氧的实际情况，具有重要意义。具有重要意义。具有重要意义。

【技术实现步骤摘要】
粉末高温合金VIDP熔炼过程中碳脱氧动力学的计算方法


[0001]本专利技术涉及金属加工
，尤其是涉及一种粉末高温合金VIDP熔炼过程中碳脱氧动力学的计算方法。

技术介绍

[0002]粉末高温合金是制造先进航空发动机涡轮盘的关键材料，随着现代航空发动机推重比的提高，人们对粉末高温合金的质量要求也愈加严格。控制合金中杂质元素含量对提升材料质量和性能有巨大影响。氧作为粉末高温合金中的主要杂质元素之一，会与合金元素形成氧化物，成为零件服役过程中的裂纹源。若要得到性能优异的粉末高温合金材料，必须将氧含量控制在极低的水平。真空感应脱气浇铸熔炼(VIDP)是冶炼粉末高温合金母合金的最基本工序。该工艺利用电磁感应在炉料中产生涡流使其加热和熔化，并通过真空脱氧、脱氮、杂质元素挥发以及控制熔体和坩埚作用等一系列物理化学反应，冶炼出化学成分准确且纯净度高的粉末高温合金母合金。
[0003]在VIDP熔炼过程中，由于气相压力低，碳具有很强的脱氧能力，碳氧反应生成的CO气体能够不断地被抽走，进而促进了反应的进行，达到显著降低合金中氧含量的目的。此外，熔炼过程中强电磁搅拌作用导致熔池流动，加速了熔池内的传质速率，具有较好的碳氧反应动力学条件。VIDP熔炼条件下，粉末高温合金母合金碳脱氧过程包括下列几个步骤：(1)溶解于合金液中的[C]和[O]原子通过对流和扩散迁移到合金液边界层；(2)[C]和[O]原子通过液体边界层内的对流或扩散传递到液
‑
气界面；(3)[C]和[O]原子通过[C]+[O]＝CO在液
‑
气界面的化学反应吸附形成分子；(4)CO通过气体边界层的扩散到气体边界层；(5)将CO从气体边界层转移到真空中，最后由真空系统去除。一般而言，化学反应速率很快，气体产物的扩散也很快，所以，碳氧反应的动力学限制性环节通常表现为[C]、[O]原子向反应界面的扩散。在粉末高温合金中，[C]含量远高于[O]。另外由于[C]的扩散系数比[O]大，所以碳氧反应的限制性环节为[O]扩散。
[0004]在VIDP熔炼过程中，在高温、高真空度条件下合金液中的总氧含量的影响因素众多。目前，国内外学者研究VIDP熔炼过程中氧含量变化时只简单地考虑了碳脱氧反应的影响，而没有考虑其它因素对氧含量分布的影响，导致结果误差较大，无法准确得出熔池中氧含量的分布。
[0005]有鉴于此，特提出本专利技术。

技术实现思路

[0006]本专利技术的目的在于提供粉末高温合金VIDP熔炼过程中碳脱氧动力学的计算方法，以解决现有技术中存在的VIDP熔炼过程中碳脱氧动力学计算结果与实际情况存在较大差异等的技术问题。
[0007]为了实现本专利技术的上述目的，特采用以下技术方案：
[0008]粉末高温合金VIDP熔炼过程中碳脱氧动力学的计算方法，包括如下步骤：
[0009]S1：利用Comsol Multiphysics建立待求解的真空感应脱气浇铸炉的几何模型；
[0010]S2：赋予所述几何模型中材料的物性参数；
[0011]S3：在所述几何模型中添加物理场并设置所述几何模型的初始值和边界条件；其中，所述物理场包括磁场、湍流、化学和稀物质传递物理场；并通过洛伦兹力将磁场与湍流耦合；
[0012]S4：对所述几何模型进行网格划分；
[0013]S5：求解计算得到VIDP熔炼过程中的脱氧过程；
[0014]其中，S3中，在稀物质传递物理场中，将坩埚壁面及底面设置为通量边界，熔池上表面为反应边界，并设置C和O的扩散系数。
[0015]在VIDP熔炼过程中，由于粉末高温合金母合金中实际溶解氧含量很低，高温、高真空度条件下极易造成炉衬分解并向熔池供氧，导致熔池中氧含量升高。因此，碳氧反应脱氧和坩埚侵蚀供氧两方面共同作用从而影响合金液中的总氧含量。
[0016]本专利技术在计算VIDP熔炼过程中碳氧反应脱氧的同时，考虑了坩埚供氧对熔池中氧含量的影响，所得结果将为改进粉末高温合金母合金VIDP熔炼过程提供有效的依据。
[0017]在本专利技术的具体实施方式中，S1中，所述几何模型的的计算域包括合金液计算域、坩埚计算域、保温材料计算域、感应线圈计算域、磁轭计算域、保护气体计算域和冷却水计算域。
[0018]在实际操作中，根据真空感应脱气浇铸炉熔炼过程以及实验装置的实物图，来建立待求解的真空感应脱气浇铸炉的几何模型。
[0019]在本专利技术的具体实施方式中，S2中：材料包括合金液、坩埚、保温材料、感应线圈、磁轭、保护气体和冷却水；物性参数包括：导热系数、恒压热容、电导率、密度、相对磁导率、相对介电常数、动力黏度和比热率。
[0020]在本专利技术的具体实施方式中，所述几何模型中，所有计算域满足安培定律。进一步的，在合金液计算域中添加力计算。
[0021]在本专利技术的具体实施方式中，所述几何模型为二维轴对称模型。进一步的，S3中，所述磁场的设置中：设置左边界为轴对称边界，假设炉壁边缘为无限元域，并设置炉壁边缘为磁绝缘；线圈采用单导线线圈组，线圈激励方式为电压或电流。
[0022]在本专利技术的具体实施方式中，S3中，所述湍流的设置中：采用k
‑
ε湍流模型对真空感应脱气浇铸炉中流体的流动进行分析，将所有流体设置为不可压缩流体，设置坩埚壁边界条件为无滑移，并用壁函数对坩埚壁进行壁处理；设置合金液速度场的初始值为0，受重力和体积力，并将磁场中计算得到的洛伦兹力作为计算结果输入体积力中；采用动网格模拟保护气体与合金液接触面随时间的变化情况。
[0023]在本专利技术的具体实施方式中，S3中，所述化学的设置中：在所述化学物理场下添加本体物质C和O、表面物质CO以及溶剂合金液，并输入相应的摩尔质量；添加碳
‑
氧反应，并设置碳
‑
氧反应的反应速率、正反应频率因子以及正反应活化能参数。
[0024]在本专利技术的具体实施方式中，S3中，所述稀物质传递的设置中：假设熔池内溶液流动速度不受O传质影响。
[0025]在本专利技术的具体实施方式中，S4中，网格类型为三角形网格。进一步的，采用多区域多尺度划分的方式进行所述网格划分。
[0026]在本专利技术的具体实施方式中，S5中，求解计算得到VIDP熔炼过程中，熔池中氧含量分布。进一步的，S5中，添加频域
‑
瞬态研究来计算VIDP熔炼过程中的湍流流场；在研究中设置时间步、频率、容差，在物理场接口中选择磁场、湍流、动网格，在多物理场中选择洛伦兹力将磁场和湍流流场进行耦合；瞬态求解器选择MUMPS，计算得到湍流速度场分布；
[0027]对所述几何模型重新进行网格划分；
[0028]添加瞬态研究来计算熔炼过程中的脱氧行为；在研究中设置时间步，在物理场接口中选择化学和稀物质传递；瞬态求解器选择MUMPS，计算得到熔池中氧含量分布。
[0029]在本专利技术的具体实施方式中，还包括结果分析。进一步的，所述结果分析包括：分别绘制不同熔炼时间的VI本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.粉末高温合金VIDP熔炼过程中碳脱氧动力学的计算方法，其特征在于，包括如下步骤：S1：利用Comsol Multiphysics建立待求解的真空感应脱气浇铸炉的几何模型；S2：赋予所述几何模型中材料的物性参数；S3：在所述几何模型中添加物理场并设置所述几何模型的初始值和边界条件；其中，所述物理场包括磁场、湍流、化学和稀物质传递物理场；并通过洛伦兹力将磁场与湍流耦合；S4：对所述几何模型进行网格划分；S5：求解计算得到VIDP熔炼过程中的脱氧过程；其中，S3中，在稀物质传递物理场中，将坩埚壁面及底面设置为通量边界，熔池上表面为反应边界，并设置C和O的扩散系数。2.根据权利要求1所述的方法，其特征在于，S1中，所述几何模型的的计算域包括合金液计算域、坩埚计算域、保温材料计算域、感应线圈计算域、磁轭计算域、保护气体计算域和冷却水计算域。3.根据权利要求1所述的方法，其特征在于，S2中：材料包括合金液、坩埚、保温材料、感应线圈、磁轭、保护气体和冷却水；物性参数包括：导热系数、恒压热容、电导率、密度、相对磁导率、相对介电常数、动力黏度和比热率。4.根据权利要求1所述的方法，其特征在于，所述几何模型中，所有计算域满足安培定律；优选的，在合金液计算域中添加力计算。5.根据权利要求1所述的方法，其特征在于，所述几何模型为二维轴对称模型；优选的，S3中，所述磁场的设置中：设置左边界为轴对称边界，假设炉壁边缘为无限元域，并设置炉壁边缘为磁绝缘；线圈采用单导线线圈组，线圈激励方式为电压或电流；优选的，S3中，所述湍流的设置中：采用k
‑
ε湍流模型对真空感应脱气浇铸炉中流体的流动进行分析，将所有流体设置为不可压缩流体，设置坩埚壁边界条件为无滑移，并用壁函数对坩埚壁进行壁处理；设置合金液速度场...

【专利技术属性】
技术研发人员：杨树峰，赵朋，张明，曲敬龙，刘明东，贾建，刘建涛，张雪良，田甜，
申请(专利权)人：北京科技大学，
类型：发明
国别省市：