01-035基于冷负荷系数法的办公建筑高大中庭的空调负荷简化算法文章

基于冷负荷系数法的办公建筑高大中庭的 空调负荷简化算法 北京市建筑设计研究院 李丹 徐宏庆 马晓钧 摘要 本文以一北京地区实际工程为例 介绍了具有高大中庭的办公建筑的空调冷负荷计算简化方法 对其中庭和与中庭相邻的办公室进行了较详细的算法介绍 本文对同类建筑的空调负计算具有参考借 鉴作用 关键词 办公建筑 中庭 冷负荷 冷负荷系数法 1 1 引言引言 随着越来越多的大型办公建筑的出现 为保持其具有良好的使用性 中庭也随之在设计中得以广 泛的应用 虽然中庭也属高大空间 但办公建筑的中庭在其热环境和空调负荷的计算上 有其特殊性 对于一般高大空间而言 往往只关注下部人活动区域的热环境状况和负荷计算 而大型办公建筑内部 的中庭通常与相当数量的办公房间相邻 因而其整个空间的热环境和垂直温度分布 即会对与其相邻 的办公室的空调负荷产生较大的影响 所以 在中庭及其周边办公室空调负荷的计算上 则需考虑它 们之间的热量转移与交换 目前在国内工程设计领域 计算冷负荷最常用的方法是冷负荷系数法 然而冷负荷系数法基于集 总参数概念 无法准确描述大空间的温度及负荷的分布状况 因而在计算大空间建筑的空调负荷是具 有一定的局限性 但是 冷负荷系数法具有简便 快速 物理意义直观的特点 在我国的工程设计中应用广泛 为 广大工程设计人员所熟知和掌握 一直以来 国内外学者和工程设计人员也不断对冷负荷系数法进行 了修正和拓展 使其适用性得以更大地增强 本文以北京地区一实际工程为例 基于冷负荷系数法原理 从热传递和负荷产生的物理过程入手 对办公建筑高大中庭的空调负荷简化算法进行探讨 2 2 中庭热传递过程的简化中庭热传递过程的简化 2 2 1 1 简化计算原则简化计算原则 太阳辐射光照射到玻璃上 入射部分的太阳能 e 由三部分组成 透射部分 反射部分 吸收部 分 三者之间的关系如下 见图 1 e e e 1 1 式中 e 太阳光直射透射比 e 太阳光直接反射比 e 太阳光直接吸收比 玻璃吸收部分的能量又以热对流 辐射方式向室外 室内热传递 二者关系如下 e q0 qi 2 式中 q0 室外侧热传递百分比 qi 室内侧热传递百分比 则太阳能总透射比 图 1 太阳能通过玻璃的物理过程示意图 图 2 进入中庭的太阳辐射分析 g e qi 3 式中 g 太阳能总透射比 e 太阳能直射透射比 qi 太阳能二次热传递系数 对双层玻璃而言 2 12 111 e ee e i ie nG q nnG 式中 e1 室外侧玻璃的太阳直接吸收比 e2 室内侧玻璃的太阳直接吸收比 G 两片玻璃间的热阻 W m2 k ni 室内侧玻璃表面的对流放热系数 W m2 k ne 室外侧玻璃表面的对流放热系数 W m2 k 根据上述对太阳能通过玻璃的传热机理的分析 对于中庭外窗得热进行简化计算 因此 简化计 算主要考虑窗玻璃的影响因素 由 于外窗玻璃所占面积很大 故玻璃 对太阳能的光热特性 即为整个外 窗的光热特性 考虑到计算的简便 性 由室内向外照射的太阳能 只 计算经玻璃反射和吸收后留在室内 的一次热量 二次及更高次反射的 能量不计 考虑到工程计算量以及 工程精度的要求 进行简化处理也 是可行的 2 2 2 2 简化计算主要过程简化计算主要过程 首先对进入到中庭的热量进行分析 中庭的热平衡方程为式 2 2 outin QQQ 0 其中 留在中庭的热量 此部分将做为中庭的得热进行计算 0 Q 进入中庭的热量 此部分包括 进入中庭的总太阳辐射热量 中庭内扰产生的热量 人 in Q 员灯光等 从中庭流出的热量 此部分包括 通过中庭外围护结构流出室外的热量 通过中庭进入 out Q 办公室的太阳辐射热 进入办公室的热量除辐射热外 还包括通过壁面的传导热 但考虑到中庭与办公室空气温差在计 算中庭温度梯度之前不可知 故不能做为已知热源条件计算中庭温度场 而且这部分热量随高度的增 加而增加 对中庭下层人活动区域的负荷计算影响几乎没有影响 另外 考虑到工程实际情况 这部 分热量在计入中庭负荷中是偏安全的 因此通过中庭进入办公室的传导热做为留在中庭的热量考虑 在计算中需考虑太阳入射位置与中庭得热关系 因此为计算中判断太阳辐射源位置和方向的简便 本文约定 太阳位置角 a 为太阳与水平面的夹角在正南方向的投影角 基于前述简化原则 将太阳 位置角 a 作为太阳能在不同时刻计算室内得热的依据 根据不同时刻太阳位置角 a 的变化 计算出投 射到中庭内各个壁面的热量按面积所占太阳入射面面积的比例 从而确定中庭各个壁面的太阳辐射负 荷值 对于室内各个壁面 将室内地面和墙面考虑成为具有相同热工特性的蓄热体 用以简化地面和墙 面阴影面积的复杂计算 同样 出于计算安全考虑 对建筑物自身遮挡的热量将不做折减 对散射光 部分 按有关设计手册 1 考虑 此处不再进行修正 在确定中庭各个时刻的得热量之后 进行中庭垂直温度分布的计算 办公建筑内部的中庭通常与 相当数量的办公房间相邻 其整个空间的热环境和垂直温度分布 即会对与其相邻的办公室的空调负 荷产生较大的影响 因而 了解中庭垂直温度分布 将是计算与中庭相邻的办公室空调负荷的前提条 件 最后 将计算得出垂直温度分布的情况 做为计算办公房间空调负荷的外部计算条件 最终计算 得出与中庭相邻的办公房间的空调负荷 3 3 中庭冷负荷计算实例中庭冷负荷计算实例 3 3 1 1 中庭冷负荷计算中庭冷负荷计算 以下将以某工程为例 介绍办公建筑高大中庭的空调负荷简化算法的计算过程 该工程中庭分为 主中庭 南庭和北庭三部分 其中主中庭面积 1700 平米 高度 52 米 北庭面积 880 平米 高度 33 米 南庭面积 880 平米 高度 33 米 根据北京地区气象资料 得出不同时刻的太阳位置角 如表 1 所示 表 1 北京地区夏季计算日的不同时刻的太阳位置角 时刻 6 0 0 7 0 0 8 0 0 9 0 0 10 0011 0012 0013 0014 0015 0016 0017 0018 00 太阳 位置 角 13 524 23647 961 47590105118 6132 1144155 8166 5 依据上文对中庭得热计算中的简化原则和简化条件 分别得出该办公建筑三个中庭的得热量 进 而依据冷负荷系数法可计算出其夏季逐时空调冷负荷 其夏季计算日不同时刻中庭空调负荷的分项统 计情况见表 2 表 4 表 2 主中庭不同时刻空调负荷分项统计表 单位 kw 时刻 6 007 008 009 0010 0011 0012 0013 0014 0015 0016 0017 0018 00 顶面 63 7 91 6 123 5 155 4 210 6 264 3 312 2 344 1 304 7 188 7 168 1 167 3 135 4 东面 112 3166 6208 8280 1263 6172 5139 0134 2129 4119 8115 0100 791 06 西面 57 5162 3167 1071 8976 6881 4881 48115 0172 5214 7221 6245 9209 6 地面 26 0348 7770 8683 93110 0135 9164 0172 9151 591 1186 50126 7104 2 人 总 259 5 369 3 470 3 591 3 660 9 654 2 696 7 766 2 758 1 614 3 591 2 640 6 540 3 表 3 北庭不同时刻空调负荷分项统计表 单位 kw 时刻 6 007 008 009 0010 0011 0012 0013 0014 0015 0016 0017 0018 00 顶面 39 0 56 0 75 5 95 0 114 5 137 0 160 8 194 7 191 9 155 2 138 3 118 5 95 9 西面 27 07 29 32 31 58 33 83 36 09 38 35 38 35 54 13 81 20 103 8 117 3 139 9 126 3 地面 人 9 40 24 45 37 75 42 28 43 01 46 75 57 87 68 13 72 68 68 61 66 60 93 39 77 44 总 75 5 109 8 144 8 171 1 193 6 222 1 257 0 317 0 345 8 327 6 322 2 351 8 299 6 表 4 南庭不同时刻空调负荷分项统计表 单位 kw 时刻 6 007 008 009 0010 0011 0012 0013 0014 0015 0016 0017 0018 00 顶面 42 9 61 7 83 1 104 6 126 0 142 1 152 8 185 0 182 3 147 5 131 4 112 6 91 2 地面 人 18 36 35 12 38 39 27 34 6 10 29 17 102 0 55 91 52 12 54 80 64 32 71 93 31 43 总 84 6 143 8 173 0 167 3 137 4 209 0 395 1 316 2 304 3 275 9 282 7 282 2 163 3 3 3 2 2 与中庭相邻的办公室负荷计算与中庭相邻的办公室负荷计算 前文提到 在确定中庭各项得热以后 以次作为计算条件 计算中庭垂直温度分布 本工程采用 CFD 模拟方法对中庭进行数值计算 得出垂直温度分布的情况 从而做为计算办公房间空调负荷的外 部计算条件 本次模拟计算采用清华大学建筑技术与科学系开发的通用三维流动与传热的数值模拟程 序 STACH 3 湍流模型为专门针对建筑环境工程提出的 MIT 零方程模型 其模拟的准确性已在有关文 献中 2 论及 在此不再赘述 计算网格采用正交的结构化网格 数值算法中的离散方法为有限容积法 动量方程采用交错网格 即 U V W 方程 差分格式选择幂指数格式 求解算法为 SIMPLE 算法 代数方程求解采用交替方 向迭代解法 ADI 法 边界条件的处理采用把其处理为各相应变量的源项的方法进行 中庭物理模型见图 3 图 5 各个壁面的热源情况按表 2 表 4 选取 开口边界条件依照实际工程 设计方案确定 由表 2 表 4 可知 主中庭计算日空调负荷的最大值出现在下午 14 00 时刻 北庭负荷的最大值 出现在下午 17 00 时刻 南庭负荷的最大值出现在中午 12 00 时刻 故选择主中庭 14 00 时刻 北 庭 17 00 时刻 南庭 12 00 时刻做为模拟计算时刻 从而得出最不利时刻的各中庭垂直温度分布 图 3 主中庭数值计算物理模型 图 4 北庭数值计算物理模型 通过数值计算 得出主中庭在模拟计算时刻下 空间 内的温度场和空气速度场的分布情况 如图 6 图 7 所示 可以看出 随着高度的增加 中庭温度增加较高 对于顶 层办公室来说 其室外相邻中庭的空气温度达到 48 6 远高于夏季室外计算温度值 因而 采用模拟计算得出的 空气温度做为办公室负荷的计算条件 将比用估算办公室 室外温度的方法更加合理 从主中庭垂直温度分布计算结 果可以看 图 5 南庭数值计算物理模型 出 沿高度的温度梯度变化 为随高度的增加呈加大的趋 势 图 6 主中庭温度垂直分布图 图 7 主中庭风速垂直分布图 与主中庭类似 分别计算出北庭和南庭在模拟计算时刻下 空间内的温度场和空气速度场的分布 情况 图 8 图 9 为北庭温度场和速度场的垂直分布情况 图 10 图 11 为南庭温度场和速度场的垂直 分布情况 从北庭 南庭垂直温度分布计算结果可以看出 沿高度的温度梯度趋势与主中庭一致 图 8 北庭温度垂直分布图 图 9 北庭风速垂直分布图 经过对中庭温度场的数值模拟 得出中庭空气在不同高度的温度值 用于计算相对应楼层的办公 房间的围护结构负荷 由于办公房间数量较多 其空调负荷计算结果在此不再赘述 图 10 南庭温度垂直分布图 图 11 南庭风速垂直分布图 4 4 结论结论 通过所介绍的工程实例计算 本文对具有高大中庭的大型办公建筑的空调负荷计算进行了一定的 探索 采用了基于冷负荷系数法与 CFD 数值模拟相结合的计算方法 应用于中庭以及与其相邻的办公 房间的空调负荷计算 以往在计算此类建筑时 通常仅采用冷负荷系数法计算出中庭总得热量 然后依负荷转移系数 得出中庭人活动区域的空调负荷 然后再依沿高度的空气温度变化系数 获得中庭垂直温度分布 本 文采用的方法与原有方法相比较 热量