凝固过程的传热

20
§1-1 凝固过程的传热特点 一热：有热源的非稳态传热过程，是第一重要的。 二迁：固相、液相间界面和金属铸型间界面，而这二个界 面随着凝固进程而发生动态迁移，并使传热现象变得更加
复杂。
三传：液态金属的凝固过程是一个同时包含动量传输、质 量传输和热量传输的三传（导热、对流和辐射传热）耦合 的三维传热物理过程。
21
§1-1 凝固过程的传热特点
在实际生产中，铸件形状和材料种类的多样性以及材料热 物性值随温度非线性变化的特点，也都使凝固的传热过程 变得十分复杂。
22
§1-1 凝固过程的传热特点
传热有热源的非稳态传热过程，导热微分方程为：
23
§1-1 凝固过程的传热特点
方程(1-1)是均质、各向同性体的传导微分方程，它表示热 传导过程的能量守恒原理。事实上，方程左侧括弧内各项， 是热流密度（单位时间、单位面积上通过的热量）在x，y和 z坐标上的分量，如 , 因此，方程前三项即是热
附近采用密网格，而远离界面处采用疏网格。
53
§1-4凝固过程的电子计算机数值模拟
二、建立数学模型
1．微分方程转变为差分方程
在凝金属才释放，因此，对于在这个温度范围以外的金属和在铸 型中，潜热值为0，方程(1-36)变为：
54
§1-4凝固过程的电子计算机数值模拟
(2）将原问题的界面热阻视为常数，即界面传热系数hi是常数；
(3)金属平面晶前沿在固定的凝固点Ti下凝固； (4)忽略液体金属的过热度和对流；
(5)铸件和铸型的热物性值视为常数。
42
§1-3金属型铸造的凝固传热
43
§1-3金属型铸造的凝固传热
44
§1-3金属型铸造的凝固传热
45
§1-3金属型铸造的凝固传热
铸件和铸型上虚拟加厚的凝固层和铸型厚度，即图中S0和- Eo上
的热阻，同时令这两个热阻上的温度降恰好等于界面上的温度降 。这样，就把一个具有界面热阻的复杂的传热问题，转变成了在
界面上理想接触因而具有共同的界面温度Ti的纯导热问题。
41
§1-3金属型铸造的凝固传热
为简化求解过程，作如下假定： (1)问题局限于一维热传导，金属型为半无限大；
图中采用了均匀网络，即令Δ X = Δ y，上下边界单元的高度和左右边 界单元的宽度各取内部单元相应尺寸的一半。
图中各单元中的点表示各该单元的温度参考点，即用这些点的温度来表
示各点所在单元的温度。 在实际计算中，经常采用疏密不均匀网格，在需要仔细计算的部位和温
度梯度较大的部位，网格划分的比其他部位密一些，一般是在铸件铸型界面
先把这一方程转变为差分方程，然后再研究凝固潜热问题。用差分来代
替微分，即可将这一微分方程转变为差分方程。 微分和差分的关系是：
55
§1-4凝固过程的电子计算机数值模拟
T关于x的向前、向后和中心差分，由图1-8可知：差分就是用函数曲线 上一个或两个单元间的割线来代替曲线上的切线，因此差分是一个近似表达 式。由图还可以看出，中心差分的准确度高于其它两种形式的差分。
32
§1-2非金属型铸造的凝固传热
33
§1-2非金属型铸造的凝固传热
于是，求温度场的问题简化成了求一维偏微分方程的问题， 求解如下：
34
§1-2非金属型铸造的凝固传热
35
§1-2非金属型铸造的凝固传热
36
§1-2非金属型铸造的凝固传热
37
§1-2非金属型铸造的凝固传热
由式(1—14)可知，金属和铸型的热物性结合起来决定凝固速度：在金属
学模型、编制程序和计算。下面着重介绍前两个步骤。
51
§1-4凝固过程的电子计算机数值模拟
一、单元剖分 有限差分法计算中，通常是将一般的铸件和铸型系统剖分为许多六面体 单元，对可以用平面二维方法处理的铸件，则是将铸件和铸型系统的某一断 面划分为许多四边形单元。
52
§1-4凝固过程的电子计算机数值模拟
46
§1-3金属型铸造的凝固传热
47
§1-3金属型铸造的凝固传热
48
§1-3金属型铸造的凝固传热
49
§1-4凝固过程的电子计算机数值模拟
在上两节中看到，即使是一维传热的简单铸件，用解析法计算温度场或 凝固时间，就已经显得相当繁杂，而实际铸件绝大多数都是具有二维或三维 传热的形体，要用解析法求解就遇到很大困难，于是产生了数值计算的方法。 常用的数值计算方法有三种：有限差分法、有限元法和边界元法，本节 只介绍较易掌握的有限差分法。这种方法将计算对象--铸件和铸型系统剖分 为许许多多有限小尺寸的单元体，假定每个单元体之向的温度梯度为常数， 在每个单元体上建立代数方程来代替以无限小单元体为基础建立的微分方程， 形成以与单元体数相等的方程组成的代数方程组，最后用计算机求解这一通 常是十分庞大的方程组。
26
§1-1 凝固过程的传热特点
27
§1-1 凝固过程的传热特点
在凝固问题的研究中，计算动态系统各点温度时间的变化即温度 场和计算凝固速度是非常重要的，因为它们直接影响金属的结晶组织、 铸件的缩孔，缩松，应力状态及许多重要的使用性能，人们已进行了 很多计算温度场和凝固速度的研究，解决的途径有解析法和非解析法。 其中解析法常受这样的限制：即使是一维传热的简单铸件，只要涉及 凝固过程，就必须作一系列假定才能求解，而且计算过程也过于繁杂， 至于形状复杂的俦件，根本无法计算。
金属凝固原理
第一章凝固过程的传热
1
前言
2
前言
3
前言
4
前言
5
前言
6
前言
7
前言
8
前言
9
前言
10
前言
11
前言
12
前言
13
前言
14
前言
15
前言
16
前言
17
§1-1 凝固过程的传热特点
凝固过程首先是从液体金属传出热量开始的，高温的液体金属浇入 温度较低的铸型时，金属所含的热量通过液体金属、巳凝固的固体金属、 金属-铸型的界面和铸型的热阻而传出；从另一个角度考察，在凝固过 程中，金属和铸型系统内发生热的传导、对流和辐射。图1-1是纯金属 浇入铸型后发生的传热模型示意图。
流密度在x、y和z轴单位长度上的增量，综合这三项就是单 位体积上的热流密度的增量，而方程的右端项，则是单位体 积的物体在单位时间内增加的内能。
24
§1-1 凝固过程的传热特点
其次，在金属凝固时存在着两个界面。即固相、液相间界面和金属
铸型间界面，而在这些界面上，通常发生极为复杂的传热现象。如一
个从宏观上看是一维传热的单相凝固的金属，当其固液界面是凹凸不 平的或生长为枝晶状时，在这个凝固前沿上，热总是沿垂直于这些界 面的不同方位从液相传入固相，因而发生微观的三维传热现象。在这 个微观区域，除了与界面垂直的热传导外，同时发生液相的对流，使 这里的传热过程十分复杂。
25
§1-1 凝固过程的传热特点
在金属与铸型的界面，由于它们的接触通常不是完全的，所以它
们之间存在接触热阻或称界面热阻，在金属凝固过程中，由于金属的
收缩和铸型膨胀，它们的接触情况也不断地在变化，在一定的条件下， 它们之间会形成一个间隙（也称气隙），所以，在这里的传热也不只 是一种简单的传导，而同时存在微观的对流和辐射传热，如图1-2所示。
30
§1-1 凝固过程的传热特点
在金属型铸造、压铸或连续铸造中，通常界面热阻Ri值远大金属和 铸型热阻，因此采用准确的hi值，是取得准确结果的关键。严格地说， hi值是随凝固时间而变化的，但是其值只是在浇注初期有较大幅度的 变化，此后较为平稳，所以常以常数处理。
31
§1-2非金属型铸造的凝固传热
18
§1-1 凝固过程的传热特点
19
§1-1 凝固过程的传热特点 凝固过程的传热有如下一些特点： 简单地说：一热、二迁、三传。 首先它是一个有热源的传热过程。金属凝固时释放的
潜热，可以看成是一个热源释放的热，但是金属的凝固潜
热，不是在金属全域上同时释放，而只是在不断推进中的 凝固前沿上释放。即热源位置在不断地移动；另外，释放 的潜热量也随着凝固进程而非线性地变化。
热流的限制环节通常不在铸型，而在铸件与铸型之间的界面， 当铸件凝固收缩和铸型受热膨胀而在铸件-铸型间形成气隙时， 界面热阻的作用将变得更为突出。
39
§1-3金属型铸造的凝固传热
40
§1-3金属型铸造的凝固传热
为解析具有这种界面温度降的传热问题，这里引进虚拟凝固
层厚度和虚拟铸型厚度的概念。即将图1-6 a分解为图1-6b和图 1-6c使后两者的组合等效于前者。这种方法将界面热阻转化成了
28
§1-1 凝固过程的传热特点
非解析法有图解法、电模拟法和数值模拟法等，自从电 子计算机问世以来，数值模拟法得到了迅速的发展。 主导方程(1-1)是均质、各向同性体的传导微分方程，
它表示热传导过程的能量守恒原理。
29
§1-1 凝固过程的传热特点
在凝固过程中，如果不计液体金属的热阻，金属的凝固速 度主要受如下三种热阻的控制：
50
§1-4凝固过程的电子计算机数值模拟
在凝固过程中，除传热现象以外，还伴随许多物理现象，如凝固潜热的
释放，液体金属的对流，金属的收缩等，因此，在计算中必须同时考虑这些
因素，采用的方法是，根据这些物理现象发生的条件，不断模拟这些现象而 变换计算过程。因此，凝固问题的数值方法，通常称为数值模拟法。
用有限差分法进行数值模拟，按如下四个步骤进行：单元剖分、建立数
非金属型的特点是，与浇注于其中的金属相比具有非常小的导热系 数，因此，金属的凝固速度主要决定于铸型的传热性能，而很少受金 属传热性质的影响。 由于铸型的导热能力差，在金属凝固的全过程中，铸型外表面的 温度变化不大，所以可以将铸型看作是半无限厚的。下面分析一个无 限大平板在这种铸型中凝固的情况。浇注的金属假定为纯金属，浇注 温度取为其熔点，即金属无过热度，这时，金属-铸型系统的温度分布 如图1-3所示。