非稳态导热的分析计算3-2.pptVIP

非稳态导热Transient Conduction 热工教研室 王晓墨 第三章 非稳态导热的分析计算 §3.1非稳态导热过程分析 §3.2一维非稳态导热过程分析 §3.3多维非稳态导热的图解法 §3.4对分析解的讨论 §3.5集总参数系统分析 二、无限长圆柱体和球体的加热（冷却）过程分析 三、半无限大固体的非稳态导热过程 半无限大系统指的是一个半无限大的空间，也就是一个从其表面可以向其深度方向无限延展的物体系统。 很多实际的物体在加热或冷却过程的初期都可以视为是一个半无限大固体的非稳态导热过程。 §3.3多维非稳态导热的图解法 §3.5集总参数系统分析 Lumped capacitance method 1.无限长圆柱 t r t∞ t∞ t0 α α 0 式中r0 为无限长圆柱体的半径 类似有 ： 和 附录P412图1、2、3。 2.球体 t t∞ α r t0 0 球体处理方法与无限大圆柱体完全相同，相应的线算图示于P414图4、图5和图6之中。 这里要注意的是特征尺寸R为球体的半径，r为球体的径向方向。 x 0 t0 tw 该微分方程的初、边值问题可以用拉普拉斯变换求解，得到的温度分布为： 式中高斯误差函数定义 而 由傅里叶定律，任意位置上的热流量为： 显然，边界表面上的热流量为 : 应用上面讨论的海斯勒线算图可以求出厚度为2?的大平板、半径为R的无限长圆柱体、及半径为R的球体的温度分布和传导的热量。 对非一维非稳态导热问题，我们能不能利用上面的一维非稳态导热线算图来进行求解呢？ 用一个无限长矩形柱为例来回答这一问题。 0 x y 2δ1 2δ2 一个无限长矩形柱，可以看成是由两个无限大平板正交而组成，它们的厚度分别为2?1和2?2。 无限长矩形柱的导热微分方程式为： 假定, 将其代入微分方程中 一个二维非稳态导热问题的解可以用两个导热方向相互垂直的一维非稳态导热问题解的乘积来表示。 同理，一个三维非稳态导热问题的解可以用三个相互垂直的一维非稳态导热问题解的乘积来表示。 2δ2 y x 0 2δ1 例如：1.矩形截面的长棱柱（正四棱柱）：可由两个大平板正交构成，因而温度分布为两个大平板对应的温度分布的乘积 y z x 2.矩形块体(立方体) 可由三个大平板正交构成，因而温度分布为三个大平板对应的温度分布的乘积 2δ x r R 0 3.短圆柱体可由一个长圆柱体和一个大平板正交构成，因而温度分布为一个长圆柱体和一个大平板对应的温度分布的乘积 r x 0 4.半长圆柱体可由一个长圆柱体和一个半无限大固体正交构成，因而温度分布为一个长圆柱体和一个半无限大固体对应的温度分布的乘积 需要强调的是，我们要确定某一点的温度时，一定要首先确定该点在对应的几个一维空间上的位置，再去确定相应的一维温度值，最终乘积得出物体在该点的温度值。 例题3-1 一块厚200mm的大钢板，钢材的密度为ρ=7790kg/m3，比热容cp=170J/(kg·K)，导热系数为43.2W/(m·K)，钢板的初始温度为20℃，放入1000℃的加热炉中加热，表面传热系数为 h=300W/(m2·K)。试求加热40分钟时钢板的中心温度。 解： 根据题意，δ=100mm = 0.1m。钢材的热扩散率为 傅里叶数为 毕渥数为 查图可得 §3.4对分析解的讨论 1. Fo准则对温度分布的影响 Fo?0.2时，进入正规状况阶段，平壁内所有各点过余温度的对数都随时间按线性规律变化，变化曲线的斜率都相等。 θ m/θ0随F0增大而见小。 Fo0.2时是瞬态温度变化的初始阶段，各点温度变化速率不同 2. Bi准则对温度分布的影响 Bi （Bi=α? /? ）表征了给定导热系统内的导热热阻与其和环境之间的换热热阻的对比关系 。 当 Bi?? 时，意味着表面传热系数 α ?? ，对流换热热阻趋于0。平壁的表面温度几乎从冷却过程一开始，就立刻降到流体温度 t? 。 当Bi?0时，意味着物体的热导率很大、导热热阻? 0（Bi= α ?/? ）。物体内的温度分布趋于均匀一致。 可用集总参数法求解. x 0 x t t∞ (b) (a) (c) 物体内热阻可以忽略，也就是 的导热系统称为集总参数系统，有时也称为充分搅拌系统或热薄物体系统 任意形状的物体假设 Bi0.1， 表明物体内部温度分布几乎趋于一致，可以近似认为物体是一个集总参数系统。 一、集总系统的能量平衡方程和温度分布