传热学第三章

1、2022-1-61主讲：陈丽湘主讲：陈丽湘华中科技大学文华学院机电学部华中科技大学文华学院机电学部2022-1-62 3-1 非稳态导热过程非稳态导热过程3-2 集总参数法集总参数法3-3 一维非稳态导热的分析解一维非稳态导热的分析解 2022-1-63定义：定义：导热系统内温度场随时间变化的导热过导热系统内温度场随时间变化的导热过程为非稳态导热。程为非稳态导热。特点：特点：温度随时间变化，热流也随时间变化。温度随时间变化，热流也随时间变化。自然界和工程上许多导热过程为非稳态自然界和工程上许多导热过程为非稳态: t = f( )例如：冶金、热处理与热加工中工件被加热或例如：冶金、热处理与热加工

2、中工件被加热或冷却；锅炉、内燃机等装置起动、停机、变工冷却；锅炉、内燃机等装置起动、停机、变工况；自然环境温度；供暖或停暖过程中墙内与况；自然环境温度；供暖或停暖过程中墙内与室内空气温度室内空气温度2022-1-64非稳态导热：非稳态导热：周期性和非周期性周期性和非周期性（瞬态导热）（瞬态导热）周期性非稳态导热：周期性非稳态导热：在周期性变化边界条件下发在周期性变化边界条件下发生的导热过程，物体温度按一定的周期发生变化。生的导热过程，物体温度按一定的周期发生变化。非周期性非稳态导热：非周期性非稳态导热：在瞬间变化的边界条件下在瞬间变化的边界条件下发生的导热过程，物体的温度随时间不断地升高发生的

3、导热过程，物体的温度随时间不断地升高（加热过程）或降低（冷却过程），在经历相当（加热过程）或降低（冷却过程），在经历相当长时间后，物体温度逐渐趋近于周围介质温度，长时间后，物体温度逐渐趋近于周围介质温度，最终达到热平衡最终达到热平衡2022-1-65 3-1 非稳态导热过程非稳态导热过程2022-1-661 温度分布温度分布一平壁初始温度为一平壁初始温度为t0，令其左侧表面的温度突，令其左侧表面的温度突然升高到然升高到t1，右侧与温度为，右侧与温度为t0的空气接触。的空气接触。首先，物体紧挨高温表面的部分温度上升很首先，物体紧挨高温表面的部分温度上升很快，经过一定时间后内部区域温度依次变化，快

4、，经过一定时间后内部区域温度依次变化，最终整体温度分布保持恒定，当最终整体温度分布保持恒定，当 为常数时，为常数时，最终温度分布为直线。最终温度分布为直线。2022-1-67(a) = 1 (b) = 2 (c) = 3 (d) = 4非稳态导热的不同时刻物体的温度分布非稳态导热的不同时刻物体的温度分布2022-1-682 两个阶段：两个阶段：非正规状况阶段（初非正规状况阶段（初始状况阶段）、正规状况阶段始状况阶段）、正规状况阶段非正规状况阶段（初始状况阶段）：非正规状况阶段（初始状况阶段）：在在 = 3时刻之前的阶段，物体内时刻之前的阶段，物体内的温度分布受初始温度分布的影响的温度分布受初始

5、温度分布的影响较大。较大。必须用无穷级数描述必须用无穷级数描述t0t1 = 3正规状况阶段：在正规状况阶段：在 = 3时刻之后，初始温度时刻之后，初始温度分布的影响已经消失，物体内的温度分布主分布的影响已经消失，物体内的温度分布主要受边界条件的影响，要受边界条件的影响，可以用初等函数描述。可以用初等函数描述。2022-1-693 热量变化热量变化：与稳态导热的另一区别：由于有温度变化要与稳态导热的另一区别：由于有温度变化要积聚或消耗热量，同一时刻流过不同界面的积聚或消耗热量，同一时刻流过不同界面的热流量是不同的。通过截面热流量是不同的。通过截面A的热流量是从的热流量是从最高值不断减小，在其它各

6、截面的温度开始最高值不断减小，在其它各截面的温度开始升高之前通过此截面的热流量是零，温度开升高之前通过此截面的热流量是零，温度开始升高之后，热流量才开始增加。始升高之后，热流量才开始增加。2022-1-610(a) = 1 (b) = 2 (c) = 3 (d) = 42022-1-6114 边界条件对温度分布的影响边界条件对温度分布的影响x0 xtt(b)(a)(c)环境（边界条件）对系统温度环境（边界条件）对系统温度分布的影响是很显著的分布的影响是很显著的,这里以这里以一维非稳态导热过程（也就是一维非稳态导热过程（也就是大平板的加热或冷却过程）为大平板的加热或冷却过程）为例来加以说明。例来

7、加以说明。图表示一个大平板的加热过程，图表示一个大平板的加热过程，并画出在某一时刻的三种不同并画出在某一时刻的三种不同边界情况的温度分布曲线（边界情况的温度分布曲线（a）、）、（b）、（）、（c）2022-1-612这实质上是表明在第三这实质上是表明在第三类边界条件下可能的三类边界条件下可能的三种温度分布。种温度分布。 按照传热关系式按照传热关系式作一个近似的分析。作一个近似的分析。 tthttqww1x0 xtt(b)(a)(c)2022-1-613曲线（曲线（a）表示平板外环境）表示平板外环境的换热热阻的换热热阻 远大于平远大于平板内的导热热阻板内的导热热阻 , 即即 h1/1h从曲线上看

8、，物体内部的温度从曲线上看，物体内部的温度几乎是均匀的，这也就说物体几乎是均匀的，这也就说物体的温度场仅仅是时间的函数，的温度场仅仅是时间的函数，而与空间坐标无关。我们称这而与空间坐标无关。我们称这样的非稳态导热系统为集总参样的非稳态导热系统为集总参数系统（一个等温系统或物数系统（一个等温系统或物体）。体）。 x0 xtt(b)(a)(c)2022-1-614曲线（曲线（b）表示平板外环）表示平板外环境的换热热阻境的换热热阻 相当于相当于平板内的导热热阻平板内的导热热阻 , 即即 h1/1h这也是正常的第三类边界这也是正常的第三类边界条件条件 x0 xtt(b)(a)(c)2022-1-615

9、曲线（曲线（c）表示平板外环）表示平板外环境的换热热阻境的换热热阻 远小于远小于平板内的导热热阻平板内的导热热阻 , 即即 h1/1h从曲线上看，物体内部温从曲线上看，物体内部温度变化比较大，而环境与度变化比较大，而环境与物体边界几乎无温差，此物体边界几乎无温差，此时可用认为时可用认为 。那么，。那么，边界条件就变成了第一类边界条件就变成了第一类边界条件，即给定物体边边界条件，即给定物体边界上的温度。界上的温度。 wttx0 xtt(b)(a)(c)2022-1-616t0t t01/h / 2 t / 1/h2 t 1/h / 2 t02022-1-617把导热热阻与换热热阻相比可得到一个无

10、把导热热阻与换热热阻相比可得到一个无因次的数，我们称之为毕欧（因次的数，我们称之为毕欧（Boit）数，）数，即即 那么，上述三种情况则对应着那么，上述三种情况则对应着Bi1。hhBi1毕欧数是导热分析中的一个重要的无因次毕欧数是导热分析中的一个重要的无因次准则，它准则，它表征了给定导热系统内的导热热表征了给定导热系统内的导热热阻与其和环境之间的换热热阻的对比关系阻与其和环境之间的换热热阻的对比关系。2022-1-618类似于类似于Bi数这种表征某一类物理现象或物数这种表征某一类物理现象或物体特征的无量纲数称为特征数，特征数中体特征的无量纲数称为特征数，特征数中的几何尺度称为的几何尺度称为特征尺

11、度特征尺度。t iBiB00Bi1223121201010000tt 0tt 0tt 2022-1-619 3-1 非稳态导热过程非稳态导热过程3-2 集总参数法集总参数法2022-1-6201 定义定义忽略物体内部导热热阻、认为物体温度均匀一忽略物体内部导热热阻、认为物体温度均匀一致的分析方法。此时，致的分析方法。此时， ，温度分布只，温度分布只与时间有关，即与时间有关，即 ，与空间位置无关，与空间位置无关，因此，也称为因此，也称为零维零维问题。问题。0Bi)(ft hhBi/1/2022-1-621以下几种情况以下几种情况Bi很小，可用集总参数法：很小，可用集总参数法：（1）导热系数相当大

12、；）导热系数相当大；（2）几何尺寸很小；）几何尺寸很小；（3）表面换热系数很小。）表面换热系数很小。2 温度分布温度分布hhBi/1/h, t AQc, c, V, t0一个集总参数系统，其体积为一个集总参数系统，其体积为V、表面积为、表面积为A、密度为密度为 、比热为、比热为c以及初始温度为以及初始温度为t0，突然放，突然放入温度为入温度为t 、换热系数为、换热系数为h的环境中。的环境中。 2022-1-622)(tthAddtcV引入过余温度：引入过余温度：tthAddcV初始条件为：初始条件为：00)0(ttdcVhAd能量守恒能量守恒：单位时间物体热力学能的变化量：单位时间物体热力学能

13、的变化量应该等于物体表面与流体之间的对流换热量应该等于物体表面与流体之间的对流换热量, t0AQc, c, V, t02022-1-623积分得：积分得：dcVhAd00cVhA0lncVhAexp0dcVhAd方程中指数的量纲：方程中指数的量纲：sJWmKkgJmkgmKmWVchA133222022-1-6243 时间常数时间常数 hAcV称为系统的时间常数，记为称为系统的时间常数，记为 s s，也称，也称弛弛豫时间。豫时间。 如果导热体的热容量（如果导热体的热容量（ Vc ）小、换）小、换热条件好（热条件好（hA大），那么单位时间所大），那么单位时间所传递的热量大、导热体的温度变化快，传

14、递的热量大、导热体的温度变化快，时间常数时间常数 ( Vc / h A) 小小2022-1-625hAcV反映了系统处于一定的环境中所表现出来的反映了系统处于一定的环境中所表现出来的传热动态特征，与其几何形状、密度及比热传热动态特征，与其几何形状、密度及比热有关，还与环境的换热情况相关。有关，还与环境的换热情况相关。可见，同一物质不同的形状其时间常数不同，可见，同一物质不同的形状其时间常数不同，同一物体在不同的环境下时间常数也是不相同一物体在不同的环境下时间常数也是不相同。同。 2022-1-626如图所示，时间常数越小，物体的温度变化就如图所示，时间常数越小，物体的温度变化就越快，物体就越迅

15、速地接近周围流体的温度。越快，物体就越迅速地接近周围流体的温度。这说明，时间常数反映物体对环境温度变化响这说明，时间常数反映物体对环境温度变化响应的快慢，时间常数小的响应快，时间常数大应的快慢，时间常数小的响应快，时间常数大的响应慢。的响应慢。用热电偶测量流体温度，总用热电偶测量流体温度，总是希望热电偶的时间常数越是希望热电偶的时间常数越小越好，时间常数越小，热小越好，时间常数越小，热电偶越能迅速地反映流体的电偶越能迅速地反映流体的温度变化，故热电偶端部的温度变化，故热电偶端部的接点总是做得很小接点总是做得很小2022-1-627/0/s0.368101当物体冷却或加热过程当物体冷却或加热过程

16、所经历的时间等于其时所经历的时间等于其时间常数时，即间常数时，即 =s，368. 010e=4s，01. 00 . 40e工程上认为工程上认为 = 4s时导时导热体已达到热平衡状态热体已达到热平衡状态2022-1-6284 瞬态热流量瞬态热流量导热体在时间导热体在时间 0 内传给流体的总热量：内传给流体的总热量：)(tthA)exp()(0VVFoBitthAdcVhAAhttdQexp)(000cVhAcVttexp1)(02022-1-629如何去判定一个任意的系统是集总参数系统如何去判定一个任意的系统是集总参数系统 ？VVFoBiAVAVhcVhAeee2)(a)(0V/A具有长度的因次

17、，称为集总参数系统的具有长度的因次，称为集总参数系统的特征尺寸特征尺寸。MBiV1 . 0为判定系统是否为集总参数系为判定系统是否为集总参数系统统 ，M为形状修正系数。为形状修正系数。 5 集总参数系统的判定集总参数系统的判定 2022-1-630厚度为厚度为2 的大平板的大平板V A1M直径为直径为2r的长圆柱体的长圆柱体V Ar 25 . 0M直径为直径为2r的球体的球体3rAV31M复杂形体复杂形体0333. 0VBi31M形状修正系数形状修正系数M2022-1-631m10953. 0)m001. 0m02. 0(2m02. 0m002. 0)5 . 0(22322RlRlRRllRA

18、V05. 01015. 136.1010953. 05 .12)/(33AVhBiV例：一温度计水银泡是圆柱形，长例：一温度计水银泡是圆柱形，长20mm，内，内径径4mm，测量气体温度，表面传热系数，测量气体温度，表面传热系数h=12.5W/(m2K),若要温度计的温度与气体的若要温度计的温度与气体的温度之差小于初始过余温度的温度之差小于初始过余温度的10%，求测温所，求测温所需要的时间。水银需要的时间。水银 =10.36 W/(mK), = 13110 kg/m3, c = 0.138 kJ/(kgK).解：解：2022-1-632故可以用集总参数法。故可以用集总参数法。%10)1015.

19、1exp()exp(300FoFoBittttVV25.2002)10953. 0(1311010138. 036.10)/(2332AVc由上式解得：由上式解得： = 333 s = 5.6 min为了减小测温误差，测温时间应尽量加长。为了减小测温误差，测温时间应尽量加长。2022-1-633一厨师在炒鸡肉丝时要品尝一下一厨师在炒鸡肉丝时要品尝一下咸淡，于是他从咸淡，于是他从100的热炒锅的热炒锅中取出一鸡肉丝，用口吹了一会，中取出一鸡肉丝，用口吹了一会，待其降至待其降至65时再放入口中。试时再放入口中。试估算厨师需要吹多长时间？估算厨师需要吹多长时间？出锅时鸡肉丝可视为平均直径为出锅时鸡肉

20、丝可视为平均直径为2mm的圆条，的圆条，厨师口中吹出的气流温度为厨师口中吹出的气流温度为30，其与鸡肉，其与鸡肉丝之间的表面传热系数为丝之间的表面传热系数为100W/m2K，鸡肉，鸡肉丝的丝的 = 810 kg/m3，c = 3.35kJ/(kg)， = 1.1 W/(m K)。2022-1-634解：首先检验是否可用集总参数法。为此计解：首先检验是否可用集总参数法。为此计算算BiV，故可以采用集总参数法。故可以采用集总参数法。05. 0045. 02)2()/(2hrrllrhAVhBiv10705. 02srchcVhA5 . 01003010065tttto)0705. 0exp(5 .

21、 0s83. 92022-1-635例：空气流过球表面的换热系数，可采用观例：空气流过球表面的换热系数，可采用观察一个纯铜制成的球的温度随时间的变化求察一个纯铜制成的球的温度随时间的变化求得。得。球的直径为球的直径为12.7mm，在把它放入温度为，在把它放入温度为27的气流中之前，温度为的气流中之前，温度为66。当球被放。当球被放进气流中进气流中69秒钟之后，球外表面上热电偶的秒钟之后，球外表面上热电偶的指示温度是指示温度是55。计算对流换热系数。计算对流换热系数。假设：假设：1）球温是均匀的；）球温是均匀的；2）辐射交换忽略）辐射交换忽略不计；不计；3）常物性。）常物性。 2022-1-63

22、6物性：查物性参数表：纯铜（物性：查物性参数表：纯铜（66）=8933kg/m3；c=389J/（kg） =398W/(m) 分析：暂且按集总参数法处理，其温度随分析：暂且按集总参数法处理，其温度随时间变化的关系式为：时间变化的关系式为： 6expexp)(iDchcVhAt2022-1-637当=69S时，有 60127. 0389893369exp27662755h解得对流换热系数解得对流换热系数h为为 h=35.3 W/(m2)说明：以上用的是集总参数法计算得结果，因说明：以上用的是集总参数法计算得结果，因此需要检验本题是否可以用集总参数法求解：此需要检验本题是否可以用集总参数法求解：

23、3/1 . 0105 . 8/1/4RhhBi2022-1-638 3-1 非稳态导热过程非稳态导热过程3-2 集总参数法集总参数法3-3 一维非稳态导热的分析解一维非稳态导热的分析解 2022-1-639当几何形状及边界条件都比较简单时可获得分当几何形状及边界条件都比较简单时可获得分析解。析解。1 无限大的平板的分析解无限大的平板的分析解2022-1-640厚度厚度 2 的无限大平壁，的无限大平壁， 、a为已知常数为已知常数； =0时温度为时温度为 t0;突然把两侧介质温度降突然把两侧介质温度降低为低为 t 并保持不变；壁并保持不变；壁表面与介质之间的表面表面与介质之间的表面传热系数为传热系

24、数为h。两侧冷却情况相同、温两侧冷却情况相同、温度分布对称。中心为原度分布对称。中心为原点。点。2022-1-641导热微分方程：导热微分方程：22xtat初始条件：初始条件： , 00tt 边界条件：边界条件： (第三类第三类)0 , 0 xtx)(- ,tthxtx2022-1-64222xtat , 00tt )(- ,0 , 0tthxtxxtx过余温度 ),(txt22xa00 , 0 -tt0 , 0 xxxhxx- , 2022-1-643采用分离变量法求解：取采用分离变量法求解：取)()(TxX22xa00 , 0 -ttxhxxxx- ,0 , 0DdxXdXddTaT221

25、1只能为常数：只能为常数：2211dxXdXddTaT只为只为 的函数的函数只为只为x x的函数的函数2022-1-644对对 积分积分DTaT1得到得到aDeCT1式中式中C1是积分常数，常是积分常数，常数值数值D的正负可以从物的正负可以从物理概念上加以确定。理概念上加以确定。当时间当时间趋于无穷大时，过程达到稳态，趋于无穷大时，过程达到稳态，物体达到周围环境温度，所以物体达到周围环境温度，所以D必须为负必须为负值值，否则物体温度将无穷增大。，否则物体温度将无穷增大。2022-1-645令令2D21ddTaT则有则有 以及以及2221dxXdX以上两式的通解为：以上两式的通解为：21aeCT

26、)sin()cos(32xCxCX于是于是)sin()cos(),(2xBxAexa2022-1-646常数常数A、B和和可由边界条件确定。可由边界条件确定。00 , 0 -ttxhxxxx- ,0 , 0(1)(2)(3)由边界条件（由边界条件（2）得）得B=0)sin()cos(),(2xBxAexa（a）边界条件（边界条件（3）代入）代入(b) 得得 (c)htg)(（a）式成为）式成为 (b)cos(),(2xAexa2022-1-647将将 右端整理成：右端整理成：htg)(Bihhy1注意，这里注意，这里Bi数的尺度数的尺度为平板厚度的一半。为平板厚度的一半。显然，显然，是两曲线交

27、是两曲线交点对应的所有值。式点对应的所有值。式(c)称为特征方程。称为特征方程。 称为特征值。分别为称为特征值。分别为1、 2 n。2022-1-648)cos(),(11121xAexa)cos(),(22222xAexa)cos(),(2xAexnnann.1)cos(),(2nnnaxAexn将无穷个解叠加：将无穷个解叠加：至此，我们获得了无穷个特解：至此，我们获得了无穷个特解：2022-1-649利用初始条件利用初始条件 求求An00 , 0 -tt1)cos(),(2nnnaxAexn)cos()sin()sin(20nnnnnA解的最后形式为：解的最后形式为：)exp()cos()

28、cos()sin()sin(2),(21n0nnnnnnaxx令令n=n)exp()cos(cossin)sin(2),(221n0nnnnnnaxx2022-1-650)exp()cos(cossin)sin(2),(221n0nnnnnnaxx2Foa傅里叶准则傅里叶准则)exp()cos(cossin)sin(2),(21n0nnnnnnFoxx2FoahBix 无量纲距离无量纲距离) ,Fo Bi,(),( 0 xfx2022-1-651定义无量纲的热量定义无量纲的热量0QQ其中其中Q为为0时间内传导的热量（内热能的改变量）时间内传导的热量（内热能的改变量） VcQ00为为 至无穷至无

29、穷时间内的总传导热量时间内的总传导热量（物体内能改变总量）（物体内能改变总量） 设从初始时刻至某一时刻设从初始时刻至某一时刻 所传递的热量为所传递的热量为Q:)(),(000ttcVdVxttcQQVdVttttttVV00)(101dVVV011dVttVV)(1是是 时刻物体的平均过余温度。时刻物体的平均过余温度。2022-1-6522 非稳态导热的正规状况阶段非稳态导热的正规状况阶段当当Fo 0.2时，采用级数的第一项计算偏差小时，采用级数的第一项计算偏差小于于1%，故当，故当Fo 0.2时时:xeFo1)(11110cos)cos()sin()sin(221其中其中 1 是第一特征值，

30、是是第一特征值，是Bi的函数。的函数。Bi0.010.050.10.51.05.0105010010.09980.22170.31110.65330.86031.31381.42891.54001.55521.5708)exp()cos(cossin)sin(2),(221n0nnnnnnaxx2022-1-653为了分析这时温度分布的特点，将式取对数得：为了分析这时温度分布的特点，将式取对数得： xFo11111210coscossinsin2ln)(ln式右边第一项是时间式右边第一项是时间 的线性函数，的线性函数， 的系数只的系数只与与Bi有关，即只取决于第三类边界条件、平壁有关，即只取决

31、于第三类边界条件、平壁的物性与几何尺寸。的物性与几何尺寸。右边第二项只与右边第二项只与Bi、x/ 有关，与时间有关，与时间 无关。无关。上式说明，当上式说明，当Fo 0.2,平壁内所有各点过余温平壁内所有各点过余温度的对数都随时间线性变化，并且变化曲线的度的对数都随时间线性变化，并且变化曲线的斜率都相等斜率都相等,这一温度变化阶段称为非稳态导热这一温度变化阶段称为非稳态导热的正规状况阶段。的正规状况阶段。2022-1-654xxm)(cos)(),(1这时比值与这时比值与 无关，仅与几何位置无关，仅与几何位置(x/ )及边界及边界条件条件(Bi数数)有关。即初始条件的影响已经消失。有关。即初始

32、条件的影响已经消失。无论初始分布如何，无量纲温度都是一样的。无论初始分布如何，无量纲温度都是一样的。这是非稳态导热的正规状态或充分发展阶段。这是非稳态导热的正规状态或充分发展阶段。 当当Fo 0.2时任一点过余温时任一点过余温度与中心过余温度之比为度与中心过余温度之比为xeFo1)(11110cos)cos()sin()sin(2212022-1-655xxm)(cos)(),(1令令x = 可以计算平壁表面温度和中心温度的可以计算平壁表面温度和中心温度的比值。比值。又由表又由表3-1可知，当可知，当Bi 0.1时，时， 1 0.95。即当。即当Bi 0.1时，平壁表时，平壁表面温度和中心温度的差别小于面温度和中心温度的差别小于5%，可以近似，可以近似认为整个平壁温度是均匀的。这就是认为整个平壁温度是均匀的。这就是3-2节集节集总参数法的界定值定为总参数法的