山西能源学院教案 授课班级能动1701-1704 授课时间 计2学时 课题（章节 及内容) 3.3典型一维非稳态导热的分析解 教学目的 理解三类边界条件下：无限大平板、无限长圆柱、球的分析解： 和要求 了解非稳态导热正规状况阶段分析解的简化： 了解诺谟图的意义。 重点 难点 三类边界条件下无限大平板的分析解。 教学进程 (含课堂 教学内容：三种几何形状物体的温度场分析解，非稳态导热正规 教学内容、 状况阶段分析解的简化：非稳态导热正规状况阶段的工程计算方 教学方法、 法。 辅助手段等) 教学方法：讲授与练习、启发讨论、诱导式、归纳总结法。 作业布置 主要 《传热学》第四版，杨世铭，陶文铨， 参考资料 高等教育出版社，2006年8月 课后自我 总结分析
山西能源学院教案 授课班级 能动 1701-1704 授课时间 计 2 学时 课题（章节 及内容） 3.3 典型一维非稳态导热的分析解 教学目的 和要求 理解三类边界条件下：无限大平板、无限长圆柱、球的分析解； 了解非稳态导热正规状况阶段分析解的简化； 了解诺谟图的意义。 重 点 难 点 三类边界条件下无限大平板的分析解。 教学进程 （含课堂 教学内容、 教学方法、 辅助手段等） 教学内容：三种几何形状物体的温度场分析解，非稳态导热正规 状况阶段分析解的简化；非稳态导热正规状况阶段的工程计算方 法。 教学方法：讲授与练习、启发讨论、诱导式、归纳总结法。 作业布置 主 要 参考资料 《传热学》第四版，杨世铭，陶文铨， 高等教育出版社，2006 年 8 月 课后自我 总结分析
山西能源学院教案 3-3典型一维非稳态导热的分析解 本节介绍第三类边界条件下：无限大平板、无限长圆柱、球的分析解及应 用。如何理解无限大物体，如：当一块平板的长度、宽度>厚度时，平板 的长度和宽度的边缘向四周的散热对平板内的温度分布影响很少，以至于可 以把平板内各点的温度看作仅是厚度的函数时，该平板就是一块“无限大” 平板。若平板的长度、宽度、厚度相差较小，但平板四周绝热良好，则热量 交换仅发生在平板两侧面，从传热的角度分析，可简化成一维导热问题。 一、无限大平板的分析解 已知：厚度2ó的无限大平板，初温o,初始瞬间将其放于温度为的 流体中，而且。>6，流体与板面间的表面传热系数为一常数。 试确定在非稳态过程中板内的温度分布。 解：如图3-5所示，平板两面对称受热，所以其内温度分布以其中心截面为对 称面。 对于x20的半块平板，其导热微分方程： 2 (00) (3-8) 定解条件：tx,0)=to(0≤x≤8) (x.0lx=0=0 8 (边界条件) 邮(6)-im】=-2(x x=6 (边界条件) 引入过余温度：8=(x,)-t 则 x2(00)(3-9) e(x,0)=8(0≤x≤6) (初始条件) a8(x,t x=0=0 Bx (边界条件)
山西能源学院教案 3-3 典型一维非稳态导热的分析解 本节介绍第三类边界条件下：无限大平板、无限长圆柱、球的分析解及应 用。如何理解无限大物体，如：当一块平板的长度、宽度 >> 厚度时，平板 的长度和宽度的边缘向四周的散热对平板内的温度分布影响很少，以至于可 以把平板内各点的温度看作仅是厚度的函数时，该平板就是一块“无限大” 平板。若平板的长度、宽度、厚度相差较小，但平板四周绝热良好，则热量 交换仅发生在平板两侧面，从传热的角度分析，可简化成一维导热问题。 一、 无限大平板的分析解 已知：厚度 的无限大平板，初温 t0，初始瞬间将其放于温度为 的 流体中，而且 > ，流体与板面间的表面传热系数为一常数。 试确定在非稳态过程中板内的温度分布。 解：如图 3-5 所示，平板两面对称受热，所以其内温度分布以其中心截面为对 称面。 对于 x 0 的半块平板，其导热微分方程： (0<x< δ, ) （ 3-8 ） 定解条件：t(x,0)= t0(0 x δ) （边界条件） （边界条件） 引入过余温度： 则 （ 0<x< δ, ） (3-9) (x,0)= (0 x δ) （初始条件） （边界条件）
h(6,)=- ae(x,t) 8x x=6 (边界条件) ae a2e -=d 对偏微分方程 a 0x2 分离变量求解得： (x,) a.0'85 =22 sm(Aeo9(A,剖 品 2=1 6,6+sin(月n6)cos(66) (3-10) 其中离散值A是下列超越方程的根，称为特征值。 B tan(B6)= 82 ，2=1,2, 年888t (3-11) 其中Bi是以特征长度为δ的毕渥数。 由此可见：平板中的无量纲过余温度 与三个无量纲数有关：以平板厚度 半δ为特征长度的傅立叶数、毕渥数及 治即： g0=,9-=fo,,爱 to-too (3-12) 二、非稳态导热的正规状况阶段 1、平板中任一点的过余温度与平板中心的过余温度的关系 前述得到的分析解是一个无穷级数，计算工作量大，但对比计算表明，当 F0>0.2时，采用该级数的第一项与采用完整的级数计算平板中心温度的误差小 于1%，因此，当F0>02时，采用以下简化结果： 8(x,x= 8+n月，0xA,eia《B)青] 2sin(318) (3-13) 其中特征值6(2=1,2)之值与Bi有关。 由上式(3-13)可知：Fo>0.2以后平板中任一点的过余温度8(x，t) (x, 与平板中心的过余温度0，t)=8m(T)之比为：8m(① =cos(u) 61 (3-14)
（边界条件） 对偏微分方程 分离变量求解得： （ 3-10 ） 其中离散值 是下列超越方程的根，称为特征值。 …… （ 3-11 ） 其中 Bi 是以特征长度为 δ 的毕渥数。 由此可见：平板中的无量纲过余温度 与三个无量纲数有关：以平板厚度 一 半 δ 为 特 征 长 度 的 傅 立 叶 数 、 毕 渥 数 及 即 ： （ 3-12 ） 二、非稳态导热的正规状况阶段 1 、平板中任一点的过余温度与平板中心的过余温度的关系 前述得到的分析解是一个无穷级数，计算工作量大，但对比计算表明，当 Fo>0.2 时，采用该级数的第一项与采用完整的级数计算平板中心温度的误差小 于 1% ，因此，当 Fo>0.2 时，采用以下简化结果： （ 3-13 ） 其中特征值 之值与 Bi 有关。 由上式（ 3-13 ）可知： Fo>0.2 以后平板中任一点的过余温度 (x ， τ) 与平板中心的过余温度 (0 ， τ)= （ τ ）之比为： （ 3-14 ）
此式反映了非稳态导热过程中一种很重要的物理现象：即当F0>0.2以后， 虽然(x,T)与8m(T)各自均与t有关，但其比值则与t无关，而仅取决 于几何位置（6)及边界条件(B)。也就是说，初始条件的影响已经消失， 无论初始条件分布如何，只要Fo>0.2,(x,之值是一个常数，也就是无 0() 量纲的温度分布是一样的。 由此可见，当Fo>0.2时，非稳态导热过程进入正规状况阶段。 2、在一个时间间隔内非稳态导热过程中传递的热量 1)从物体初始时刻平板与周围介质处于热平衡，这一过程中传递的热量： 20=V(to-too) (3-15) 此值为非稳态导热过程中传递的最大热量。 2)从初始时刻到某一时间τ，这段时间内所传递的热量0→： Q→r=a[-(x,W (3-16) 90→x 3）之比： 答0r 0-too =1- 1 (3-17) 其中：a=(是时刻T物体的平均过余温度， 8-ll-todv 对于无限大平板，当Fo>0.2,将式(3-13)代入的定义式，可得： 2sm4—g(4网)im4 4+sin4cos玛 巧 (3-18) 对圆柱体、球体 2>0.2时，无穷级数的解也可用第一项近似代替
此式反映了非稳态导热过程中一种很重要的物理现象：即当 Fo>0.2 以后， 虽然 (x ，τ) 与 （τ）各自均与τ有关，但其比值则与τ无关，而仅取决 于几何位置（x ）及边界条件（ Bi ）。也就是说，初始条件的影响已经消失， 无论初始条件分布如何，只要 Fo>0.2 ， ( ) ( , ) m x 之值是一个常数，也就是无 量纲的温度分布是一样的。 由此可见，当 Fo>0.2 时，非稳态导热过程进入正规状况阶段。 2 、在一个时间间隔内非稳态导热过程中传递的热量 1 ） 从物体初始时刻平板与周围介质处于热平衡，这一过程中传递的热量： （ 3-15 ） 此值为非稳态导热过程中传递的最大热量。 2 ） 从初始时刻到某一时间 τ ，这段时间内所传递的热量 ： （ 3-16 ） 3 ） 之比： ( 3-17 ) 其中： 是时刻 τ 物体的平均过余温度， 。 对于无限大平板，当 Fo>0.2 ，将式（ 3-13 ）代入 的定义式，可得： （ 3-18 ） 对圆柱体、球体 >0.2 时，无穷级数的解也可用第一项近似代替
8(x,) 并且品 及武)可表示为： (x,=Ae即(-4好Fo)f(47) (3-19) e=Aop(-Li Ro)B (x)/ (3-30) 其中：n为无量纲几何位置，对平板刀=x16,对柱体及球体刀=r1R,R 为外表面半径，系数A、B及函数()的表达式取决于几何形状，见教材 表3-2所示。 三、正规阶段状况的实用计算方法 当Fo>0.2时，可采用上述计算公式求得非稳态导热物体的温度场及交换的 热量，也可采用简化的拟合公式和诺谟图求得。 1、诺模图：工程技术中，为便于计算，采用按分析解的级数第一项绘制的一些 图线，叫诺谟图。 2、海斯勒图：诺漠图中用以确定温度分布的图线，称海斯勒图。 首先根据（3一13)式给出8m/随Fo及Bi变化的曲线（此时x/6 =0),然后根据(3一14)式确定8/8m的值，于是平板中任意一点的日/ 日_8m8 值便为： 668m (3-21) 同样，从初始时刻到时刻τ物体与环境间所交换的热量，可采用(3一 15)、(3一17)作出20 =f(Fo.Bi) 曲线。 3、诺漠图法评述 优点：简洁方便。 缺点：准确度有限，误差较大。 目前，随着计算技术的发展，直接应用分析解及简化拟合公式计算的方法受 到重视
并且 及 可表示为： （ 3-19 ） （ 3-30 ） 其中： η 为无量纲几何位置，对平板 ，对柱体及球体 ， R 为外表面半径，系数 A 、 B 及函数 的表达式取决于几何形状，见教材 表 3-2 所示。 三、正规阶段状况的实用计算方法 当 Fo>0.2 时，可采用上述计算公式求得非稳态导热物体的温度场及交换的 热量，也可采用简化的拟合公式和诺谟图求得。 1 、诺模图：工程技术中，为便于计算，采用按分析解的级数第一项绘制的一些 图线，叫诺谟图。 2 、海斯勒图：诺谟图中用以确定温度分布的图线，称海斯勒图。 首先根据（ 3—13 ）式给出 随 Fo 及 Bi 变化的曲线（此时 x/δ =0 ），然后根据（ 3 — 14 ）式确定 的值，于是平板中任意一点的 值便为： （ 3-21 ） 同样，从初始时刻到时刻 τ 物体与环境间所交换的热量，可采用（ 3 — 15 ）、（ 3 — 17 ）作出 曲线。 3 、诺谟图法评述 优点：简洁方便。 缺点：准确度有限，误差较大。 目前，随着计算技术的发展，直接应用分析解及简化拟合公式计算的方法受 到重视
四、分析解应用范围的推广及讨论 1、推广范围 1)对物体被冷却的情况也适用： 2)也适于一侧绝热，另一侧为第三类边界条件的厚为δ的平板： 3)当固体表面与流体间的表面传热系数h→∞时，即表面换热热阻→0时， 所以五→时分析解就是固体表面温度发生一突然变化然后保持不变时的解， 即第一类边界条件的解。 2、讨论Bi与Fo对温度场的影响： 1)傅立叶数Fo: 由(3-10)、(3-13)式及诺模图可知：物体中各点的过余温度随时间τ的 增加而减小：而Fo与成正比，所以物体中各点过余温度亦随Fo的增大而 减小。 2)毕渥数Bi Bi对温度的影响从以下两方面分析： 一方面，从教材图3一6可知，F0相同时，Bi越大，8m/品越小。因 为，B越大，意味着固体表面的换热条件越强，导致物体的中心温度越迅速地 接近周围介质的温度：当Bi→∞时，意味着在过程开始瞬间物体表面温度就达 到介质温度，物体中心温度变化最快，所以在诺模图中1B=0时的线就是壁面 温度保持恒定的第一类边界条件的解。 另一方面B的大小决定于物体内部温度的扯平程度。如：对于平板，从诺 模图3一7中可知： 1 当>10（即B<0.1)时，截面上的过余温度差小于5% 当Bi下限一直推到0.01时，其分析解与集总参数法的解相差极微。 综上可得如下结论：介质温度恒定的第三类边界条件下的分析解：当B →∞时，转化为第一类边界条件下的解，B→0时，则与集总参数法的解相同
四、分析解应用范围的推广及讨论 1 、推广范围 1 ）对物体被冷却的情况也适用； 2 ）也适于一侧绝热，另一侧为第三类边界条件的厚为 δ 的平板； 3 ）当固体表面与流体间的表面传热系数 h 时，即表面换热热阻 0 时， 所以 时分析解就是固体表面温度发生一突然变化然后保持不变时的解， 即第一类边界条件的解。 2 、讨论 Bi 与 Fo 对温度场的影响： 1 ）傅立叶数 Fo ： 由 (3-10) 、 (3-13) 式及诺模图可知：物体中各点的过余温度随时间 τ 的 增加而减小；而 Fo 与 成正比，所以物体中各点过余温度亦随 Fo 的增大而 减小。 2 ）毕渥数 Bi Bi 对温度的影响从以下两方面分析： 一方面，从教材图 3 — 6 可知， Fo 相同时， Bi 越大， 越小。因 为， Bi 越大，意味着固体表面的换热条件越强，导致物体的中心温度越迅速地 接近周围介质的温度；当 Bi 时，意味着在过程开始瞬间物体表面温度就达 到介质温度，物体中心温度变化最快，所以在诺模图中 1/Bi=0 时的线就是壁面 温度保持恒定的第一类边界条件的解。 另一方面 Bi 的大小决定于物体内部温度的扯平程度。如：对于平板，从诺 模图 3 — 7 中可知： 当 >10 （即 Bi<0.1 ）时，截面上的过余温度差小于 5 % 当 Bi 下限一直推到 0.01 时，其分析解与集总参数法的解相差极微。 综上可得如下结论：介质温度恒定的第三类边界条件下的分析解；当 Bi 时，转化为第一类边界条件下的解， Bi 0 时，则与集总参数法的解相同