在线工博会
返回

铸造合金凝固模拟潜热处理综合模型研究
徐宏 侯华 赵宇宏 陈铮
为节省流量，手机版未显示文章中的图片，请点击此处浏览网页版
摘要：铸造合金潜热释放模型对铸件凝固过程温度场模拟有很大影响。本文对合金潜热处理温度回升法进行修正，建立了一种适用于任意合金的凝固综合潜热处理方法。并经实际应用进行了考核。
关键词：铸造合金；潜热；温度场；模拟
凝固过程伴随着潜热释放是铸件成形过程的一大特点，铸件在向铸型所传递的热量中，金属过热的热量仅占20%左右，凝固潜热约占80%[1]。潜热因素对铸件温度场模拟影响很大，必须加以考虑。通常对潜热处理有如下几种方法：温度回升法[2]、等效比热法[3]及热焓法[4]。
温度回升法适用于纯金属或共晶成分合金的潜热释放过程， 等效比热法和热焓法适用于处理有一定结晶温度范围合金。采用热焓法需要知道材料热焓与温度的关系，由于这方面资料较少，因此这方法在国内基本没使用。chen 和Tsai[5]研究指出不同合金的热释放模型是不一致的，合金潜热释放模型对合金凝固过程温度场计算分布会造成重大影响。处理结晶潜热的好坏对铸件凝固数值计算精度起着非常关键的作用，研究适用于任意合金的凝固潜热处理方法是目前凝固模拟研究重点之一。
1 铸件凝固传热数学模型
铸件凝固过程是一个不稳定导热过程，温度、时间和空间的关系可用傅立叶（Fourier）导热偏微分方程描述：

（图片）

式中T为温度；C为比热；ρ为液相和固相的平均密度；Q为潜热释放量；λ为导热系数；u，v，w为枝晶间液体流速，fL为液相率。
式（1）中左边第二项是枝晶间液体流动的热对流，由于潜热释放、两相区温度梯度及液相率比较小，可以忽略不计，公式（1）可描述为：

（图片）(2)

在铸件凝固过程中，从液相到固相的相变过程中释放的结晶潜热，其单位时间单位体积内的释放量就是方程中的内热源Q，如何处理结晶潜热是铸件凝固模拟成功的关键之一。
2 潜热处理综合模型
铸件凝固冷却过程实质上是铸件内部显热（温度下降）和潜热（相变结晶）不断向外散失的过程。
2.1 合金固相率变化规律
实际上，合金固相率的变化规律在凝固理论研究中一直受到高度重视。它与凝固潜热的模型化处理、缩孔（松）判据等机理性研究密切相关。已提出的固项率fs表达式有[6-8]平衡凝固模型、Scheurer - Scheil模型、Brody - Flemings模型、Clyne - Kurz模型及介万齐 -- 周尧和模型等。这些理论计算式以溶质质量守恒为基础，对溶质传输规律作了不同的假设，对于二元固溶体型合金的凝固模拟及其机理性研究起着积极的指导作用。但对于多元复杂相变合金，在理论上，目前还难以研究。由于复杂合金的固相率变化规律研究目前仍存在困难，在当前凝固模拟中， fs主要按数学简化模型处理，如式（3）（4）。

（图片）（3） （4）

2.2 修正温度回升法
本文对温度回升法进行修正，从而使之能适应于有结晶温度范围及发生共晶发应的合金。
对于有结晶温度范围的合金，根据上一时刻单元温度和当前时刻温度回升前后值的不同情况，将潜热释放分为六种模型[9]，温度回升方法如下：
当某一体积的单元温度进入固液两相区（固液两相区由合金液相线温度和固相线温度确定），在某一时间步长dt，单元温度低于液相线温度时释放的能量通过下式来确定。

（图片）（5）

式中，V是单元的体积，L为合金的结晶潜热， 为dt时间内单元固相率增加量 。
假设凝固潜热释放是线性的，即凝固过程释放的潜热与生成固相率成正比，释放的潜热使单元本身温度回升 ，相应的内能增加为

（图片）（6）

上式中Cp为比热。
根据能量守恒原理：

（图片）（7）

由式（5）（6）（7）得：

（图片）（8）

假设结晶潜热在液固两相区内成线性释放，则：

（图片）（9）

根据公式（8）和（9），可以计算六种模型的不同情况：
模型（1）：上一时刻的单元温度Tn在液相线Tl以上，当前时刻的单元温度回升前后都在固液相线之间：

（图片）（10）

（图片） （11）

模型（2）：上一时刻的单元温度Tn在液相线Tl以上，当前时刻的单元温度回升前在固相线Ts以下，回升后的温度也在固相线Ts以下：

（图片）（12）

（图片）（13）

模型（3）：上一时刻的单元温度Tn在液相线Tl以上，当前时刻的单元温度回升前在固相线Ts以下，回升后的温度在固相线Ts以上：

（图片）（14）

（图片）（15）

模型（4）：上一时刻的单元温度Tn在液相线Tl与固相线Ts之间，当前时刻的单元温度回升前后也在液相线l与固相线之间：

（图片）（16）

（图片）（17）

模型（5）：上一时刻的单元温度Tn在液相线Tl与固相线Ts之间，当前时刻的单元温度回升前在固相线以下，回升后在固相线以下：

（图片）（18）

（图片） （19）

模型（6）：上一时刻的单元温度Tn在液相线Tl与固相线Ts之间，当前时刻的单元温度回升前在固相线以下，回升后也在固相线以上：

（图片）（20）

（图片） （21）

对于在恒温下发生的共晶反应问题，全部释放潜热量QT为：

（图片）（22）

在某一时间步长dt，温度低于液相线温度时释放的能量为：

（图片）（23）

当Q>QT，此时实际温度为：

（图片）(24)

式中（ρc）s是固相线温度下的合金热容量。
当Q（图片）(25)

3 实际模拟考核
实例一：材质为ZG35的铸件三维实体如图1(a)，实验解剖缺陷结果如图1(b)所示，模拟预测结果如图1(c)所示。

（图片） （图片） （图片）

（a）三维实体 （b）解剖结果 （c）缩孔模拟预测结果
图1 实例一应用结果

实例二： 图2为一兵器壳体铸件，材质为ZL104，铸型为金属型，砂芯材质为水玻璃砂，浇注温度为700°，铸件产生缩孔缩松缺陷（图4所示）。模拟温度及缩松预测结果如图3和图5所示。实际结果与模拟结果吻合。

（图片） （图片）

图2 铸件三维图 图3 温度模拟结果

（图片） （图片）

图4 缩松缺陷位置 图5 缩松模拟预测结果

通过实际应用对比考核可以看出，利用该潜热处理综合模型模拟结果与实际检测缺陷结果一致。这验证了本文提出的潜热处理综合模型的正确性。
4．结论
对温度回升法修正建立的铸造合金潜热处理综合模型能处理任意结晶温度范围合金的潜热释放计算，通过实际验证，证明该综合数学模型能普遍用于铸件凝固模拟计算中。
参 考 文 献
[1] 程军. 计算机在铸造中的应用. 北京： 机械工业出版社， 1993
[3] 陈海清等. 铸件凝固过程数值模拟. 重庆大学出版社，1991
[8] 介万齐，周尧和. 二元固溶体型合金凝固固相分数变化规律. 西安：西北工大科技资料， 1986.
[9] 贾良荣. 压力铸造过程数值模拟技术的研究. 清华大学博士论文，2000 2/20/2005

电脑版 客户端 关于我们
佳工机电网 - 机电行业首选网站