计算机在材料科学中的应用-第六章.pptVIP

第六章 材料温度场的计算 主要参考书 1）传热学，曹玉璋，北京航空航天大学出版社。 2）计算机在材料科学中的应用，曾令可，叶卫平，武汉理工大学出版社，2004。 主要内容 热传递的基本方式及表述方法 傅立叶定律及热传导微分方程 边界条件的引入 有限差分法求解微分方程 平板样品的热传导问题 材料温度场的数值计算及其应用 温度场与材料的结晶，热处理，相变等工艺密切相关。 * * 一、材料热传导基本原理 1.1 热传递的基本方式： ① 热传导：由于温度不同，在导体内存在温差和温度梯度引起自由电子移动，导致热量传递。一般来说，好的导电体也是好的导热体。 ② 热对流：由于运动的流动质点发生相对位移，而引起的热能转移现象。 ③ 热辐射：物体受热后，内部原子振动而出现的一种电磁波能量传递，一切物体只要其温度高于绝对零度，就会从表面放出辐射能。 1.2 热传导的热量传递： 热流量：单位时间内通过某一表面所传递的热量称为热流量，用f来表示。 热流密度：单位时间，单位面积所传递的热量称为热流密度，用q来表示。 根据实验结果，热流量可以表示为： l为比例系数，热传导率，DT 为温差，A为传热面积，d为温差距离。 1.3 热对流的表示： A为面积，a为对流换热系数， T－传热流体温度，Tf为流体温度，l为导热系数。 1.4 热辐射 T1，T2分别为两个辐射体的温度，d为玻尔兹曼常数，Rr为辐射换热热阻。 二、热传递的基本公式 温度场：在一个热传递系统中，温度通常可以是空间和时间的函数，在直角坐标系中，可用函数的形式表示： t为时间，x，y，z为空间坐标，在某一瞬间，空间各点温度的集合，称为温度场。如果温度场不随时间变化而变化，称为稳态温度场。 又可按照空间维数分为：一维，二维和三维温度场。 2.1 傅立叶公式 在dt时间内通过微元等温面dA的热量dQ，正比于温度梯度 ，其方向与温度梯度方向相反。 可用如下公式表示： 热流密度 l可以用 ，l0为0°C时导热系数， b为实验常数。 温度梯度 时间 面积 2.2 导热微分方程 X Y Z Qz+ dz Qz 根据傅立叶定律： 温度梯度 导热系数 时间 根据能量守恒定律：导入微元体的热量与导出微元体的热量之差等于微元体内能的增量，所以有： r为密度，C比热容。 根据傅立叶定律： 则： 所以： a－等温系数或热扩散率。 稳态扩散时： 如果内部有热源时： qi－单位时间，单位体积导热物体生成热。 稳态导热方程为： 3. 导热微分方程的单值条件：即导热问题的初始条件和边界条件。 初始条件：物质开始导热的瞬间温度分布。 即：t＝0时的T＝T（x，y，z），最简单的情况是： t ＝0，T＝常数 边界条件可以分为： 1）第一类边界条件：直接给出物体表面温度Tw随时间的变化情况，即： 2）第二类边界条件：给定物体边界上的热流密度qw的分布及随时间的变化规律，即： 如果qw＝0，表示没有热量的出入，是绝热条件。 3）第三类边界条件是对流边界条件。 导热问题的数值求解 稳态温度场中，有一个厚度为d的平板，板的两侧温度分别为Tw1，Tw2，求其温度场的分布即热流。 根据已知条件：稳态温度场为： 平板样品只沿x方向传热，y，z方向无导热。因此可将导热微分方程写成： 边界条件： 求解上述的微分方程可得： 涓埃那嘎边界条件代入求解系数：可得到： 玻尔兹曼常数 导热微分方程得有限差分解法 一阶微商： 向前差商 向后差商 中心差商 热流密度 二阶微商： 区域离散化： 二维问题用离散点Pij得集合代替连续区域R。对于非稳态问题，取第三个坐标t，另其方向从板面向上，并将二维网格在Dt间隔内向外重复而形成一个三维网格系，网格间距为Dx，Dy，Dt。 x y Pij i j 在差分网格中，一维问题用线段，二维问题用矩形，三维问题用长方体，温度变化大的地方要划分细一些。 节点有内部节点和外部节点。每个节点都集中了它周围的热容，节点的温度代表着它周围的平均温度，节点温度T（xi，yj，zk，tm）代表连续区域的温度分布。 ，简写成 在tm时刻，（xi，yj，zk）处的温度值。 内部节点差分方程 厚度为l的无限大板材，温度分布T＝f（x，t），假定初始条件 为： ，l，C，r均为常数。 x t i m 1 i＋1 1 m＋1 根据导热差分方程：建立差分方程的x沿空间的步长为Dx，Dt为时间步长。 在i节点m时刻有： 用微商代替微分： 有： 整理可得： 取： 整理得： 本差分方程的收敛与否与温度项的系数是否大于0有关，为了保证上述递推过程收敛必须保证温度项系数大于零。即： 2F01 思考题：有一厚度为200mm的平板工件，一段的温度恒定为300C，另一