墙体传热反应系数的计算.PDF

第28卷第1期 东 南 大 学 学 报 Vol28 No1 1998 年1月 JOURNAL OF SOUTHEAST UNIVERSITY Jan. 1998 何嘉鹏 徐明立 ( 南京建筑工程学院城建系,南京210009) (东南大学建筑设计研究院, 南京210096) 王东方 (南京建筑工程学院机 系,南京210009) 研究空调房间围护结构的传热动态特性, 寻找空调房间墙体结构传热传递函数的求解方 法, 是进行动态 荷计算和空调系统运行最优控制所必须解决的一个重要问题.求解墙体非稳 态传热的过程, 就是要求解它的温度场和热流场随时间的变化规律.60 年代末、70 年代初,加 拿大的D. .Stephenson 和 .P.Mitalas提出的反应系数法就是利用了拉普拉斯变换和反变 换, 并将现代控制论引入其中, 将墙体作为一个热力系统, 将任意随时间变化的温度边界条件 分解成一系列的三角脉冲, 预先求出热力系统反应系数,认为各脉冲线性无关, 热力系统本身 是线性系统.应用叠加原理,将每一脉冲引起的热流变化按卷积原则进行叠加,给出系统的热 反应. 1 用传热反应系数法计算墙体的非稳态传热, 是把系统的热特性与外部干扰分开, 认为反应 系数是墙体本身的热特性, 而与外部干扰无关. 只要知道系统的反应系数以及干扰的瞬时值和 历史值, 就可以求出该时刻的热流. 1. 1 求解墙体内部温度场分布和热流随时间的变化规律, 即求解式( 1) 所示偏微分方程组 t ( x , ) t 2( x , ) = a 2 , 0 x l , 0 x ( 1) t ( x , ) 10 , 0 x l , 0 q ( x , ) = - x 文献[ 1] 利用拉氏变换法对式( 1) 进行变换, 墙体的温度分布和热流为 T ( x , s ) = ch( s/ a x ) T ( 0, s ) - sh( s/ a x ) Q ( 0, s) / ( s/ a ) ( 2) Q ( x , s ) = - s/ a sh( s/ a x ) T ( 0, s) + ch( s/ a x ) Q ( 0, s) 收稿 日期:1997- 04- 02, 修改稿收到 日期: 1997- 11- 04. 第1期 何嘉鹏等:墙体传热反应系数的计算 137 式中, T ( x , s) 为墙体 x 处温度对时间的拉氏变换; Q ( x , s) 为 x 处热流对时间的拉氏变换; s 为 对应时间 的拉氏变换象时间的参变量; a 为热扩散系数; x 为计算面坐标值; l 为墙体厚度. 由式( 1) 和式( 2) 可以看出, 只要知道边界条件 t ( 0, ) 和 q( 0, ) , 就可以通过拉氏变换求 得墙体任一部位的温度响应和热流响应, 再对 T ( x , s) 和 Q ( x , s) 进行拉氏反变换, 就能得到 最终解. 1. 2 式( 2) 中, 令 x = l , 并且用矩阵的形式表达可得 T ( l, s