CN105945270B - 一种定量浇铸曲线的获取方法及装置 - Google Patents

一种定量浇铸曲线的获取方法及装置 Download PDF

Info

Publication number
CN105945270B
CN105945270B CN201610424666.9A CN201610424666A CN105945270B CN 105945270 B CN105945270 B CN 105945270B CN 201610424666 A CN201610424666 A CN 201610424666A CN 105945270 B CN105945270 B CN 105945270B
Authority
CN
China
Prior art keywords
casting
copper water
tundish
pouring package
curve
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201610424666.9A
Other languages
English (en)
Other versions
CN105945270A (zh
Inventor
蒋朝辉
刘昕
桂卫华
阳春华
谢永芳
周家达
李中美
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Central South University
Original Assignee
Central South University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Central South University filed Critical Central South University
Priority to CN201610424666.9A priority Critical patent/CN105945270B/zh
Publication of CN105945270A publication Critical patent/CN105945270A/zh
Application granted granted Critical
Publication of CN105945270B publication Critical patent/CN105945270B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • B—PERFORMING OPERATIONS; TRANSPORTING
    • B22—CASTING; POWDER METALLURGY
    • B22D—CASTING OF METALS; CASTING OF OTHER SUBSTANCES BY THE SAME PROCESSES OR DEVICES
    • B22D37/00—Controlling or regulating the pouring of molten metal from a casting melt-holding vessel

Landscapes

  • Engineering & Computer Science (AREA)
  • Mechanical Engineering (AREA)
  • Manufacture And Refinement Of Metals (AREA)
  • Continuous Casting (AREA)

Abstract

本发明提供一种定量浇铸曲线的获取方法及装置，其中，该方法包括:建立铜水流动过程的数学模型，以及建立浇铸工具对应的物理模型；对物理模型和数学模型进行离散化处理；确定定量浇铸过程的初始条件和边界条件；通过基于Fluent软件的有限体积算法、物理模型、初始条件及边界条件对数学模型进行求解，得到铜水流动过程中浇铸包和中间包内铜水流动的数据，并生成定量浇铸曲线。实现通过建立铜水流动过程的数学模型，准确描述了高温铜水在浇铸过程中的运动状态和温度变化，通过对该数学模型进行求解并对求解得到的数值进行仿真，得到浇铸精度很高的定量浇铸曲线，根据该定量浇铸曲线进行铜水的定量浇铸，大大提高了生产出的铜阳极板的质量。

Description

一种定量浇铸曲线的获取方法及装置
技术领域
本发明涉及铜冶炼及定量浇铸技术领域，尤其涉及一种定量浇铸曲线的获取方法及装置。
背景技术
目前，铜冶炼过程中，定量浇铸技术直接影响着铜阳极板质量的好坏，而定量浇铸曲线作为定量浇铸技术的核心，定量浇铸曲线设置的合理与否直接影响阳极板的质量和产量。
当前，在定量浇铸过程中，技术人员首先根据经验设定定量浇铸过程所需的浇铸包和中间包的参数，当设定好各个参数后，圆盘浇铸机根据设定的参数生成定量浇铸曲线，并根据该定量浇铸曲线控制铜水的定量浇铸过程，当定量浇铸曲线不满足浇铸精度要求时，从而生产出的铜阳极板的质量不符合要求，此时，技术人员需要修改定量浇铸过程所需的浇铸包和中间包的参数，以使生成的定量浇铸曲线符合满足浇铸精度要求，以使生产出的铜阳极板符合质量要求。由于定量浇铸现场受高温、高尘、浓烟等恶劣环境的干扰，且定量浇铸过程需要很多的人工干预，如此，导致根据设置的浇铸包和中间包的参数确定出的定量浇铸曲线的准确性很低，根据准确性很低的定量浇铸曲线生产出的阳极板的质量很低。
发明内容
鉴于上述问题，提出了本发明以便提供一种克服上述问题或者至少部分地解决上述问题的一种定量浇铸曲线的获取方法及装置。
依据本发明的一个方面，提供了一种定量浇铸曲线的获取方法，所述方法包括:
分析高温铜水在浇铸工具内的流动状态特征，建立铜水流动过程的数学模型；
获取浇铸工具的实际尺寸，等比例建立浇铸工具对应的物理模型；
对所述物理模型和所述数学模型进行离散化处理；
确定定量浇铸过程的初始条件和边界条件；
通过基于Fluent软件的有限体积算法、所述物理模型、所述初始条件及所述边界条件对所述数学模型进行求解，得到铜水流动过程中浇铸包和中间包内铜水流动的数据，所述Fluent软件是基于计算流体动力学方法的求解流动问题的通用软件；
根据所述浇铸包和中间包内铜水流动的数据生成定量浇铸曲线。
可选地，所述分析高温铜水在浇铸工具内的流动状态特征，建立铜水流动过程的数学模型，包括：
根据高温铜水在浇铸工具内的流动状态特征，确定流动状态特征对应的参数和约束条件；
根据所述参数和所述约束条件，确定浇铸工具的流动状态特征对应的连续性方程、动量方程、能量守恒方程、流体控制VOF方程、湍流模型k-ε方程和用户自定义UDF函数，所述VOF方程表示体积分数分布和自由表面运动，k-ε方程表示流体的湍流流动，UDF函数方程用于确定浇铸包与中间包的运动状态；
通过所述连续性方程、所述动量方程、所述能量守恒方程、所述VOF方程、所述k-ε方程和所述UDF函数实现铜水流动过程的三维仿真；
将三维仿真后得到的各个模型，确定为定量浇铸过程的数学模型。
可选地，所述对所述物理模型和所述数学模型进行离散化处理，包括：
通过Meshing对所述物理模型进行网格划分，得到空间网格区域所述Meshing为ANSYS软件中的划分网格功能单元；
通过控制方程离散技术将所述数学模型在所述空间网格区域上进行离散。
可选地，所述确定初始条件和边界条件，包括：
通过UDF函数设置浇铸工具的运动曲线；
通过动网格技术控制浇铸工具实现所述运动曲线的运动过程；
将所述通过动网格技术控制浇铸工具实现所述运动曲线的运动过程，作为浇铸过程的初始条件和边界条件。
可选地，所述通过UDF函数设置浇铸工具的运动曲线，包括：
根据所述浇铸工具的运动状态，通过Profile文件设置各个时间点中间包和浇铸包的运动速度，所述Profile文件是Fluent软件中一种指定边界运动的描述方式；
判断中间包速度、倾角与浇铸铜水重量的关系是否满足预设浇铸条件，以及分析浇铸包不同位置、速度与浇铸铜水重量的关系是否满足预设浇铸条件；若中间包以及浇铸包均满足预设浇铸条件，通过Profile文件生成浇铸过程的离散时间点与速度的对应关系；若否，不断增加时间点个数，修正中间包或浇铸包各个时间点的运动速度。
可选地，根据所述浇铸包和中间包内铜水流动的数据生成定量浇铸曲线，包括：
对所述浇铸包和中间包内的铜水流动的数据进行分析，得到铜水的体积分数和速度场；
根据所述体积分数和所述速度场，确定中间包和浇铸包中铜水液面流动情况以及铜阳极质量；
根据所述铜水液面流动情况及所述铜阳极质量判断浇铸精度是否满足预设浇铸精度；若是，根据所述浇铸包和中间包内铜水流动的数据确定出浇铸过程的离散时间点与速度的对应关系，根据所述离散时间点与速度的对应关系拟合浇铸包和中间包的运动曲线，将所述浇铸包和中间包的运动曲线确定为浇铸工具的定量浇铸曲线，若否，则修正UDF函数，重新确定边界条件。
依据本发明的另一个方面，提供了一种定量浇铸曲线的获取装置，所述装置包括:
第一建立模块，用于分析高温铜水在浇铸工具内的流动状态特征，建立铜水流动过程的数学模型；
第二建立模块，用于获取浇铸工具的实际尺寸，等比例建立浇铸工具对应的物理模型；
处理模块，用于对所述物理模型和所述数学模型进行离散化处理；
确定模块，用于确定定量浇铸过程的初始条件和边界条件；
计算模块，用于通过基于Fluent软件的有限体积算法、所述物理模型、所述初始条件及所述边界条件对所述数学模型进行求解，得到铜水流动过程中浇铸包和中间包内铜水流动的数据，所述Fluent软件是基于计算流体动力学方法的求解流动问题的通用软件；
生成模块，用于根据所述浇铸包和中间包内铜水流动的数据生成定量浇铸曲线。
可选地，所述确定模块包括：
设置单元，用于通过用户自定义UDF函数设置浇铸工具的运动曲线；
控制单元，通过动网格技术控制浇铸工具实现所述运动曲线的运动过程；
第一确定单元，用于将所述通过动网格技术控制浇铸工具实现所述运动曲线的运动过程，作为浇铸过程的初始条件和边界条件。
可选地，所述设置单元包括：
设置子单元，用于根据所述浇铸工具的运动状态，通过Profile文件设置各个时间点中间包和浇铸包的运动速度，所述Profile文件是Fluent软件中一种指定边界运动的描述方式；
判断子单元，用于判断中间包速度、倾角与浇铸铜水重量的关系是否满足预设浇铸条件，以及分析浇铸包不同位置、速度与浇铸铜水重量的关系是否满足预设浇铸条件；
生成子单元，若判断子单元判断中间包以及浇铸包均满足预设浇铸条件，用于通过Profile文件生成浇铸过程的离散时间点与速度的对应关系；
修正子单元，若判断子单元判断中间包或浇铸包不满足预设浇铸条件，用于不断增加时间点个数，修正中间包或浇铸包各个时间点的运动速度。
可选地，生成模块包括：
分析单元，用于对所述浇铸包和中间包内的铜水流动的数据进行分析，得到铜水的体积分数和速度场；
第二确定单元，用于根据所述体积分数和所述速度场，确定中间包和浇铸包中铜水液面流动情况以及铜阳极质量；
判断单元，用于根据所述铜水液面流动情况及所述铜阳极质量判断浇铸精度是否满足预设浇铸精度；
第三确定单元，若判断单元判断浇铸精度是否满足预设浇铸精度时，根据所述浇铸包和中间包内铜水流动的数据确定出浇铸过程的离散时间点与速度的对应关系；
拟合单元，用于根据所述离散时间点与速度的对应关系拟合浇铸包和中间包的运动曲线；
所述第三确定单元，用于将所述浇铸包和中间包的运动曲线确定为浇铸工具的定量浇铸曲线；
修正单元，用于判断单元判断浇铸精度不满足预设浇铸精度时，修正UDF函数，以使重新确定边界条件。本发明提供的定量浇铸曲线的获取方法及装置，通过建立铜水流动过程的数学模型，准确描述了高温铜水在浇铸过程中的运动状态和温度变化，通过对该数学模型进行求解并对求解得到的数值进行仿真，可以得到浇铸精度很高的定量浇铸曲线，根据该定量浇铸曲线进行铜水的定量浇铸，大大提高了生产出的铜阳极板的质量。
附图说明
通过阅读下文优选实施方式的详细描述，各种其他的优点和益处对于本领域普通技术人员将变得清楚明了。附图仅用于示出优选实施方式的目的，而并不认为是对本发明的限制。而且在整个附图中，用相同的参考符号表示相同的部件。在附图中：
图1是本发明一种实施方式的一种定量浇铸曲线的获取方法的流程图；
图2是本发明一种实施方式的确定定量浇铸过程的数学模型的流程图；
图3是本发明一种实施方式的设置浇铸工具的运动曲线的流程图；
图4是本发明一种实施方式的生成定量浇铸曲线的流程图；
图5是本发明一种实施方式的一种定量浇铸曲线的获取装置的示意图。
具体实施方式
下面结合附图和实施例对本发明的实施方式作进一步详细描述。以下实施例用于说明本发明，但不能用来限制本发明的范围。
图1是本发明实施方式的一种定量浇铸曲线的获取方法的流程图。参见图1，该方法可通过以下S101-S106步骤实现。
步骤S101:分析高温铜水在浇铸工具内的流动状态特征，建立铜水流动过程的数学模型。
上述分析高温铜水在浇铸工具内的流动状态特征，建立铜水流动过程的数学模型，具体可通过以下S1011-S1014步骤建立，如图2所示。
步骤S1011:根据高温铜水在浇铸工具内的流动状态特征，确定流动状态特征对应的参数和约束条件。
在本发明实施例中，首先需要根据高温铜水在浇铸工具内的流动状态特征，剖析铜水在浇铸工具内流动的基本物理过程。剖析的过程如下：
根据高温铜水在浇铸工具内的流动状态特征，还原炉的铜水，让经过溜槽流向中间包，存储并将铜水按一定量交替倾倒在两边定量浇铸包中，以便进行铜阳极板的定量浇铸；浇铸包承接中间包倾倒的一定量的铜水，存储并将铜水按一定量浇铸在铸模中，进行阳极板精确浇铸。浇铸过程通常被认为是伴随着热量散失的不可压缩粘性流体带有自由表面的非稳态流动。浇铸过程具体的问题为浇铸工具的主动运动导致的高温铜水的运动变化，只考虑浇铸工具对高温铜水的作用力，需要预先设定浇铸工具的运动轨迹，进而计算得出高温铜水的运动规律。
在本发明实施例中，上述流动状态对应的参数，主要是与流体流动过程、中间包与浇铸包运动过程等相关的已知因素，包括溜槽中铜水流动速度，铜水温度，铜水成分，铜水粘性系数，铜水密度、工作时铜水容量和中间包与浇铸包运动参数等。上述流动状态对应的约束条件为上述各个参数的取值范围。
步骤S1012:根据参数和约束条件，确定浇铸工具的流动状态特征对应的连续性方程、动量方程、能量守恒方程、流体控制VOF方程、湍流模型k-ε方程和UDF函数(UserDefined Function，用户自定义函数)，VOF方程表示体积分数分布和自由表面运动，k-ε方程表示流体的湍流流动，UDF函数方程用于确定浇铸包与中间包的运动状态。
在本发明实施例中，上述确定的质量守恒方程如下式(1)所示，上述确定的动量守恒方程如下式(2)所示，上述确定的能量守恒方程如下式(3)所示，对于不可压缩的高温铜水，上述确定的k-ε方程如下式(4)和(5)所示。
式(1)中为，ρ为流体密度(kg/m3)；x为位移(m)；i为三维直角坐标系下的坐标方向；ui是xi方向上的速度分量(m/s)；t为时间(s)。
式(2)中，t是时间(s)；x为位移(m)，i和j为三维直角坐标系下的坐标方向，ui、uj分别是xi、xj方向上的速度分量(m/s)；p是压强(Pa)；ρ是流体密度(kg/m3)；μeff是有效粘度，μeff为μ和μt之和；μ是分子动力学粘度(Pa·s)，μt是湍流动力粘度(Pa·s)。
式(3)中，ρ为流体密度(kg/m3)；x为位移(m)；i为三维直角坐标系下的坐标方向；ui是xi方向上的速度分量(m/s)；t为时间(s)；keff为导热率(w/m·k)，H为显焓(J/mol)。
湍流动能方程
湍流耗散率方程
式(4)和式(5)中，G是由于平均速度梯度引起的湍动能k的产生项，G用等式可以表述为：
式(4)和(5)中，t是时间(s)；ui分别是xi方向上的速度分量(m/s)；ρ是流体密度(kg/m3)；μ是流体的分子动力学粘度(Pa·s)；k是湍动能(m2/s2)；ε是湍动耗散率(m2/s3)；p是压力(Pa)；σk和σε分别是k和ε湍流Prandtl数；C1ε，C2ε，σk和σε为模型常数，在标准k-ε模型中，根据Launder等的推荐值及后来的实验验证，C1ε的值为1.44，C2ε的值为1.92，σk的值为1.0，σε的值为1.33。
步骤S1013:通过连续性方程、动量方程、能量守恒方程、VOF方程、k-ε方程和UDF函数实现铜水流动过程的三维仿真。
上述步骤S1014:将三维仿真后得到的各个模型，确定为定量浇铸过程的数学模型。
要实现铜水流动过程的三维仿真，需将整个流动过程描述为若干基本过程，根据研究问题的性质和要求的精度，分别建立起基本方程中各变量之间的数学式，即上述连续性方程、动量方程、能量守恒方程、VOF方程、k-ε方程和UDF函数。
将上述连续性方程、动量方程、能量守恒方程、VOF方程、k-ε方程和UDF函数综合起来，得到浇铸过程的方程组，将该方程组及该方程组对应的约束条件组成的模型作为定量浇铸过程的数学模型。
在本发明实施例中，建立浇铸过程的数学模型的同时，需要通过以下S102步骤建立浇铸工具对应的物理模型。
步骤S102:获取浇铸工具的实际尺寸，等比例建立浇铸工具对应的物理模型。
上述浇铸工具的实际尺寸为生产现场提供的浇铸工具详细的几何尺寸，可以通过Ansys中DM等比例建立浇铸工具对应的实体模型，该实体模型即为浇铸工具对应的物理模型，其中，Ansys软件是一种有限元分析软件，DM是该Ansys软件中的建模功能单元。
步骤S103:对物理模型和数学模型进行离散化处理。
在本发明实施例中，上述对物理模型和数学模型进行离散化处理时，首先通过Ansys中的Meshing对物理模型进行网格划分，得到空间网格区域；然后通过控制方程离散技术将数学模型在空间网格区域上进行离散，其中，Meshing为Ansys软件的划分网格功能单元。
在本发明实施例中，对数学模型进行离散后，该数学模型中方程组中的各个方程得到简化，从而可以快速、准确地对简化后的各个方程进行求解，大大提高了模型的求解速度以及求解得到的数据的准确性。
步骤S104:确定定量浇铸过程的初始条件和边界条件。
在进行瞬态问题数值模拟时，需要设置初始条件。初始条件是所研究对象在过程开始时刻各个求解变量的空间分布情况，是确定偏微分方程具有唯一解的初值。
在本发明实施例中，在设置定量浇铸过程的初始条件时，需要设置初始时刻整个计算域内的压力场、气-液两相体积分数及高温铜水的温度和初始流动速度，同时由于流场和温度场的耦合计算，还需要设定初始状态下计算域内的温度场等。在Fluent软件中初始条件的设定通过计算域初始化功能选项完成。
边界条件是指在求解域边界上所求解的变量或其一阶导数随地点及时间变化的规律。
在本发明实施例中，在设定定量浇铸过程的边界条件时，除了设置基本边界条件，如进出口边界条件、给定压力边界和wall壁面边界条件等。还要根据浇铸的目标，使用UDF函数设置浇铸工具运动曲线，用动网格技术实现浇铸工具的运动。具体的，可通过预先设定的UDF(User Defined Function，用户自定义函数)函数设置浇铸工具的运动曲线；通过动网格技术控制浇铸工具实现运动曲线的运动过程；将通过动网格技术控制浇铸工具实现运动曲线的运动过程，作为浇铸过程的初始条件和边界条件。其中，通过UDF函数设置浇铸工具的运动曲线，具体可通过以下S1041-S1044步骤设置，如图3所示。
在设置浇铸工具的运动曲线时，需要首先分析浇铸工具的运动状态。浇铸工具包括中间包和浇铸包。中间包与浇铸包的运动状态不同，浇铸的精度也不相同。中间包的运动相对比较简单，主要是旋转运动，有最大旋转角度、最大旋转速度和浇铸时间的限制。浇铸包的运动十分复杂，为后滑移运动。由于浇铸包浇铸的精度要求很高，浇铸包的运动一般分解为五个运动阶段，快速启动浇铸包-慢速浇铸-快速浇铸-慢速浇铸-恒速浇铸-快速收起浇铸包。浇铸包的运动也有极限位、运动速度和浇铸时间的限制。
步骤S1041:根据浇铸工具的运动状态，通过Profile文件设置各个时间点中间包和浇铸包的运动速度。
在本发明实施例中，上述设置各个时间点的中间包的运动速度时，包括设置中间包各个时间点内的y方向速度、z方向速度和旋转角速度速度；上述设置各个时间点的浇铸包的运动速度时，包括设置浇铸包的旋转角速度、y方向速度和z方向速度、位置变化与浇铸铜水重量的关系。
步骤S1042:判断中间包速度、倾角与浇铸铜水重量的关系是否满足预设浇铸条件，以及分析浇铸包不同位置、速度与浇铸铜水重量的关系是否满足预设浇铸条件；若中间包以及浇铸包均满足预设浇铸条件，则执行以下S1043步骤，若中间包或浇铸包不满足预设浇铸条件，则执行以下S1044步骤。
步骤S1043:通过Profile文件生成浇铸过程的离散时间点与速度的对应关系。
步骤S1044:不断增加时间点个数，修正中间包或浇铸包各个时间点的运动速度。
在本发明实施例中，当修正中间包或浇铸包各个时间点的运动速度后，需要通过以上S1042步骤继续进行判断，直至满足预设浇铸条件。
步骤S105:通过基于Fluent软件的有限体积算法、物理模型、初始条件及边界条件对数学模型进行求解，得到铜水流动过程中浇铸包和中间包内铜水流动的数据，该Fluent软件是基于计算流体动力学方法的求解流动问题的通用软件。
步骤S106:根据浇铸包和中间包内铜水流动的数据生成定量浇铸曲线。
在本发明实施例中，上述根据浇铸包和中间包内铜水流动的数据生成定量浇铸曲线，具体可通过以下S1061-S1067步骤生成，如图4所示。
步骤S1061:对浇铸包和中间包内的铜水流动的数据进行分析，得到铜水的体积分数和速度场；
在本发明实施例中，对浇铸包和中间包内的铜水流动的数据进行分析时，首先将该浇铸包和中间包内的铜水流动的数据转化为可视化的浇铸过程图像，从该可视化的浇铸过程图像中获取铜水在浇铸工具内的速度场和温度场，根据温度场确定出铜水的体积分数。
步骤S1062:根据体积分数和速度场，确定中间包和浇铸包中铜水液面流动状态以及铜阳极质量；
在本发明实施例中，上述铜水液面流动状态可以是铜水流动速度。其中，可以通过速度场和铜水的体积分数，分析浇铸工具中不同容量的铜水、不同的运动速度、不同角度下的铜水流动速度。根据该铜水流动速度与铜阳极板的质量的对应关系，确定出铜阳极板的质量。
步骤S1063:根据铜水液面流动状态及铜阳极质量判断浇铸精度是否满足预设浇铸精度；若浇铸精度是否满足预设浇铸精度时，则执行以下S1064-S1066步骤，若浇铸精度是不满足预设浇铸精度，则执行以下S1067步骤。
在本发明实施例中，上述预设浇铸精度为生产出的铜阳极板质量符合质量要求时所需的浇铸精度，具体的浇铸精度值在浇铸前设定。
步骤S1064:根据浇铸包和中间包内铜水流动的数据确定出浇铸过程的离散时间点与速度的对应关系。
步骤S1065:根据离散时间点与速度的对应关系拟合浇铸包和中间包的运动曲线。
步骤S1066:将浇铸包和中间包的运动曲线确定为浇铸工具的定量浇铸曲线。
步骤S1067:修正UDF函数，重新确定边界条件。
在本发明实施例中，当需要重新确定边界条件时，具体可通过上述S104步骤提供的方法确定。
在本发明实施例中，修正UDF函数时，首先需要将上述拟合得到浇铸包和中间包的运动曲线用函数形式表示，编写成UDF函数。然后指定该运动曲线上特定动态区域的运动，产生随时间变化的浇铸工具运动的线速度、角速度及重心运动。最后，利用这些速度更新动态区域的位置。当重新确定边界条件后，将得到的运动曲线函数加载到Fluent软件中，通过该重新确定的边界条件以及上述S105步骤，以及S106步骤中的S1061-S1063步骤重新确定浇铸精度是否满足预设浇铸精度，若不满足再增加时间点的个数，然后再根据上述拟合得到的浇铸工具的运动曲线，重新确定UDF函数，直至浇铸精度满足预设浇铸精度时，通过S1064-S1066步骤生成定量浇铸曲线。
本发明实施提供的一种定量浇铸曲线的获取方法中，通过建立铜水流动过程的数学模型，准确描述了高温铜水在浇铸过程中的运动状态和温度变化，通过对该数学模型进行求解并对求解得到的数值进行仿真，可以得到浇铸精度很高的定量浇铸曲线，根据该定量浇铸曲线进行铜水的定量浇铸，大大提高了生产出的铜阳极板的质量。
图5是本发明实施方式的一种定量浇铸曲线的获取装置的示意图，参见图5，该装置包括:
第一建立模块S1，用于分析高温铜水在浇铸工具内的流动状态特征，建立铜水流动过程的数学模型；
第二建立模块S2，用于获取浇铸工具的实际尺寸，等比例建立浇铸工具对应的物理模型；
处理模块S3，用于对物理模型和数学模型进行离散化处理；
确定模块S4，用于确定定量浇铸过程的初始条件和边界条件；
计算模块S5，用于通过基于Fluent软件的有限体积算法、物理模型、初始条件及边界条件对数学模型进行求解，得到铜水流动过程中浇铸包和中间包内铜水流动的数据，该Fluent软件是基于计算流体动力学方法的求解流动问题的通用软件；
生成模块S6，用于根据浇铸包和中间包内铜水流动的数据生成定量浇铸曲线。
在本发明实施例中，上述第一建立单元建立铜水流动过程的数学模型时，首先根据高温铜水在浇铸工具内的流动状态特征，确定流动状态特征对应的参数和约束条件。然后根据参数和约束条件，确定浇铸工具的流动状态特征对应的连续性方程、动量方程、能量守恒方程、流体控制VOF方程、湍流模型k-ε方程和UDF(User Defined Function,用户自定义函数)函数，VOF方程表示体积分数分布和自由表面运动，k-ε方程表示流体的湍流流动，UDF函数方程用于确定浇铸包与中间包的运动状态。通过连续性方程、动量方程、能量守恒方程、VOF方程、k-ε方程和UDF函数实现铜水流动过程的三维仿真。最后，将三维仿真后得到的各个模型，确定为定量浇铸过程的数学模型。其中，上述第一建立单元建立铜水流动过程的数学模型的具体过程，可通过上述实施例提供的定量浇铸曲线的获取方法中提供的数学模型的建立方法来建立，在此不再赘述。
上述浇铸工具的实际尺寸为生产现场提供的浇铸工具详细的几何尺寸，上述第二建立模块S2确定浇铸工具对应的物理模型时，可以通过Ansys中DM等比例建立浇铸工具对应的实体模型，该实体模型即为浇铸工具对应的物理模型，Ansys软件是一种有限元分析软件，DM是该Ansys软件中的建模功能单元。
在本发明实施例中，上述处理模块S3对物理模型和数学模型进行离散化处理时，首先通过Ansys中的Meshing对物理模型进行网格划分，得到空间网格区域；然后通过控制方程离散技术将数学模型在空间网格区域上进行离散，其中，Meshing为Ansys软件中的划分网格功能单元。
在本发明实施例中，处理模块S3对数学模型进行离散后，该数学模型中方程组中的各个方程得到简化，从而可以快速、准确地对简化后的各个方程进行求解，大大提高了模型的求解速度以及求解得到的数据的准确性。
在进行瞬态问题数值模拟时，需要设置初始条件。初始条件是所研究对象在过程开始时刻各个求解变量的空间分布情况，是确定偏微分方程具有唯一解的初值。边界条件是指在求解域边界上所求解的变量或其一阶导数随地点及时间变化的规律。上述确定模块S4确定浇铸过程的初始条件和边界条件时，可通过以下设置单元、控制单元和第一确定单元确定出。
上述确定模块S4包括设置单元、控制单元和第一确定单元。
设置单元，用于通过UDF函数设置浇铸工具的运动曲线；
控制单元，通过动网格技术控制浇铸工具实现运动曲线的运动过程；
第一确定单元，用于将通过动网格技术控制浇铸工具实现运动曲线的运动过程，作为浇铸过程的初始条件和边界条件。
在本发明实施例中，在确定定量浇铸过程的初始条件时，上述设置单元，需要设置初始时刻整个计算域内的压力场、气-液两相体积分数及高温铜水的温度和初始流动速度，同时由于流场和温度场的耦合计算，还需要设定初始状态下计算域内的温度场等。在Fluent软件中初始条件的设定通过计算域初始化功能选项完成。
在本发明实施例中，在确定定量浇铸过程的边界条件时，上述设置单元除了需要设置基本边界条件，如进出口边界条件、给定压力边界和wall壁面边界条件等，还要根据浇铸的目标，使用UDF函数设置浇铸工具运动曲线，用动网格技术实现浇铸工具的运动。
上述设置单元，通过预先设定的UDF(User Defined Function，用户自定义函数)函数设置浇铸工具的运动曲线时，具体可通过以下设置子单元、判断子单元、生成子单元和修正子单元进行设置。
上述设置单元包括：
设置子单元，用于根据浇铸工具的运动状态，通过Profile文件设置各个时间点中间包和浇铸包的运动速度，该Profile文件是Fluent软件中一种指定边界运动的描述方式；
判断子单元，用于判断中间包速度、倾角与浇铸铜水重量的关系是否满足预设浇铸条件，以及分析浇铸包不同位置、速度与浇铸铜水重量的关系是否满足预设浇铸条件；
生成子单元，若判断子单元判断中间包以及浇铸包均满足预设浇铸条件，用于通过Profile文件生成浇铸过程的离散时间点与速度的对应关系；
修正子单元，若判断子单元判断中间包或浇铸包不满足预设浇铸条件，用于不断增加时间点个数，修正中间包或浇铸包各个时间点的运动速度。
在本发明实施例中，上述设置子单元、判断子单元、生成子单元和修正子单元浇铸工具的运动曲线时，具体可通过上述实施例提供的定量浇铸曲线的获取方法中提供的对应方法进行处理。
上述设置单元通过UDF函数设置浇铸工具的运动曲线后，通过动网格技术控制浇铸工具实现运动曲线的运动过程；将通过动网格技术控制浇铸工具实现运动曲线的运动过程，作为浇铸过程的初始条件和边界条件。
当建立上述数学模型、上述物理模型，以及确定上述初始条件和边界条件后，上述计算模块S5对该数学模型进行求解，得到铜水流动过程中浇铸包和中间包内铜水流动的数据。
当上述计算模块S5计算得到铜水流动过程中浇铸包和中间包内铜水流动的数据后，可通过上述生成模块S6生成定量浇铸曲线。
上述生成模块S6包括：
分析单元，用于对浇铸包和中间包内的铜水流动的数据进行分析，得到铜水的体积分数和速度场。
在本发明实施例中，对浇铸包和中间包内的铜水流动的数据进行分析时，首先将该浇铸包和中间包内的铜水流动的数据转化为可视化的浇铸过程图像，从该可视化的浇铸过程图像中获取铜水在浇铸工具内的速度场和温度场，根据温度场确定出铜水的体积分数。
第二确定单元，用于根据体积分数和速度场，确定中间包和浇铸包中铜水液面流动状态以及铜阳极质量。
在本发明实施例中，上述铜水液面流动状态可以是铜水流动速度。其中，可以通过速度场和铜水的体积分数，分析浇铸工具中不同容量的铜水、不同的运动速度、不同角度下的铜水流动速度。根据该铜水流动速度与铜阳极板的质量的对应关系，确定出铜阳极板的质量。
判断单元，用于根据铜水液面流动状态及铜阳极质量判断浇铸精度是否满足预设浇铸精度；
第三确定单元，若判断单元判断浇铸精度是否满足预设浇铸精度时，根据浇铸包和中间包内铜水流动的数据确定出浇铸过程的离散时间点与速度的对应关系；
拟合单元，用于根据离散时间点与速度的对应关系拟合浇铸包和中间包的运动曲线；
第三确定单元，用于将浇铸包和中间包的运动曲线确定为浇铸工具的定量浇铸曲线；
修正单元，用于判断单元判断浇铸精度不满足预设浇铸精度时，修正UDF函数，以使重新确定边界条件。
在本发明实施例中，当修正单元需要重新确定边界条件时，具体可通过上述确定模块S4确定。
在本发明实施例中，修正UDF函数时，上述修正单元首先需要将上述拟合得到浇铸包和中间包的运动曲线用函数形式表示，编写成UDF函数。然后指定该运动曲线上特定动态区域的运动，产生随时间变化的浇铸工具运动的线速度、角速度及重心运动。最后，利用这些速度更新动态区域的位置。
当确定模块S4重新确定边界条件后，上述设置单元将得到的运动曲线函数加载到Fluent软件中，通过该重新确定的边界条件以及上述计算模块S5，以及生成模块S6中的第二确定单元、判断单元和第三确定单元重新确定浇铸精度是否满足预设浇铸精度，若不满足，上述修正单元再增加时间点的个数，然后再根据上述拟合得到的浇铸工具的运动曲线，重新确UDF函数，直至浇铸精度满足预设浇铸精度时，通过上述拟合单元和上述第三确定单元生成定量浇铸曲线。
本发明实施例提供的定量浇铸曲线的获取装置，通过建立铜水流动过程的数学模型，准确描述了高温铜水在浇铸过程中的运动状态和温度变化，通过对该数学模型进行求解并对求解得到的数值进行仿真，可以得到浇铸精度很高的定量浇铸曲线，根据该定量浇铸曲线进行铜水的定量浇铸，大大提高了生产出的铜阳极板的质量。
综上所述，本实施例提供了一种定量浇铸曲线的获取方法及装置，通过建立铜水流动过程的数学模型，准确描述了高温铜水在浇铸过程中的运动状态和温度变化，通过对该数学模型进行求解并对求解得到的数值进行仿真，可以得到浇铸精度很高的定量浇铸曲线，根据该定量浇铸曲线进行铜水的定量浇铸，大大提高了生产出的铜阳极板的质量。
本发明的实施例是为了示例和描述起见而给出的，而并不是无遗漏的或者将本发明限于所公开的形式。很多修改和变化对于本领域的普通技术人员而言是显而易见的。选择和描述实施例是为了更好说明本发明的原理和实际应用，并且使本领域的普通技术人员能够理解本发明从而设计适于特定用途的带有各种修改的各种实施例。

Claims (10)

1.一种定量浇铸曲线的获取方法，其特征在于，所述方法包括:
分析高温铜水在浇铸工具内的流动状态特征，建立铜水流动过程的数学模型；
获取浇铸工具的实际尺寸，等比例建立浇铸工具对应的物理模型；
对所述物理模型和所述数学模型进行离散化处理；
确定定量浇铸过程的初始条件和边界条件；
通过基于Fluent软件的有限体积算法、所述物理模型、所述初始条件及所述边界条件对所述数学模型进行求解，得到铜水流动过程中浇铸包和中间包内铜水流动的数据，所述Fluent软件是基于计算流体动力学方法的求解流动问题的通用软件；
根据所述浇铸包和中间包内铜水流动的数据生成定量浇铸曲线。
2.根据权利要求1所述的方法，其特征在于，所述分析高温铜水在浇铸工具内的流动状态特征，建立铜水流动过程的数学模型，包括：
根据高温铜水在浇铸工具内的流动状态特征，确定流动状态特征对应的参数和约束条件；
根据所述参数和所述约束条件，确定浇铸工具的流动状态特征对应的连续性方程、动量方程、能量守恒方程、流体控制VOF方程、湍流模型k-ε方程和用户自定义UDF函数，所述VOF方程表示体积分数分布和自由表面运动，k-ε方程表示流体的湍流流动，UDF函数方程用于确定浇铸包与中间包的运动状态；
通过所述连续性方程、所述动量方程、所述能量守恒方程、所述VOF方程、所述k-ε方程和所述UDF函数实现铜水流动过程的三维仿真；
将三维仿真后得到的各个模型，确定为定量浇铸过程的数学模型。
3.根据权利要求1所述的方法，其特征在于，所述对所述物理模型和所述数学模型进行离散化处理，包括：
通过Meshing对所述物理模型进行网格划分，得到空间网格区域，所述Meshing为Ansys软件中的划分网格功能单元；
通过控制方程离散技术将所述数学模型在所述空间网格区域上进行离散。
4.根据权利要求1所述的方法，其特征在于，所述确定初始条件和边界条件，包括：
通过UDF函数设置浇铸工具的运动曲线；
通过动网格技术控制浇铸工具实现所述运动曲线的运动过程；
将所述通过动网格技术控制浇铸工具实现所述运动曲线的运动过程，作为浇铸过程的初始条件和边界条件。
5.根据权利要求4所述的方法，其特征在于，所述通过UDF函数设置浇铸工具的运动曲线，包括：
根据所述浇铸工具的运动状态，通过Profile文件设置各个时间点中间包和浇铸包的运动速度，所述Profile文件是Fluent软件中一种指定边界运动的描述方式；
判断中间包速度、倾角与浇铸铜水重量的关系是否满足预设浇铸条件，以及分析浇铸包不同位置、速度与浇铸铜水重量的关系是否满足预设浇铸条件；若中间包以及浇铸包均满足预设浇铸条件，通过Profile文件生成浇铸过程的离散时间点与速度的对应关系；若否，不断增加时间点个数，修正中间包或浇铸包各个时间点的运动速度。
6.根据权利要求2所述的方法，其特征在于，根据所述浇铸包和中间包内铜水流动的数据生成定量浇铸曲线，包括：
对所述浇铸包和中间包内的铜水流动的数据进行分析，得到铜水的体积分数和速度场；
根据所述体积分数和所述速度场，确定中间包和浇铸包中铜水液面流动情况以及铜阳极质量；
根据所述铜水液面流动情况及所述铜阳极质量判断浇铸精度是否满足预设浇铸精度；若是，根据所述浇铸包和中间包内铜水流动的数据确定出浇铸过程的离散时间点与速度的对应关系，根据所述离散时间点与速度的对应关系拟合浇铸包和中间包的运动曲线，将所述浇铸包和中间包的运动曲线确定为浇铸工具的定量浇铸曲线，若否，则修正所述UDF函数，重新确定边界条件。
7.一种定量浇铸曲线的获取装置，其特征在于，所述装置包括:
第一建立模块，用于分析高温铜水在浇铸工具内的流动状态特征，建立铜水流动过程的数学模型；
第二建立模块，用于获取浇铸工具的实际尺寸，等比例建立浇铸工具对应的物理模型；
处理模块，用于对所述物理模型和所述数学模型进行离散化处理；
确定模块，用于确定定量浇铸过程的初始条件和边界条件；
计算模块，用于通过基于Fluent软件的有限体积算法、所述物理模型、所述初始条件及所述边界条件对所述数学模型进行求解，得到铜水流动过程中浇铸包和中间包内铜水流动的数据，所述Fluent软件是基于计算流体动力学方法的求解流动问题的通用软件；
生成模块，用于根据所述浇铸包和中间包内铜水流动的数据生成定量浇铸曲线。
8.根据权利要求7所述的装置，其特征在于，所述确定模块包括：
设置单元，用于通过用户自定义UDF函数设置浇铸工具的运动曲线；
控制单元，通过动网格技术控制浇铸工具实现所述运动曲线的运动过程；
第一确定单元，用于将所述通过动网格技术控制浇铸工具实现所述运动曲线的运动过程，作为浇铸过程的初始条件和边界条件。
9.根据权利要求8所述的装置，其特征在于，所述设置单元包括：
设置子单元，用于根据所述浇铸工具的运动状态，通过Profile文件设置各个时间点中间包和浇铸包的运动速度，所述Profile文件是Fluent软件中一种指定边界运动的描述方式；
判断子单元，用于判断中间包速度、倾角与浇铸铜水重量的关系是否满足预设浇铸条件，以及分析浇铸包不同位置、速度与浇铸铜水重量的关系是否满足预设浇铸条件；
生成子单元，若判断子单元判断中间包以及浇铸包均满足预设浇铸条件，用于通过Profile文件生成浇铸过程的离散时间点与速度的对应关系；
修正子单元，若判断子单元判断中间包或浇铸包不满足预设浇铸条件，用于不断增加时间点个数，修正中间包或浇铸包各个时间点的运动速度。
10.根据权利要求7所述的装置，其特征在于，所述生成模块包括：
分析单元，用于对所述浇铸包和中间包内的铜水流动的数据进行分析，得到铜水的体积分数和速度场；
第二确定单元，用于根据所述体积分数和所述速度场，确定中间包和浇铸包中铜水液面流动情况以及铜阳极质量；
判断单元，用于根据所述铜水液面流动情况及所述铜阳极质量判断浇铸精度是否满足预设浇铸精度；
第三确定单元，若判断单元判断浇铸精度满足预设浇铸精度，则根据所述浇铸包和中间包内铜水流动的数据确定出浇铸过程的离散时间点与速度的对应关系；
拟合单元，用于根据所述离散时间点与速度的对应关系拟合浇铸包和中间包的运动曲线；
所述第三确定单元，用于将所述浇铸包和中间包的运动曲线确定为浇铸工具的定量浇铸曲线；
修正单元，用于判断单元判断浇铸精度不满足预设浇铸精度时，修正UDF函数，以使重新确定边界条件。
CN201610424666.9A 2016-06-16 2016-06-16 一种定量浇铸曲线的获取方法及装置 Active CN105945270B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610424666.9A CN105945270B (zh) 2016-06-16 2016-06-16 一种定量浇铸曲线的获取方法及装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610424666.9A CN105945270B (zh) 2016-06-16 2016-06-16 一种定量浇铸曲线的获取方法及装置

Publications (2)

Publication Number Publication Date
CN105945270A CN105945270A (zh) 2016-09-21
CN105945270B true CN105945270B (zh) 2017-12-12

Family

ID=56906507

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610424666.9A Active CN105945270B (zh) 2016-06-16 2016-06-16 一种定量浇铸曲线的获取方法及装置

Country Status (1)

Country Link
CN (1) CN105945270B (zh)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107194033B (zh) * 2017-04-21 2020-11-24 中国人民解放军后勤工程学院 一种在开敞环境中石油泄漏蒸发量计算方法
CN107309419B (zh) * 2017-06-29 2018-10-12 中南大学 一种铜圆盘浇铸机定量浇铸过程建模方法
CN108637230B (zh) * 2018-04-23 2019-08-23 中南大学 一种铜圆盘浇铸控制方法及装置
CN111360215B (zh) * 2019-11-13 2021-06-18 浙江东阳东磁稀土有限公司 一种改善钕铁硼磁体甩片厚度分布和微观组织结构的方法
CN118682090A (zh) * 2024-08-28 2024-09-24 北京科技大学 一种连铸机钢包至中间包钢水温度的稳定控制方法和系统

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1388444A (zh) * 2001-05-25 2003-01-01 中国科学院金属研究所 一种铸件充型过程模拟方法
JP2004506515A (ja) * 2000-07-01 2004-03-04 エイ・イー・エム・ピー・コーポレーション 鋳造／成形課程用の熱流れシミュレーション
CN101727530A (zh) * 2010-01-07 2010-06-09 冶金自动化研究设计院 一种实现连铸过程仿真系统及其方法
CN101833605A (zh) * 2010-04-29 2010-09-15 浙江工业大学 基于流体体积模型模具微细流道磨料流精密加工控制方法
CN105302985A (zh) * 2015-11-12 2016-02-03 哈尔滨工业大学 一种基于fluent软件的合金微铸造成形过程的仿真方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2004506515A (ja) * 2000-07-01 2004-03-04 エイ・イー・エム・ピー・コーポレーション 鋳造／成形課程用の熱流れシミュレーション
CN1388444A (zh) * 2001-05-25 2003-01-01 中国科学院金属研究所 一种铸件充型过程模拟方法
CN101727530A (zh) * 2010-01-07 2010-06-09 冶金自动化研究设计院 一种实现连铸过程仿真系统及其方法
CN101833605A (zh) * 2010-04-29 2010-09-15 浙江工业大学 基于流体体积模型模具微细流道磨料流精密加工控制方法
CN105302985A (zh) * 2015-11-12 2016-02-03 哈尔滨工业大学 一种基于fluent软件的合金微铸造成形过程的仿真方法

Also Published As

Publication number Publication date
CN105945270A (zh) 2016-09-21

Similar Documents

Publication Publication Date Title
CN105945270B (zh) 一种定量浇铸曲线的获取方法及装置
Demirdžić et al. Numerical method for coupled fluid flow, heat transfer and stress analysis using unstructured moving meshes with cells of arbitrary topology
Boomsma et al. Simulations of flow through open cell metal foams using an idealized periodic cell structure
Macklin et al. An improved geometry-aware curvature discretization for level set methods: application to tumor growth
JP7408276B2 (ja) 中子造型装置および中子造型装置の制御方法
CN105302985A (zh) 一种基于fluent软件的合金微铸造成形过程的仿真方法
US5572434A (en) Method for simulating mold filling of semi-solid material
Baltussen et al. Numerical simulation of the fountain flow instability in injection molding
CN113947003A (zh) 一种面向热流耦合场景的粒子型无网格仿真系统
Ng et al. Effect of solder bump shapes on underfill flow in flip-chip encapsulation using analytical, numerical and PIV experimental approaches
Pang et al. Numerical simulation two phase flows of casting filling process using SOLA particle level set method
CN108197072B (zh) 一种基于加权守恒变量阶跃的高精度间断Galerkin人工粘性激波捕捉方法
Grucelski et al. Lattice Boltzmann simulation of fluid flow in porous media of temperature-affected geometry
Binding et al. Experimental and numerical simulation of dough kneading in filled geometries
McBride et al. Complex free surface flows in centrifugal casting: Computational modelling and validation experiments
Yan et al. Numerical simulation of AZ91D alloy automobile plug in pressure die casting process
CN111931392B (zh) 等离子加热中间包底吹氩参数优化方法和装置
Young Thermal behaviors of the resin and mold in the process of resin transfer molding
JP4006316B2 (ja) 樹脂流動解析方法及びその装置
CN111695307A (zh) 显式考虑动态摩阻的水锤有限体积模拟方法
Friedrichs et al. A novel hybrid numerical technique to model 3-D fountain flow in injection molding processes
CN104368800A (zh) 模拟铝铸件中的氧化物的方法
Vakhrushev et al. Experimental verification of a 3-phase continuous casting simulation using a water model
CN110013331B (zh) 一种蜡制颌堤模拟牙齿移动过程中动态阻力模型建立方法
CN113505544A (zh) 基于有限体积法自行车运动虚拟数值风洞系统

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant