第1章传导传热(免费阅读).pptVIP

图1-15平板在恒温介质中加热时的函数 图1-16平板在恒温介质中加热时的函数 称为毕欧准数，为物体内部的导热热阻与物体外部的换热热阻之比. Bi值较小说明物体内热阻较小，外热阻较大，物体断面温差较小；当Bi0.1时，物体内温度趋于一致，可忽略物体断面温差，此时物体的加热称为“薄材”加热，“薄材”加热不必使用三类边界条件下求解。 Bi值较大说明物体内热阻大，外热阻小，物体断面温差大；当Bi 0.1时，物体内断面温差较大，此时物体的加热称为“厚材”加热，“厚材”加热问题需采用三类边界条件求解。 例题2-4 直径为500mm的碳钢圆锭，在炉温为1250℃的恒温炉内均匀加热，问中心温度达到1180℃时所需加热时间和这时钢锭表面的温度。已知钢锭初始温度t0=20℃，平均给热系数α∑=350W?(℃?m2)-1，钢的导热系数λ=32.4 W?(m?℃)-1，导温系数a=0.022m2?h-1。 解： 查图1-18得到Fo=1.15，则 根据Bi=2.7和Fo=1.15，查图1-17得中心相对温度 可得加热3.3h后表面温度为 1.4 导热问题的数值解法 数学分析法的优点是可以直接给出温度场的分析解，便于分析各种因素对温度分布的影响，但缺点是仅限于经典问题的近似求解，即主要用于求解几何形状和边界条件比较简单的问题，实际情况下由于假设条件过多会引起较大误差。 数值解法计算简单，结果更符合实际要求，但计算工作量大。 常用的方法包括有限差分法、有限元法等 . 1.4.1 二维稳定态导热问题的差分解法 有限差分解法是将连续变化的物理过程用不连续的阶跃过程代替，把解微分方程变为解差分方程。即把连续的定解区域用有限个离散点构成的网格来代替，这些离散点称作网格节点；把定解区域上连续变量的函数用网格上定义的离散变量函数来近似；把原方程和定解条件中的微商用差商来近似，于是原微分方程和定解条件就近似地以代数方程组代之，此即有限差分方程组，解此方程组就可以得到原问题在离散点上的近似解，最后再利用插值方法从离散解得到定解问题在整个区域上的近似解。 将二维温度场（z向无温度变化）按等距离分成若干网格（离散化），网格的交叉点为节点。从中取出一部分，得到若干相等的单元体0、1、2、3、4，单元体的边长分别为△x和△y，厚度均为△z，节点上相应温度分别为t0、t1、t2、t3、t4。稳定态导热情况下，温度不随时间而变，即传热过程中流向任意节点的热量总和为零。 根据热量平衡关系，对单元体0得到 根据傅立叶导热定律，上式表达成差分形式为 若划分网格时使△x=△y，则上述差分方程简化为 设所分析的物体共有n个节点，同理可列出n个节点的差分方程式，然后联立求解，可得所有节点的温度值。张弛法是求解联立方程的方法之一，网格划分得越细，结果就越接近实际的温度分布值。 例题2-4 设炉墙转角截面如图1-20所示，内外壁温度各保持400℃和100℃，求炉墙中心的温度分布。 解 取单元边长为墙厚的一半，确立节点1、2、3、…，墙角两边对称，故t4=t4a。 400 100 先假定各节点温度为t1=t2=t3=t4=t4a=250℃，t5=200℃（节点5两边散热，温度必低于其它各点），将假设温度值代入式（1-32）， 求出各节点的余数Q′，得到 由于节点5的余数最大，为-100℃，于是，第二次迭代时将节点5的温度设定为 t5=200+(-100/4)=175℃，将节点4的温度设定为 t4=250+(-50/4)=237℃，将此值重新代入式（1-32），得到节点5的余数最大，为-26℃，再重新设定节点5的温度为t5=175+(-26/4)≈168℃，重新计算，直到达到精度要求为止，如表1-1所示 1.4.2 一维不稳定态导热问题的差分解法 对一维不稳定态导热问题，用有限个微小时间段内的稳定态等温过程代替连续变化的不稳定态导热过程，可得不稳定态导热过程的解。解的精度取决于差分值△τ、△x的值，即一般情况下△τ、△x的值取的越小，精度越高。 以平板为例用差商代替微商建立有限差分方程 如图，将一维导热大平板沿x方向以等间距△x分成N段；时间从τ=0开始，等间隔△τ分成M段。用i表示内节点坐标x的位置，用k表示第k?△τ时刻，则 表示第k?△τ时刻节点i的温度。 对内部节点：用一阶向前差分表示内节点（i, k）的温度对时间的偏导数 用二阶中心差分表示内节点（i, k）的温度对坐标x的偏导数 代入一维导热微分方程，可得 整理后得内部节点差分方程 令 上式简化成 一维不稳定态导热问题内节点的差分方程 ，具有显函数特征