D0I:10.13374/i.issm1001053x.2003.03.043 第25卷第3期 北京科技大学学报 Vol.25 No.3 2003年6月 Journal of University of Science and Technology Beijing Jun.2003 全氢罩式退火炉退火热过程的研究① 一对流换热系数和钢卷径向等效导热系数的分析 林林”张欣欣”左燚”向顺华”刘铁树)钮龙英) 黄夏兰) 1)北京科技大学机械工程学院,北京1000832)宝钢研究院设备研究所,上海201900 摘要介绍影响全氢罩式退火炉内换热的两个重要参数一对流换热系数和钢卷径向等 效导热系数,详尽分析了两个参数的影响因素,并对比了氮气和氢气气氛下两参数的不同, 从而在机理上阐明了全氢罩式炉相对传统混氢罩式炉的优越性,为优化炉内换热提供了理 论的依据. 关键词全氢罩式退火炉;对流换热系数:等效导热系数 分类号TG155;T℉061.2:T℉062 对流换热系数和钢卷径向等效导热系数是 400 影响全氢罩式退火炉(以下简称全氢炉)退火钢 350 300 卷传热的两个最主要参数,对这两个参数的研究 250 有助于了解炉内换热的机理,为全氢炉内传热的 200 强化提供理论的依据. 150 100 1对流换热系数 501 全氢炉内罩内的换热以对流换热为主,图1 0 10 2030405060 t/h 给出了全氢炉对流换热的热流密度与辐射换热 图1炉内对流换热与辐射换热热流密度之比随退火过 的热流密度之比 程的变化 由图可知,在整个退火过程开始和将近结束 Fig.1 Flux ratio of convective and radiative heat transfer 时,由于钢卷和内罩温度都较低,辐射换热条件 during the annealing process 很差,炉内对流换热的热流密度可达辐射换热热 由式(1)可知炉内对流换热系数与保护气体 流密度的250倍以上,在均热阶段上述比值也保 的流速(循环流量)和热物性密切相关.图2给出 持在20倍以上.所以,对流换热在炉内换热中起 了氢气与钢卷四个表面的对流换热系数以及辐 着决定性的作用. 射等效换热系数随整个退火过程的变化规律, 全氢炉内保护气体与钢卷表面的对流换热 图中整个退火过程氢气的循环流量取为70×10 系数可由下式四计算: mh(在101.3kPa大气压下).由图可知,钢卷内外 -oms资台"2门 (1) 表面的对流换热系数最大,上下表面次之,等效 式中,"为气体流速,m/s;为气体导热系数,W/ 辐射换热系数最小,与对流换热系数相差1至2 (m·℃):,为气体动力粘度,Pas;p。为气体密度, 个数量级. kgm';D,为气体通流直径,m;L,为气体流动路线 氢气的密度仅为氮气的1/14,使得气氛循环 的长度,m;为由于漩流引起的传热强化系数. 风扇直径得以增大,气流循环量大大增加,加之 氢气的动力粘度仅是氨气的50%,因此,使用纯 收稿日期20020902林林男,35岁,讲师,博士研究生 氢气作为炉内保护气氛,无疑会大大提高炉内的 *国家自然科学基金资助项目No.50136020) 对流换热系数,见图3
第 卷 第 期 年 月 北 京 科 技 大 学 学 报 一 。 全氢罩式退火炉退火热过程的研究 — 对流换热系数和钢卷径 向等效导热系数的分析 林 林 ” 张 欣欣 ‘, 左 数 ‘, 向 顺 华 ” 刘 铁树 ” 钮龙 英 , 黄 夏 兰 ” 北京科技 大 学 机械工程学 院 , 北京 宝 钢研究 院设备研究所 , 上 海 摘 要 介绍 影 响全氢罩式退 火炉 内换 热 的两 个重 要 参数— 对 流换 热系数 和 钢 卷径 向等 效导 热系数 , 详尽 分析 了两个参数 的影 响 因素 , 并对 比 了氮气和 氢气气 氛下 两参数 的不 同 , 从而 在机理上 阐明 了全氢罩式炉 相 对传统 混氢罩 式炉 的优越性 , 为优化 炉 内换热提供 了理 论的依据 关 键词 全氢罩式退火炉 对流 换热 系数 等效导 热系数 分 类号 , 对 流 换 热 系 数 和 钢 卷径 向等效 导 热 系数 是 影 响全氢罩 式退 火 炉 以 下 简称全 氢炉 退 火 钢 卷传热 的两个最 主要参数 , 对这 两个参数 的研究 有 助 于 了解 炉 内换 热 的机理 , 为 全氢炉 内传 热 的 强 化 提供理论 的依据 对流 换热 系 数 全氢炉 内罩 内的换热 以对流 换 热 为 主‘ , 图 给 出 了 全 氢炉 对 流 换 热 的热 流 密 度 与 辐 射 换 热 的热 流 密度 之 比 由图可 知 , 在 整 个退 火过 程 开 始 和将 近 结束 时 , 由于钢 卷 和 内罩温度 都较低 , 辐 射换 热 条件 很差 , 炉 内对 流换 热 的热 流 密度 可 达 辐 射换 热热 流 密 度 的 倍 以 上 , 在 均 热 阶段 上 述 比值也 保 持在 倍 以上 所 以 , 对 流换 热 在 炉 内换 热 中起 着 决定 性 的作 用 全 氢 炉 内保 护 气体 与钢 卷 表 面 的 对 流 换 热 系 数可 由下 式 仪,计算 , , 一鳄二 。 ‘五、 一。 ,「,二扩夕 。 , 一 ’ 一 , 酬“ “ 户 ‘ ’ 式 中 , 气为气体流 速 , 耐 又为气体导 热 系数 , · ℃ 为 气体动 力 粘 度 , · 为气体 密度 , ‘ 为气体通 流 直径 , 为 气体 流 动路线 的长 度 , 古为 由于漩流 引起 的传 热 强 化 系 数 收稿 日期 刁 一 林林 男 , 岁 , 讲师 , 博 士 研究生 国家 自然 科学基金资助项 目困认 、弓乙 匕 一 一 一‘ 一一 一目一一一一一山一一 曰一一一 一山一 习 图 炉 内对 流换热 与 辐射换热热流密 度之 比 随退 火 过 程 的变化 卜 由式 可 知 炉 内对 流 换 热 系数 与保护气体 的流 速 循 环 流量 和 热物性 密切 相 关 图 给 出 了 氢 气 与 钢 卷 四 个表 面 的对 流换 热 系数 以及辐 射 等 效 换 热 系 数 随 整 个 退 火 过 程 的 变 化 规 律 , 图 中整 个 退 火 过 程 氢 气 的循环 流 量 取 为 恤 在 沙 大气压 下 由图可 知 , 钢卷 内外 表 面 的对 流 换热 系 数最 大 , 上 下 表 面 次 之 , 等效 辐 射换 热 系 数最小 , 与对 流 换 热 系数相 差 至 个数 量 级 氢 气 的密 度 仅为 氮气 的 , 使得气 氛循环 风扇 直 径 得 以 增 大 , 气 流 循环 量 大 大增 加 , 加 之 氢气 的动 力 粘度 仅 是 氮 气 的 , 因此 , 使 用 纯 氢气作 为炉 内保护 气 氛 , 无疑会大大提高炉 内的 对 流 换热 系 数 , 见 图 DOI :10.13374/j .issn1001—053x.2003.03.043
Vol.25 No.3 林林等:全氢罩式退火炉退火热过程的研究(山) ·255. 300 卷芯对流换热系数 250 200 卷边对流换热系数 M)/4 150 卷底对流换热系数 带 100 卷顶对流换热系数 50 等效辐射换热系数 0 0 10 2030 40 5060 t/h 图2退火过程中对流与等效辐射换热系数的变化规律 Fig.2 Convection heat transfer coefficient and the radi- ation heat transfer coefficient during annealing 350 300 —H21.900℃2.700℃ A一两层带钢间保护气体导热;B一两层带钢间辐射热 250 N23.500℃4.400℃ 交换;C一接触点导热;D一带钢内导热 5.100℃ 日200 图4钢卷径向传热模化单元体 M)/4 Fig.4 Schematic diagram of the modeled unit of a coil in the 150 radial direction 4 100 0 0 10152025303540 V/(m-s-) 图3氯气与氨气介质对流换热系数的比较 图5钢卷径向等效热阻网络模拟图 Fig.3 Convection heat transfer coefficient as a function of Fig.5 Schematic diagram of thermal resistances for the gas velocity for hydrogen and nitrogen at different tem- modeled unit of a coil in the radial direction peratures 经推导各个热阻的表达式为: 2钢卷径向等效导热系数 R2= 2λ 钢卷由上千层带钢组成,相邻两层带钢间气 2- 隙的存在使钢卷的径向导热率降低.因此,径向 R=4(1-0)eoTk b (2) 传热是制约全氢退火炉内钢卷热交换的主要 R。-0-pi 环节. 钢卷径向接触表面的换热形式参见图4叭.由 Ro-1.134tang 图可知,通过交界面的钢卷径向传热是通过相邻 式中,b和S分别为气隙与带钢厚度,m;s为带钢 两层带钢间气隙内传热的综合效应,包括通过保 导热系数,Wm·K):p为界面上点接触与表观面 护气体的导热、带钢间辐射换热及通过接触点的 积的比;e为带钢黑度;σ为斯蒂芬·波耳兹曼常 导热等.为简化计算,假设涉及的物性仅为单元 数, W(m2K);T。为相邻两层带钢的平均温度, 体平均温度的函数,且将圆柱面等效为平面,并 K;tan0为带钢表面形状的平均斜度;,为带钢表 按稳态传热的概念讨论热阻问题,于是得到如图 面粗糙度,m.研究表明,p和b的值可分别由下 5所示的钢卷内部导热热阻网络模拟图.图中Rs 式给出: 为12厚度带钢的导热热阻,R为两层带钢间的 P 9=P (3) 辐射热阻R。为保护气体的导热热阻,R为通过 b=42.7×10-exp(-5×10-2P) (4) 带钢接触点的导热热阻(实际上R,由通过表面 式中,H为接触固体中,材料较软者的硬度,MPa; 粗糙形成的点接触的导热阻和由热流流线变形 P为钢卷打卷张力,MPa. 形成的收缩热阻组成).以上所有热阻的单位均 单元体内总热阻R:可表示为带钢导热热阻 为m2.℃/W R与交接面接触热阻Re川的和:
、 林林等 全 氢罩 式退火炉退 火 热 过程 的研 究 · 卷芯对流换热系数 卷边对流换热系数 卷底对流换热系数 卷顶对流换热系数 等效辐射换热系数 乙,‘伟盛门以盈 ︵ 一园日 · ︾灵沙 坛 二二山一一一一一二一一一一司一 图 退 火过 程 中对流 与等效辐射换热 系数 的变化 规律 卜 介 “ 尸 一‘ 护 曰 , 二 一两层带钢间保护气体导热 一两层 带钢 间辐射热 交换 一接触点 导 热 一 带钢 内导 热 图 钢 卷径 向 传热模化 单元体 · 此 ︸,︶工”、,︸ ”八︷︸ · 甲日 · ︾去沐 川 · 一 ’ 圈 氮气与 氮气 介质对 流换热 系数 的 比较 ” ‘ 加 心 免 图 钢 卷径 向等效 热 阻 网络模拟 图 · 水 溉 日 ” 比 阳 经 推 导 各个 热阻 的表 达 式 为 ‘ 钢 卷 径 向等效导热 系 数 钢卷 由上 千 层 带 钢组 成 , 相 邻两 层 带钢 间气 隙的存在使钢卷 的径 向导 热率降低 因 此 , 径 向 传 热 是 制 约 全 氢 退 火 炉 内 钢 卷 热 交 换 的 主 要 环 节 钢卷径 向接触表 面 的换 热形 式参见 图 由 图可 知 , 通 过交界 面 的钢卷径 向传热 是通 过相 邻 两 层 带钢 间气 隙内传热 的综 合效应 ,包 括 通 过 保 护气体的导热 、 带钢 间辐 射换热及通 过 接触点 的 导 热 等 为 简化计算 , 假 设 涉及 的物 性 仅为 单元 体平均 温 度 的 函数 , 且 将 圆柱 面 等效为平 面 , 并 按稳 态传热 的概念讨论热阻 问题 , 于是得 到如 图 所示 的钢 卷 内部导 热热 阻 网络模拟 图 图 中 。 为 厚度 带 钢 的导 热 热阻 , 为 两层 带钢 间 的 辐 射 热 阻 。 为 保 护气 体 的 导 热 热 阻 , 。 为 通 过 带 钢 接触 点 的导 热 热 阻 实 际 上 。 由通过 表 面 粗糙 形 成 的 点接 触 的导 热阻 和 由热 流 流 线 变形 形 成 的收缩热 阻组 成 以上 所 有 热 阻 的单位 均 为 · ℃ 一心 一名 一﹃‘一刀毕 叮沙一几, 一了、 一 丁、了 跳’ 。 “ 外 又 势 一公 式 中 , 和 分别为气隙与带钢 厚 度 , 砚为带钢 导 热 系数 , 汉 切 为界 面上 点 接触 与表 观 面 积 的 比 。 为带钢黑度 , 为斯蒂芬 · 波耳 兹曼 常 数 , 珊 , · 几 为相 邻 两 层 带 钢 的平 均 温度 , 为带 钢 表 面 形 状 的平均 斜度 价 为带 钢表 面粗糙度 , 卿 研究表 明 , 价和 的值可 分别 由下 式 给 出’, 尸 尹节 二 二含万 ‘ 刀十尸 、 “ , 二 一‘ 一 一 式 中 , 为接触 固体 中 , 材料较软者 的硬度 , 尸为钢卷 打 卷张力 , 单元体 内 总热 阻 可 表示 为带 钢 导 热 热 阻 与交 接 面接 触 热阻 ‘叨 的和
·256· 北京科技大学学报 2003年第3期 R:=Rs+Rc (5) 图6(a)所示,当TRs;打卷 如前所述,钢卷径向的传热是带钢金属导热 张力大于13MPa时,R<R;显然钢卷径向导热在 及相邻两层带钢间气隙内传热的综合效应.由式 低打卷张力的情况下,主要制约因素是交界面的 (2)可知,单元体平均温度T和钢卷打卷张力P是 传热,且Re随P的变化比较显著. 影响径向热阻的主要因素.图6给出了带钢导热 交接面上的接触热阻是由R,R,R组成的, 热阻R和接触热阻R随T和P的变化规律.如 图7分别给出了以上三个热阻随T和P的变化规 4.5 4.5 4.0 -Rc 3 -Rs 4.0 Rc 3.5 1.4MPa 3.5 -Rs 2.8MPa 1.300℃ 3.0 3.12MPa 3.0 2.500℃ 2 3 3.700℃ 2.5 正 2.5 3 2.0 123 良 2.0 2 1.5 1.5 1.0 1.0 0100200300400500600700 2 46810 121416 T/℃ P/MPa (a)R和Rc随单元体平均温度的变化规律 (b)R,和R随钢卷打卷张力的变化规律 图6带钢导热热阻和气隙内并联热阻的变化规律 Fig.6 Variation of thermal resistances,Rs and Rc,with average temperature and compressive stress 0.25 2.0 一RR 1.4Ma 2.0 --Rp 2.8 MPa 1 R 1.300℃ 0.20 1.6 …RG 3.12 MPa 1.6 RD 2.500℃ 日 3.700℃ 0.15≥ 1.2 2 1.2全 123 0.10 0.8 3 2 2 3 0.05 0.4 0.4 123 0 ☐0.00 0.0 0.0 0 2 6 810 1214 16 0 100200300400500 600700 P/MPa T/℃ (a)R.,R,Rn随单元体平均温度的变化规律 (b)R,Ra,Rn随钢卷打卷张力的变化规律 图7缝隙内各个热阻的变化规律 Fig.7 Variation of thermal resistances,R,Re and Ro,with average temperature and compressive stress 律.如图7(a)所示,当T从0-700℃变化时,R比 0.60 Ra,Rn大3~4个数量级;R随温度升高而降低;R。 0.50 即随温度升高而升高,且Rc始终大于R.当P从 0.40 H: 1-15MPa变化时,R同样比R,R,大3-4个数量 0.30 级;R和R都随压力升高而降低,且R始终大于 0.20 R。.由以上分析可知,缝隙内的传热主要依靠保 N2 0.10 护气体的导热.因此,增加保护气体的导热系数 0.00 可以加大钢卷径向等效导热系数.由于氢气的导 0100200300400500600700 T/℃ 热系数大约是氮气的7倍,因此,使用100%氢气 图8氢气与氨气介质中钢卷径向等效导热系数与钢导 作为保护气氛,可大大提高钢卷径向导热率,加 热系数的比值 快钢卷的径向传热.氢气与氨气气氛下钢卷径向 Fig.8 Equivalent radial thermal conductivity of a coil for 等效导热系数的比较见图8. hydrogen and nitrogen atmosphere
一 北 京 科 技 大 学 学 报 年 第 期 勺石︸ 、尹、 ‘ 了、 其 中 , 如前所述 , 钢卷径 向的传热是带 钢 金属 导 热 及 相邻 两层 带 钢间气 隙 内传热 的综合效应 由式 可 知 , 单元体平 均 温度 和 钢 卷打卷 张力 尸是 影 响径 向热 阻 的主要 因素 图 给 出 了带 钢 导 热 热 阻 和 接触热 阻 。 随 和 尸 的变化 规律 如 图 所示 , 当 ℃ 时 , 。 , 交 界 面上 的传 热 是 制 约 钢卷 径 向传热 的 主要 因素 如 图 所 示 , 当钢卷打卷 张力 小 于 时 , 打卷 张 力 大于 时 , , 显 然钢卷径 向导热在 低 打卷 张力 的情况 下 , 主要 制 约 因素是交界 面 的 传 热 , 且 。 随 尸 的变 化 比较显 著 交 接 面上 的接触 热 阻是 由 凡 , , 凡组 成 的 , 图 分别给 出 了 以上 三 个热阻 随 和 尸 的变化规 ︵ 一沐 · 丫日 。 — 乙、叫咐 — , 一 一 一尹洲 ︵ 一一沐 · 岁,、叫。一︶日 仍 城 ℃ — , — ℃ 次 ℃ 迢 , ℃ 、 之耳卜、 一 、 、 乏泛遥 ’ 和 随单元体平 均 温度 的变化规律 尸 和 随 钢卷打卷 张力 的变化规律 图 带钢 导热 热 阻 和 气 隙 内并联热 阻 的变化 规 律 · 比 , , , ︵ ︸毕 · ︶岁、日叫‘ , ︶‘︹ 二 刁州川 一 幸工… 平不布蔫幕票票拿 —— 。 一 王 ︵沐 · ︶岁、日叫“ 乡 ’ · “ 窝 卜 乙 口 丫口 叫 省 丫军 ︵ 一气孟︸、 一毕 · 乙飞,、园 。 气 。 叫 尸 一一一一 一一一 一习一 一 ℃ 凡 , , 随单元体平均 温度 的变 化规律 凡 , , 。 随钢 卷打卷张 力 的变化规律 图 缝 隙 内各 个 热 阻 的 变化 规 律 , , 。 , 律 如 图 所示 , 当 了从 一 ℃ 变 化 时 , , 比 。 , 凡大 一 个数量 级 随温 度 升 高而 降低 凡 即 随温度 升 高而 升 高 , 且 始终 大 于 凡 当尸从 一 变 化 时 , 同样 比 。 , 尺。 大 一 个 数量 级 。 和 都 随压力 升 高而 降低 , 且 。 始终大 于 由 以上 分析 可 知 , 缝 隙 内的传热 主要 依 靠保 护 气体 的导 热 因 此 , 增 加 保护 气体的 导 热 系数 可 以加 大 钢卷径 向等效 导 热系 数 由于 氢气 的导 热 系数大约是氮 气 的 倍 , 因 此 , 使用 氢气 作 为保护 气 氛 , 可 大 大提 高钢 卷径 向导 热 率 , 加 快钢 卷 的径 向传热 氢气 与氮气气 氛 下 钢 卷径 向 等效导 热 系数 的 比较见 图 省 ︸内、 ﹄ 叹 ℃ 图 氮气与 氮气 介质 中钢 卷径 向等效导热 系数 与钢 导 热 系 数 的 比值
Vol.25 No.3 林林等:全氢罩式退火炉退火热过程的研究(山) ·257· 将式(2(6)带入热阻的定义式,可得钢卷径 参考文献 向等效导热系数的计算公式(⑦),其中各参数的取 1 Zuo Yi,Wu Wenfei,Zhang Xinxin,et al.A study of heat 值如下:S=0.806mm;H=1133.86MPa;P=8MPa; transfer in high-performance hydrogen bell-type annea- T.=573K;=42.0W/mK):入=0.307WmK);e ling furnaces[J].Heat Transfer-Asian Research,2001.3 =0.2;om=3.2215m;tan0=0.08554;g=5.67×10 (8:615 =W/(m2.K'). 2 PerrinAR.GuthrieRIL,Stonehill BC.The process tech- 元=4b= nology of batch annealing [J].Iron and Steel Making, Rs 1988(10):27 14s2-6 b 3 Wolfgang Potke,Jeschar Rudof,Kehse Georg.Warmetra- nsport beim gluhen von coils in haubenofen mit,wasser- 3结论 stoff als schutygas []Stahl u Eisen,1988,108(12):581 4 Pullen J,Williamson JB P.On the plastic contact of rough (1)分析了全氢罩式退火炉内对流换热系数 surface [A].Proc R Soc [M].London,1972.159 和钢卷径向等效导热系数的影响因素,对比了氢 5 Park Seong-Jun,Hong Byung-Hee,Baik Seung-Chul,et al. 气和氨气气氛下以上两参数的不同,从机理上阐 Finite element analysis of hot rolled coil cooling [J],ISIJ 明了全氢气氛相对于传统混氢气氛的优越性. International,1998,38(11):1262 6殷跷静,Degievanni A.同心园柱套筒间接触热阻研 (2)使用网络模拟的方法给出了钢卷径向等 究(0:数学模型U.北京科技大学学报,1996,184: 效导热系数的计算公式,并分析了单元体平均温 383 度和钢卷打卷张力对径向热阻的影响,从而为研 7殷晓静,DegievanniA.同心园柱套简间接触热阻研 究钢卷径向传热机理和优化炉内换热提供了科 究():参数影响分析北京科技大学学报,1996,18 学的依据. (4):464 Study on Annealing Thermal Process in HPH Furnace(II): Analysis of Convection Heat Transfer Coefficient and Equivalent Radial Thermal Conductivity LIN Lin",ZHANG Xinxin",ZUO Y,XIANG Shunhua,LIU Tieshu,NIU Longying,HUANG Xialan 1)Mechanical Enginecring School,University of Science and Technology Beijing,Beijing 100083,China 2)Equipment Research Dept,Shanghai Baosteel Research Institute,Shanghai 201900,China ABSTRACT Convection heat transfer coefficient and equivalent radial thermal conductivity are discussed,which have great effect on heat transfer in a HPH furnace.Using the approach of mathematical modeling,the convection heat transfer coefficient on each surface of a coil are studied for the complete cyele of annealing.In view of the analysis of thermal resistances,the equivalent radial thermal conductivity is investigated for different temperatures, compressive stresses,and atmosphere gases,respectively.It is shown that the HPH technique can significantly en- hance the heat transfer compared with the traditional HHX method. KEY WORDS HPH furnace;convection heat transfer coefficient;equivalent radial thermal conductivity
心 林林 等 全 氢 罩 式退 火 炉 退 火 热 过 程 的研 究 · · 将 式 卜 带 人 热 阻 的定 义 式 , 可 得 钢 卷 径 向等效 导热 系数 的计算公式 , 其 中各参数 的取 值 如下 尸 二 又, · 又 , · 。 拌 “ 一 , · 又 竺泛二£ 令坐龚畔 十 卜, 含 十少鬓呵 卯 仓 ’ ‘ 一 ’ 一 ’ , 结 论 分析 了全 氢罩 式 退 火 炉 内对 流换 热 系数 和 钢卷径 向等效导 热 系数 的影 响 因素 , 对 比 了氢 气 和 氮气 气氛 下 以 上 两参数 的不 同 , 从 机理 上 阐 明 了全 氢 气氛相 对 于传统 混 氢气 氛 的优越 性 使用 网络模 拟 的方 法 给 出 了钢 卷径 向等 效导 热 系数的计算公式 , 并 分析 了单元体平均 温 度 和钢卷 打卷 张力对径 向热 阻 的影 响 , 从而 为研 究 钢 卷 径 向传热 机 理 和 优 化 炉 内换 热 提 供 了科 学 的依据 参 考 文 献 砚 , 认 , , 一 一 【 」 一 , , , , , 场 邵 , , “ , , , , 【 , 一 , 一 , 一 , 』 , , , 殷 晓静 , 同心 园柱套筒 间接触热 阻研 究 数学模型 【 匕京科技大学学报 , , 殷 晓静 , 同心 园柱 套筒 间接触热 阻研 究 参数影 响分析 」北 京科技 大学学 报 , , ” , ” , 万 ,’ , 尤翻 刃 ,气 , 心 ,, 介 ,, , , , , , , , · , , , ,