1、第三章 非稳态导热,unsteady state heat conduction,3-1 非稳态导热过程的特点,一、定义导热体的温度随时间而变化的导热过程称为非稳态导热工程中的许多过程都是非稳态导热：冶金、热处理和热加工；工件被加热和冷却；锅炉、内燃机、燃气轮机等装置启动、停机、变工况；自然环境温度，供暖和停暖过程中墙内与室内空气温度。原因：边界条件，内热源发生了变化,二、非稳态导热过程及分类非稳态导热过程：周期性和非周期性（瞬态导热）周期性非稳态导热：物体温度按照一定的周期发生变化。例如：墙体温度的昼夜变化非周期性非稳态导热（瞬态导热）：物体温度随时间的推 移趋近于恒定的值 。例如：热铁块投
2、入凉水中,非稳态导热过程温度变化曲线,三、非稳态态导热过程的特点,物体表面层最先受热，然后由表及里地逐渐传播到物体内部。理论上经过 无限长时间后，物体内各处温度趋于一致，此时达到新的热平衡状态。,钢棒和木棒的次端哪一个温度变化剧烈？,如果两者的导热系数相同，但是，Cp不同，又有何影响？,思考：,影响平壁内温度变化快慢的因素是什么？,反映了物体的导热能力和储热能力之比，它可以用来衡量物体在加热或冷却时内部温度变化传播速度。越大，说明在相同的温度梯度下可以传递更多的热量； 越小，单位体积的物体温度升高1所需要的热量越小，可以剩下更多的热量继续向物体内部传递，使得物体内部各点温度趋于一致的能力提高。
3、,四、导温系数（热扩散系数）,五、毕渥准则对温度分布的影响,六、物体加热、冷却时其内部温度场变化的三个阶段,不规则状况阶段过程开始时，初始温度的影响,正规阶段物体内各点的温度变化具有一定规律,新的稳态阶段,七、非稳态导热的研究目的,温度分布和热流量分布随时间和空间的变化规律非稳态导热的导热微分方程式：求解方法：分离变量法、积分变换法 近似分析法：集总参数法、积分法、瑞利-里兹法 数值解法：有限差分法、有限元法,3-2 集总参数法,基本思想：对任意形状的物体，忽略物体内部的导热热阻，认为物体温度均匀一致的分析方法。,温度t只是时间的函数，与坐标无关 零维问题。,使用条件：,体积V 表面积F 物性
4、、cp、 温度,导热体,流体,Q,已知：任意形状的物体,求：,Q,体积V 表面积F 物性、cp、 温度,导热体,流体,Q,解：,任意时刻物体表面与周围流体 对流换热的热流量为：,此时物体单位时间内能的变化为：,由能量平衡，物体单位时间与流体的换热量 应等于物体内能的变化：,负号是因为：,(a),(c),(b),引入过余温度 ：,(d),(e),将式(d)分离变量，并应用初始条件，对时间从0到积分，有:,(3-1),(3-1),(3-2),毕渥准则 ， 具有长度的量纲,傅立叶准则，为无因次时间,无因次的过余温度,物体内温度随时间呈指数衰减,令,时间常数，记为,热电偶的时间常数是说明热电偶对流体温
5、度变化响应快慢的指标。,时间常数,过余温度降至初始过余温度 的36.8%时所需的时间。 时间常数越大过余温度变化越慢。,由式(3-2)可求出任意瞬时的热流量：,(3-3),从 0 到 时刻之间所交换的总热量为：,(3-4),使用条件：,对于无限大平板、柱体和球体，若满足下面条件：,(3-5),采用集总参数法计算的误差不超过5%（与精确解比较而知）,对于厚度为2 的大平板： ，对于半径为r的长圆柱体： ，对于半径为r的球体： ，,毕渥准则采用一般的表达形式 ：,应用集总参数法计算的条件为 ：,例题3-1 一直径为50mm的钢球，初始温度为450，突然 被置于温度为30的空气中。设钢球表面与周围空
6、气的对流 换热系数为24W/(m2K)，试计算钢球冷却到300时所需要 的时间。已知钢球的,解： 首先检验是否可以采用集总参数法。,可以采用集总参数法计算。,由式(3-1),例题3-2 热电偶是用来测量温度的一种感受件。如果用其测量 流体的温度，要求此感受件在温度变化较大的范围内反应时滞 （即从测温开始至到达真实温度的时间）不大于1.5s，问在暂 不考虑沿热电偶丝的导热影响时热电偶的热接点应该多大。已知 热电偶的 热接点表面的对流换热系数 ，而且不随直径变化。,解：任何感受件从测温开始到指示出真实温度所需的时间，理论上 为无限长。在工程实际中可以按照允许误差为4来设计热电偶的 热接点。假定热接
7、点的形状为球形。,设初始的流体温度为t0，变化后的流体温度为tf。流体温度 变化1.5s后，热接点的温度将迅速变化t。,t0,t,tf,t0,t0,t0,tf,按题意，测量误差为|t-tf|=tf4%，即有,按题意，因温度变化范围较大，故可认为 ，则,故热接点的直径为：,例题3-3 一直径为5cm、长为30cm的钢圆柱体，初始温度为 30将其放入炉温为1200的加热炉中加热，升温到800方可 取出，设钢圆柱体与烟气间的表面对流换热系数（其中包括烟气 与表面之间的辐射换热）为 ，钢的物理性质参数为。问需要多长 时间才能达到要求。,解 首先检查是否可以采用集总参数法。,可采用集总参数法。,由式(3
8、-1)得：,S,无限大平板的坐标系统,3-3 无限大平板的加热或冷却,问题： 无限大平板，、cp 为 常数，厚度为2，初始温度 均匀为t0，将其突然置于温 度恒为t f的流体中，t ft0. 流体与板侧面之间的对流换 热系数为。求:板内温度t=f(,x)及物体所吸收的热量Q()。 解： 一维问题。,t0,无限大平板的坐标系统,数学模型： 导热微分方程：初始条件：边界条件：（第二类）（第三类）,，,t0,引入过余温度 ，则数学模型变为：,(a),(b),(c),(d),用分离变量法进行求解。令,(e),则式(a) 变为,(f),Why minus,即,(g),两个常微分方程的通解为：,(h),导
9、热微分方程的通解为：,(i),由边界条件式(c)得：,式(i)变为：,(j),代入到边界条件式(d)中得：,(k),关于自变量的超越方程(k)有无穷多个根，记为n。,o,对应于每一个 ，式(j)中都有一个待定常数 与之对应，这样 可以得到无穷多个特解，这些特解的和即为微分方程的通解：,(l),根据初始条件式(b)：,(l),三角函数的正交性 ：,式(l)两边同乘 并积分：,即,（m）,令 无因次过余温度 的函数傅立叶准则，无因次时 间； 无因次坐标,(3-6),(3-7),板内温度场，为时间、位置、板的物性及换热特性的函数,能量守恒原理：到时刻平板所吸收的热量物体焓值的变化 每平方米平板吸收的
10、热量为：,(3-8),式中，Q0为平板从初始温度t0加热到流体温度tf时每平方米平板所 吸收的热量。,当 后，只采用级数的第一项进行简化计算，误差小于1% ：,（Why： 是减函数 , ）,(3-9),(3-10),由式(3-9)，当 后，平板中任意一点的过余温度与平板 中心的过余温度之比为：,比值与时间无关，初始条件的影响已经消失。,正规阶段,(3-1),(3-2),工程上常利用正常情况阶段的特点来测定材料的 热物性参数，如导温系数、比热容或导热系数。测同一点不同时刻的 值，利用关系式求物性值。,圆柱体和球体的计算：作为一维非稳态问题也可以用 分离变量法求解。参考有关文献。,诺模图： 按分析
11、解的级数第一项而绘制的。使用条件：,由式：,板内某点、某时刻的温度,平板单位面积的换热量计算：,(3-10),由式：,画出曲线图如下：,查出 的值，计算出 值后，可以求得 的值。,诺模图：简捷方便，计算的准确度受到有限的图线的影响。通过编程直接应用分析解来计算 或 直接采用数值计算。,3-4 半无限大物体的非稳态导热,“半无限大物体”：研究热流方向仅沿着x流动的导热问题。一维非稳态导热问题 工程应用：大部件的表面淬火。四周被绝热的长杆，初始温度均匀一致，一端被冷 却或加热时，也是一种等效的半无限大物体。,已知：半无限大物体，初始温度为t0，表面突然与温度为tf的流体 相接触，表面温度从t0升高
12、到tf，并保持不变。假物性参数为常数 求：此非稳态导热一维问题。,解：,导热微分方程式：,初始条件：边界条件：,(b),(c),(a),由拉普拉斯变换法，此问题的解为：,(3-12),(3-13),是由拉氏变换引入的虚变量。 积分是其上限的函数。,温度分布：,(3-14),半无限大物体中的温度分布,任意点的热流量：,对式(3-14)求偏导数得：,通过表面的热流量：,误差函数的数值可参看有关数学书籍及附录中附表。,3-5 有限大物体的非稳态导热,有限大物体： 两方向或三方向尺寸相差不大的物体。它们的非稳态导热是二维或三维非稳态导热问题,求解方法： 对于这些形状规则的有限大物体的非稳态导热问题，可以利用一维非稳态导热问题分析解的组合求解。,无限长方柱体的温度场：等于厚度分别为21及22的两块无限大平板的温度场的乘积。无限长方柱体的初始条件及边界条件分别为两块无限大平板的 对应边界上的初始条件和边界条件。,(3-17),无限长方柱体的无量纲过余温度厚度为21无限大平板的无量纲过余温度厚度为22无限大平板的无量纲过余温度,图3-18 无限长方柱体,对于短圆柱体 和立方体的温度场问题，方法类似。,本章习题：1、4、8、9,
