PDF《建筑环境设计模拟分析软件 DeST》

建筑环境设计模拟分析软件 DeST

第二章、建筑动态热过程模型 清华大学 谢晓娜 宋芳婷 燕达江亿谢晓娜,女,1979 年8月生,大学,在读博士研究生

100084 北京清华大学建筑学院建筑技术科学系 (010)62789761 xxna01@mails.

tsinghua.edu.cn 摘要 本文介绍了 DeST 中基于状态空间法建立的建筑动态热过程模型.建筑热过程模型是一 切建筑热环境及建筑环境控制系统模拟分析的基础, 因此本章首先介绍了建立建筑热过程模 型中的主要问题和基本方程, 指出建立建筑热动态模型必须求解壁体传热的偏微分方程、 必 须同时考虑构成房间的各个部件的影响以及房间之间的相互影响. 以此为前提, 详细介绍了 DeST 对这些问题的处理方法和具体算法,并以实例阐明 DeST 对建筑热环境的模拟分析.为 了突出状态空间法的特点, 本章最后从算法和建立动态热过程模型的方式上对该方法与其他 几种方法进行了比较. 关键字 状态空间法 动态热模型 建筑热平衡 Abstract The dynamic building thermal process model of DeST based on State Space solution method is explained in this paper. The dynamic thermal model is the foundation to describe the physics of building thermal environment and its control systems, the problems encountered and the equations to be solved when constructing such a model is, therefore, presented first. It is necessary to solve the heat transfer equations of opaque envelopes correctly and to consider the interaction among components in a room and rooms in a building in details. The detail procedures on how to solve these problems and how to carry out a building thermal environment simulation by DeST are explained with examples. Finally, a comparison between State Space methods and other solution methods is made. Key Word: State-Space method, dynamic thermal process model, building heat balance

一、建筑热过程模拟的主要问题 建筑热过程模拟是一切建筑热环境及建筑环境控制系统模拟分析的基础. 它的基本问题 就是给出在不同的气象条件下、不同的使用状况以及不同的环境控制系统(采暖空调系统) 作用下, 建筑物内温度的变化情况. 在此基础上就可以分析预测不同的建筑设计与围护结构 形式可能形成的室内温度状况, 了解为了维持要求的建筑物热环境所需要消耗的冷热量, 以 及不同形式的采暖空调系统作用在这一建筑中的运行品质. 要解决建筑热过程模拟的基本问题, 必须根据实际的建筑热过程建立准确的物理模型和 数学模型. 首先, 应全面准确地归纳出影响建筑热过程的各种因素以及建筑热过程满足的物 理规律,并在此基础上建立相应的物理模型;

然后,在此基础上建立准确表述建筑热过程的 数学方程.

1、建筑热过程模拟的物理模型 物理模型的示意如图

1 所示: 图1建筑热过程的物理模型 图1中各符号含义如表

1 所示: 表1建筑热过程中的各种热扰 o t 室外空气温度 out h 建筑围护外表面与室外空气的对流换热 系数 adjacent t 邻室空气温度 in h 建筑围护内表面与房间内空气的对流换 热系数 env t 建筑物周围环境的平均温度 , r o h 建筑围护外表面与周围环境的辐射换热 系数 ground t 地下土壤的温度 _ _ sun radiation wall Q 太阳辐射中被建筑不透明围护外表面吸 收的部分 _ _ sun radiation window Q 太阳辐射中被窗户吸收的部分 _ sun trans Q 太阳辐射中透过窗户进入室内的部分 ocupants Q 室内的人员扰量 lights Q 室内的灯光扰量 equipments Q 室内的设备扰量 longwave Q 各围护内表面之间的长波辐射换热 r h 各围护内表面之间的辐射换热 系数 o G 房间与室外的通风换气量 adjacent G 房间与邻室的通风换气量 hvac Q 房间空调设备的供热(冷)量_adjacent radiation Q 墙体靠近邻室侧接收到的邻室 的太阳透窗辐射和室内内扰的 辐射部分 对于建筑中的一个房间而言, 其热过程主要包括四个方面: 外扰通过围护结构的热传递 过程,内扰的热传递过程,室内外通风和空调投入热量. 一个孤立房间的外扰包括室外气象条件 (如外温、 辐射、 风向风速等) 和环境热状况 (如 周围环境表面的温度) .而对于一栋建筑中的某一个房间,其邻室的热状况也可看作一种外 扰.这些扰量主要是通过围护结构热传递影响房间内表面温度,并通过对流、辐射影响室内 空气.但是,太阳辐射还可透过窗户直接影响房间内表面温度. 内扰对房间的热作用,包括潜热和显热两方面.人体和设备的散湿伴随着潜热散热,它 们直接作用到室内空气,立刻影响室内空气的焓值.而照明、人体和设备的显热散热则以两 种方式在室内进行热交换, 一种是以对流方式直接传给室内空气, 另一种则以辐射形式向周 围各表面传递,之后再通过各表面和室内空气之间的对流换热,逐渐传递给室内空气. 室内外通风直接与室内空气混合,因此它们的热量和湿量会立刻影响室内空气的热状 态,改变室内空气的温度和湿度[1] . 空调热量可能以直接送风的方式进入房间, 这种情况与室内外通风类似. 但随着各种新 的空调方式(如辐射板)的出现,空调热量影响室内热环境的方式也不再单一.

2、建筑热过程模拟的数学模型 大多数建筑都由大量的具有很大热惯性的钢筋混凝土或砖砌体组成.当外温或室内温度 变化时,这些砌体内的温度也会相对缓慢变化,同时吸收或放出热量.因此,必须准确描述 这一蓄热放热的动态过程,才能准确模拟和预测建筑热过程. 对于一面由多层材料组成的墙体,通常其沿表面方向结构均匀,且厚度远小于表面长宽 尺度, 因此可以忽略其内部沿平行于表面方向的导热, 按一维过程分析其沿厚度方向的导热. ( ) P t t t C k x x ρ τ ? ? ? = ? ? ? (1) 式中,t 为壁体内的温度;

P C 、 ρ 分别为壁体材料的比热(kJ/kg.K)和密度(kg/m3 ),它们 随沿厚度方向的不同材料层取不同的值;

k 为沿厚度方向的导热系数(W/m・K),其值亦随不 同材料层而异;

x 为厚度方向. 在室内一侧,式(1)的边界条件为: , ( ) ( ) x l in a r j j r in j t k h t t q hr t t q x = ? ? ∑ (2) 式中 in h 为壁体内表面与空气的对流换热系数(W/m2 .K);

a t 为房间室温;

r q 为壁体内表 面吸收的透过窗户的太阳辐射热量(W/m2 );

, r in q 为室内其它热源(人员、灯光、设备等)以 辐射方式传至该表面的热量;

j hr 为温度为 j t 的另一表面 j 与该表面的长波辐射换热系数, 该系数的量纲与对流换热系数相同. 将长波辐射换热表示成与温差的线性关系从而简化动态 模拟方程求解, 这种方法是热环境动态模拟中处理长波辐射问题的常用方法, 辐射换热系数 j hr 的确定将在本章第二部分进行讨论.求和号∑指对所有可见表面求和. 当壁体另一侧也是室内时(即壁体为内隔断或楼板) ,其边界条件同式(2),只是变为

0 x = ,且k前面无负号.当壁体另一侧为室外时,则,

0 , ( ) ( ) x out o r o out env t k h t t q hr t t x = ? ? (3) 式中 out h 为壁体外表面与室外空气的对流换热系数(W/m2 .K);

o t 为室外空气温度;

, r o q 为壁体外表面吸收的太阳辐射热量;

out hr 为壁体外表面与周围环境表面的长波辐射换热系 数,这里周围环境表面主要是指其他建筑表面、地面和天空,由于各个环境表面的温度不同 [2] ,需要根据壁体外表面对环境表面的角系数确定长波辐射换热计算的几何关系,但是实际 建筑的周围环境表面情况复杂,要全面准确地描述这种换热关系十分困难,因此,将周围环 境表面作为一个表面看待;

env t 代表该表面的综合温度,这样壁体外表面与该表面的角系数 总是等于 1;

out hr 的确定则与壁体内表面辐射换热系数 j hr 的确定方法类似,至于 env t 的取 值,将在本连载的

第五章进行讨论. 对建筑物内所有壁体(楼板、墙体、门窗)均可列出式(1)到(3)的动态传热方程.这些 围护结构在建筑物内围合成许多建筑空间, 对每个建筑空间内的空气温度变化可列出如下方 程: cov

1 n a Pa a a j in j a f vent hvac j dt C V F h t t q q q q d ρ τ τ τ = ∑ (4) 式中 Pa a a C V ρ 为建筑空间内空气的热容(KJ/℃), j F 为建筑空间的壁体内表面 j 的面积, j t 为j表面的温度,n 为内表面个数, cov q 为室内热源(人员、灯光、设备等)以对流方式 传递给空气的热量[3] , f q 为室内家具放出的热量, vent q 为由于室内外空气交换或与邻室的 空气交换带入室内的热量, hvac q 为供热空调系统送入建筑空间的冷热量. 由于壁体的蓄热放热作用不能忽视,导致建筑热过程的模拟必须求解式(1)的偏微分方 程;

而式(2)中()jjjhr t t ? ∑ 一项描述的各壁体内表面之间的长波辐射换热以及式(4)中给出 的室内空气温度与各表面之间的热耦合关系, 导致建筑热过程的分析求解必须同时考虑整个 房间的各个部分,而不能简化地仅基于对外墙及屋顶的分析计算;

另一方面,当壁体的另一 侧也是室内时,壁体外表面与邻室空气温度的热耦合关系以及式(4)中vent q 包含的与邻室的 空气交换, 因而建筑热过程的分析求解必须同时考虑建筑中的所有房间. 这些问题构成了建 筑热环境动态模拟的主要难点, 因而也是各种模拟方法的主要不同所在. 本章着重介绍 DeST 中建筑热过程的求解方法,尤其是对上述问题的处理方法与具体算法. 此外,建筑热环境的动态模拟还将涉及如下问题:

9 外窗的遮阳[4] ,太阳光的透过[5] ,窗的导热,及式(2)中的 r q ,即透过窗户的太阳辐 射被各内表面吸收的计算,这将在本连载的

第四章做详细讨论;

9 室内家具的处理,这将在本章第二部分进行讨论.

9 对于紧接地面的建筑空间和半地下建筑空间而言,地面以下的围护(如地板、地下 的........