二维和三维不稳定态导热的时间计算.pdfVIP

年 月 了 钢 铁 研 究 乙 第 期 总第 期 二维和三维不稳定态导热的时间计算 王 春 霖 武钢职工大学 邮编 【摘要】 荃 于一维不稳态导热分析解 , 提出一个在第三类边 界条件下二维和三维不稳态导热 二、 。 , , , 一 二 。二 、李 , 、一二 声 。吃刃 月 二、 二 ‘ 二 斗 “ 、 二 , 廿 二 的时间计算方法 。 其中 , 在求超越方程 苦 一 群和举召头 一 合 的根时 , 应用牛顿插值法求解甚为“ “ ’一 , 犷 ’ 汁 尹 ’ 执 “ 片 ” 协 , 、 彬喝 ” 一 “ 尸 ‘” 月 “ , , “ 一 ’ 一 ’ ’ ” , ‘川 , “ “ 、 一 方便 。 关链词 二维和三维不稳态导热 牛顿插值 ! 入 ! ? ! ! ! , 一 一 一 、 ‘‘ · ‘ · , ‘ ‘ 、 · 。, ‘ 会 一 ‘ 。 ‘ ,。· ·‘ · 。 月 , 月 一 音 ‘。二 , “‘一 ,’ 。 一 一 前 言 关 于物体的不稳态加热 或冷却 下同 计算一为 已知加热时间时求解物体内的温 度分布 另为 既定物体内某点 通 常是 中心 点 温度时求解所需加热时间 。 姑且称前者为 加热的 ’‘ 直接 ” 计算 后者为 ‘’ 间接 ” 计算 。 在工程技术上 对一维 问题的直接 计算 或间接计算 都可利用其 分析 解的线算 图方 便解出 。 对于二维和 三维问题的直接计算 , 也 可用一维分析解的组 合直接求解 但是对于 间接计算加热时间 则较 费周折 例如文献 〔 ’于 介绍了一种 ’ 试凑 ” 的方法 。 然而 一 卜算公二践表 明 试凑法颇有盲 目性 尤其是对三维问题很 难经二 、 三次就能试凑成功 。 也可以 用数值法 直接计算来递推加热时间仁 , · 但是数值法运 算涉及数据汲为繁多因而不便于手算 且每 一时间步长的离散温度分布 只用作递推时间 的过渡 并非是间题 的直接答案 。 所以 , 本文 提 出一种应用一维不稳态分析的组合来计算 二维和三维问题加热时间的方法 。 基本方法和公式 传热学 已 阐明 对干 厚 短 方 柱体 三 维 和圆柱体 准二 维)等规则物体困 其在第 三类边界条件下加 热 . 且有均匀的初始温度 . 33 . 分布时 , 可利 用一维 问 题的组 合求解一般 地 , 对于三维问题其关系式 为[’, 5· 。 2 。 日 一 日(x , y , 二 , r ) 一 0. : 、 ( 二 . : ) · 即又, , , : ) · 灸(z , : ) ( l ) o r = 0 , · ( r , 了 ) 艺 ” 一 1 2 J I ( 月 , : ) 月 , : 〔J :(月 , : ) + J 若(月 , : ) 〕 · J 。 ( 月, ‘灭 )e一 人产 ( 4 ) —物体 内某点的相 对过一式中 : 日 ~ 余温度 ; 夕 : = 夕 , ( 了 · r ) - 式中 :R — 圆柱 体半径 ; t(j · 约 一 勺 to 一 tf —作为 一 口 . 一2 维问题处理时 . 某坐标方 向j(’ 一 x , y , : ) 的相对 过余温 度 , 其有 式 (2) 或式 (4 )的分析解关系式 ; t一 t仕 , y · “ , 约 — 在 : 时刻物体坐标 点 (x r , 少 , 二 ) 的 温 度 ; t(j · 丁 ) — 在 r 时刻 . 作为一维问题处 理 时某坐标 j(j 一 二 , y . : ) 的温度 ; ‘ — 加热介质温度 ; t。 — 物体初始温度 。 对于物体物性参量 为常量 (对于实际计 算 . 通常取加热温度范围内的平均值)的一维 加热 . 其在第三类边界条件下的分析解为 : 平壁 一 筹一 径 向傅立叶准数 , 一 梦一径衅 欧准数 , 月是方程式 (5)的根 . 亦有无穷多个 , 并 一一对应于式 (4)和式中的 月 , : , 其 中 J 。( 月) _ 三 (5) 不万不 一 尽 J 。 ( 夕) 一零阶贝塞尔函数 . 其展开式为 , 。 ( , ) 一 1 一 (粤) 2 长子示 (粤) 4 ‘ 气乙 ! j 乙 1 , 月 、 。 万二万二二不 戈下万 少 州卜 仁j ! ) 一 乙 ( 6 ) J . (月) 一一阶贝塞尔函数 . 其展开式 为 J 、( 召) = l 月 、 : 贬1 — 又下 几万 少 川片乙! 乙 Z si n 从毖一一生里卫全 2 一, 一 1 产,。 一卜 Sln 产, C 0 5 产 co s (拜 和 。一 ‘% (2) 式中 :S — (热流方向 )平 壁厚度 ;当两侧对 称加热时为l/2厚度 ; — 傅立 叶准数 (a 为导温 系 数) ; 竺了 一 F 产 为