建筑环境设计模拟工具包DeST

1、建筑环境设计模拟工具包DeST             摘要建筑环境设计模拟工具包DeST是基于功能的模拟软件，用于对建筑、方案、系统及水力计算进行模拟，以校核设计，保证设计的可靠性。介绍了DeST的结构、用户界面，并结合工程实例说明了DeST的应用。 关键词建筑环境模拟 模拟软件 DeST 1 开发DeST的目的近30年来，建筑和空调系统的模拟被广泛地应用在学术研究领域，出现了很多的模拟模型、模拟方法以及模拟应用工具。模拟技术已经相对成熟，但在实际设计中，采用模拟技术辅助分

2、析的仍然很少，大多数设计仍然仅考虑最不利工况，而没有考虑全年的运行过程，这导致了诸如设备选择不合理、过渡季无法满足要求等问题。最近，国内外的研究机构和设计公司开始投入越来越多的力量以缩小模拟和设计之间的应用鸿沟。例如在国际能源组织（IEA）最近结束的研究子项ANNEX 301 "Bring simulation into application"中，设计过程中的模拟技术应用研究是其最主要的章节。虽然，采用模拟分析的手段可以提高设计的可靠性，但只有在明确了设计和模拟之间的关系之后，才能制造出在实际设计过程中能够被有效利用的模拟分析工具。而现有国外几种主流的模拟分析工具由于在开

3、发时没有充分考虑设计过程的阶段性、延续性等特点，只能用于学术研究领域。 设计过程是一个阶段化的过程，包括初步设计、方案设计、详细设计以及后设计阶段。这是一个复杂的、不断反馈的过程，目前是否存在一个模拟分析工具能够服务于整个过程呢？ 目前的模拟工具可以划分成两大类：一类是基于功能的，以DOE为主要代表；另一类是基于模声的，以TRNSYS，HVACSIM+,MATLAB为代表。基于功能的模拟工具从满足某种功能要求（如计算建筑全年能耗）出发来设计模拟系统；基于模块的模拟工具注重于构建系统的灵活性，其特点是采用通用的模块接口和统一的非线性求解核心。 当所有的条件都已知时，使用模块化的模拟系统可以很方便

4、地建立起整个系统的框架并进行模拟计算，这使得它们非常适用于研究领域，通过模拟计算去了解在系统的各个部件确定后系统的运行状态。但在设计过程中，尤其是在设计的初步阶段，设计人员无法掌握所有的信息，某些数据是模拟的输出而不是输入。例如，在方案设计阶段，当设计人员试图比较不同的系统形式时，因为空调机组的选择应该在方案确定之后进行，他无法了解空调机组的具体信息。为了使模块化的模拟程序能够运行，使用者不得不采用"缺省"的部件，选择某种"缺省的空调机组"来构建模拟系统。计算机模拟软件一个很重要的特性是输出严重地依赖于输入。因此，模拟由诸多"缺省"的

5、部件组成的系统，其结果对实际的设计不具备指导意义，当其"缺省"的是后续设计的目标时，这样的模拟计算无助于设计人员进行决策。因此，模块化的模拟工具虽然适用于学术研究，但并不适用于工程设计。 与模块工具相比，基于功能的模拟工具灵活性较低，但是更接近于设计人员思路，因此较容易被采用。为了把握全年的运行特性，设计人员通常用其来计算建筑物全年的能耗要求。作为其代表的DOE由建筑模拟、系统模拟和机组模拟三大部分组成。但其中各模块之间的关联存在着缺陷。例如，"在空调系统的模拟中，假设送风温度是已知的。这样的处理对于简单的运行方式（例如已经确定每一个时刻的送风温度）是有效的，但对

6、于较复杂的系统则无法工作，需要采用前一时刻的数据，或者建筑物的模拟必须进行两次（对于最热控制模式和最冷控制模式）"。这意味着设计人员在方案分析阶段，进行方案比较时，不得不回到概念性阶段，再次进行建筑模拟计算。建筑模拟和系统模拟之间的联系无法体现出设计过程中两两个阶段之间的关联。而采用房间负荷作为各个阶段之间的联系导致建筑模拟、系统模拟和机组模拟等各模块过度紧密地耦合在一起，这使得DOE被限制在建筑物全年能耗分析，而不能胜任设备选择以及管网系统校验等工作。 由此可见，上述两种模拟分析工具都存在着某些缺点而无法有效地应用在设计过程中。因此，制作一个适用于设计的模拟分析工具，必须充分考虑设

7、计过程的阶段性；处理好各个设计阶段中的已知、未知关系；设计过程应考虑全年的运行状态，因此必须采用另一种运行方式来替代实际的小步长控制方式模拟。 作为ANNEX30的一个参加者，清华大学提出了"分阶段设计，分阶段模拟"的思路，在充分考虑上述3个要素的基础上，开发出了建筑环境控制系统模拟分析工具包（DeST），并应用在若干实际工程中。DeST是基于功能的模拟软件，对应设计的不同阶段，提供相应的功能性模块。其任务是在设计的整个过程中，通过建筑模拟、方案模拟、系统模拟、水力模拟等手段对设计进行校核，并根据模拟数据结果对设计进行验证，从而保证设计的可靠性。 2 DeST的结构DeST

8、在设计时充分考虑了"设计的阶段性"这一特点。相应于设计的不同阶段，DeST由不同的功能性模块组成，并根据阶段之间的联系在模块之间建立其相应的关联。图1是DeST的整体框架结构示意图。矩形表示各个模块，箭头表示各模块之间的关系。 图1 DeST结构框架示意图 DeST所需要的气象数据由Medpha产生，其基础是20年的实测数据和随机气象数学模型。目前Medpha可以生成各格式的、193个中国城市的逐时气象参数。计算机辅助建筑描述程序CABD是一个基于ACAD平台的建筑描述界面，设计人员通过它描述建筑物的围护结构（几何尺寸，热特性参数）以及各种内扰的变化情况。在进行详细的建筑模

9、拟时，需要输入各种经验系数（例如热量在空间内的分布等），这通过经验系数维护程序ECM完成。CABD是DeST的主控界面，它把用户绘制的建筑物的相关数据自动传输给建筑分析模拟模块BAS。BAS的任务对建筑物进行详细的逐时模拟，其数学模型是增强的状态空间法6，7。BAS是一个精确的多空间建筑模拟程序，它负责计算逐时的房间基础室温（RBT，在没有任何空调系统影响下的房间温度）。逐时的基础室温反映了房间在被动热扰影响的下的热特性。在初步设计阶段，建筑师可以通过基础室温来比较各种因素的影响，如围护结构的材料、朝向、建筑物的形状等等。当建筑设计确定之后，方案模拟程序Scheme可用来计算建筑物在各种空调方

10、案（分区，系统类型，运行方式）下的热特性，在方案设计阶段，设计者可以通过模拟结果对不同的空调系统方案进行比较取舍。在方案确定之后，方案模拟程序计算出对机组或者末端的详细要求，通过逐时系统要求的送回风参数以及风量，空调机组选择程序ACSel对选择的设备进行全工况满足要求，另一方面可以检验各设备在全年工况下是否能完全满足需求，另一方面得到对冷热源的水温、水量要求。当对冷热源的需求明确后，类似的方法可以用于冷热源的需求明确后，类似的方法可以用于冷热源的优化选择，通过冷热源优化程序CPO对用户选择的冷冻机类型、台数和运行方式进行校核，保证机组在整个运行周期内保持最高的能效比。通过方案模拟得到全年逐时要

11、求的风量后（对于变风量系统），再通过送风管网可及性分析DNA计算出风机全年的工况点，从而可以根据其全运行要求选择风机，使其大部时间工作在高效率区间内。同时，通过可及性分析也可计算出各末端要求的压差，以此通过NLA对变风量末端进行噪声分析。类似的策略可以用于水管分析，通过PNA来实现。图2列出了设计的各阶段及其对应的DeST模拟程序。 图2 DeST各模块及其对应的设计阶段 通过这样的结构设计，设计人员在每一个阶段都能利用相应的模拟模块来计算不同设计中系统的性能（满意度、能耗要求），并通过比较确定较佳的方案。同时，本阶段模拟的一部分结果也是下一阶段设计的输入（对下一阶段的需求。）作为一个服务于设

12、计者的工具，DeST根据设计者的要求进行繁复的计算，而设计得通过分析模拟结果对设计进行比较取舍。 3 对已知和未知条件的处理设计过程包含各种不同的设计阶段，每个阶段的已知和未知条件不同，随着设计的展开，各阶段的已知和未知条件也随时之相互转化，前一阶段的未知因素通过设计成为本阶段的已知条件。例如，在初步设计阶段，内外扰是已知条件，在这些扰运作用下建筑物的热特性是未知的；而到了方案设计阶段，建筑物的热特性成为已知因素，设计者需要在详细的建筑物热特性的基础上对空调方案进行比较、取舍，并为进一步的设备选择提供依据。建筑物的热特性是初步设计和方案设计之间的重要桥梁，通过设计分析，它从前一阶段的未知条件变

13、为后一阶段的已知条件。 在每一个设计阶段，DeST采用详细的数学模型来表述已知的部分，而"理想化"的部件来表述未知的部分。假设"理想化"的部件能满足任何的要求（冷热量、水量等）。这样的处理与设计过程相当吻合，并且避免了"缺省的部件"对模拟结果的不利影响。因为有些未知的部分往往是到下一阶段才能解决，无论采用何种的"缺省部件"都不能保证与下一阶段最终选择的部件一致。采用"理想化"模型具有两个优点： ·基于"理想化"模型的模拟结果具有可比性，因为它们采用了相同的输入和假设

14、。 ·可以得到对下一阶段的需求。"理想化"模块的输出便是对实际设备的要求，而"缺省设备"则无法为下阶段选择提供有益的信息。 以方案设计阶段为例，表1中给出了本阶段的任务、已知和未知条件。 表1 方案分析阶段的任务和条件 任务选择适当的分区方式、空调系统形式（变风量、定风量等）以及运行方式（新风利用策略等） 已知建筑物的热特性（基础室温、各种热扰对房间温度的响应） 未知送风管道系统空调机组的详细信息控制手段 方案模拟数学模型为 ，每一时间步长 式中，tk为房间 k的空气温度；dk为房间k的含湿量；Qk为通过风机盘管或末端再热器投入到房间k的冷热量

15、；Gk为房间k的送风量；tk,base为房间 k的基础室温；A，Aj，Bj为状态空间法中的系数；tsupply为机组的送风温度；d supply为机组的送风含湿量；t k,set,min，t k,set,max分别为房间 k的温度设定最小值和最大值； k,set,min， k,set,max为房间 k的相对温度最小值和最大值，%；G k, min，G k, max 为房间k的最小和最大送风量（变风量系统）；Q k, min，Q k, max 为房间k的风机盘管或者末端热器可投入的冷热量最小值和最大值；t out，d out 为房间k分别为室外的空气温度和含湿量；t return，d retur

16、n 分别为机组的回风温度和含湿量；min：J为空气处理室及各末端的最小能耗。 通过详细的建筑模拟，得到了各房间的基础室温tk,base和空调系统（冷热量）作用在房间温度上的影响系数A和Aj。它们是本阶段的已知因素。各个时刻下送风量（Gk）、送风状态（tsupply，d supply）以及风机盘管或者末端再热器的输出（Qk，如果系统中有相应的设备）需要通过本阶段的模拟来确定，其目标是使系统内所有的房间都能满足其设定值要求，所有的参数都必须满足方程组（1）的限制。通过下列步骤，可以计算出各种空调方案下系统的性能： 计算每一个房间所要求的送风状态区域（SSRk）； 由于在任一时刻，系统只能存在一个统

17、一的送风状态，因此需要求出所有房间送风状态区域的交集（公共的送风状态区间，CSSR）； 以公共的送风状态区域内的任何一点作为送风参数，都能使分区内所有的房间满足其设定值要求，而不同的送风状态点对应的空气处理能耗是不同的，在此通过某种优化算法计算出公共送风状态区域中的最优点，同时确定产生此最优送风状态点的相应的空气处理过程： 在确定最优的送风状态之后，依据能耗最小的原则，可以确定系统的送风量及风机盘管或者末端再热热器投入的冷热量。 计算出各个房间的温湿度。 根据计算出来的逐时的各个房间的温湿度，统计全年内各房间满足设定值要求的小时数或者比率，并以此来比较不同空调方案的性能。如果不同的方案都能满足

18、各房间要求，则通过各方案要求能耗值来进行比较。 在此阶段，空气处理室被发作一个"理想化"的设备，假设它能够产生要求的任意送风温湿度。在进行方案模拟时，只需要确定空气处理室的类型，而不需要确定各组成的详细参数。图3是一个"理想化"空气处理室的示意图，该空气处理室的示意图，该空气处理室包括一个混合段、表冷盘管、再热盘管以及蒸汽加湿器。t1，t2代表室外和室内的空气状态，tW是冷冻水进口温度。当空气处理室的类型确定后，类似图3的焓湿图分区也随这确定。空气处理过程的方向受到设备类型的限制，例如加热过程必然是一根垂直向上的线。而设备的容量是无限的（"理想

19、化"假设），空气处理过程线的长度可以任意长。采用"理想"化的空气处理室进行方案模拟，可以计算出全年中各个时刻对各个部件的要求。当方案确定后，对整个空气处理室的要求随之确定。根据逐时要求的各段空气处理过程线，在详细设计阶段可以对机组的各个组成部件进行详细的校核。 图3 一个“理想化”的空气处理室 同样的假设也被用在送风管网分析上。在方案分析阶段，假设送风管网可以提供任意要求的送风量。当方案设计完成后，也得到了要求的逐时风量分布数据，而这些数据正是进行风机和管网的详细校核所需要的。 4 控制和逆向计算过程通常控制都是以小步长进行的，但在空调系统设计时，需要考虑建筑物和

20、系统全年的运行情况。如何将这两种不同类型的过程结合在一起呢？实际上，设计可以划分成两个层次：空调系统设计和控制设计。DeST注重于解决前一层次的问题。无论选择何种系统，采用何种设备，系统设计的目的是要产生一个完全可控的、能够满足用户要求的系统。DeST对系统进行模拟时，以1h为时间步长进行长时间（全年）的计算。为了避免不同控制器特性的影响，没有采用小步长的控制方法，而用逆向的计算过程。例如，在详细设计阶段，当对变风量系统的送风管网进行分析时，设计者的任务是校验管网能否满足各个时刻的风量分布要求，并选择适当的风机。 为解决此问题，采用传统小步长控制的模拟过程是： 选择一个风机，设定控制参数（比如

21、PID参数）； 计算管网各处的流量； 如果管网的流量与要求的流量没，通过某种控制策略调整变风量末端； 在下一个小的时间步长内，重复，以期达到要求的流量分布。 由于风道的惯性非常小，此模拟必须以相当小的时间步长进行计算（1s），能否达到要求的流量分布与控制策略有很大的关系，因此此方法不适用于对全年各种工况的校验。从另一个角度考虑此问题，假想在一种理想化的控制下，各变风量末端可以满足要求的送风量，从而可以计算出各管段上的流量；假设各房间的压力为零，则当定压点压力能够维持时，各变风量前后的压差可以计算出来，对风机的流量和压力要求也可以确定。相对于传统的模拟过程，这是一种逆向的求解过程。因此，校验送风

22、量分布是否能够实现，应该首先通过在各工况下地风机的要求来判断，而不应立即着眼于确定详细的控制策略。只要风机和管网能够在理想控制下满足要求，则必然可以通过某种具体的控制方式满足流量分布。在设计时，应先解决"可控性"的问题，然后再解决"如何控制"的问题。图4给出了在两种定静压控制方式和一种变静压控制方式下风机的工况点。在此基础上可以选择风机以满足全年运行，结合风机的性能参数，确定风机全年能耗，也可以根据各变风量末端的压差计算出各个时刻下的末端噪声。 图4 不同静压控制方式下的风机工况点定压点位于风机出口处浮动静压控制定压点位于风道上距风机出口2/3处 采用逆

23、向的求解思路，避免了采用小步长的反馈控制，逆向的求解过程可以看成一种开环的"理想化"控制方法。用此方法在进行设备（例如表冷盘管）校核计算时，可以计算出已知出口和入口参数时对冷冻水侧的要求。如果校核发现该设备无法达到要求的出口状态，则无论采用何种控制方式，该设备都无法满足运行的要求。从此意义上说，通过逆向的求解算法，DeST着重研究系统的可控制性，即： 该设备能否通过某种控制方式满足要求？ 如果可以，该设备的最佳运行效果是什么？ 通过校核回答以上两个问题后，设计人员可以进一步研究具体的控制方法，并通过与最佳的运行效果进行比较以确定控制方法的优劣。 5 DeST的用户界面DeS

24、T在WINDOWS95/98/NT下运行。所有的模块都集成到CABD中。CABD是一个基于AutoCAD R14开发的用户界面（见图5）。用户在此界面上进行建筑物的描述，通过选单调用其它模拟模块，与建筑物相关的各种数据（材料、几何尺寸、内扰等）通过数据库接口与CABD相连。采用CABD绘制如图5所示的建筑物大约需要13h。建筑模拟所需要的逐时气象参数，通过利用Medpha（图6）选取相应的城市即可获得。图5 CABD的界面形式 图6 Medpha(逐时气象数据生成程序) 完成建筑物内外扰的描述之后，BAS负责对建筑物进行详细的模拟计算，并同步在如图7的示的界面中显示逐时的基础室温以及月平均温度

25、。在方案模拟的进程中，空气处理过程可动态地在焓湿图上显示，让用户了解在DeST中各种参数（送风温度、送风量等）是如何确定的。 图7 BAS(建筑分析及模拟程序) 各种模块以ActiveX、DLLs通讯ARX的形式集成在一起，使得DeST成为一个高度集成化的软件工具，其目的是最大限度地减少用户花费在输入数据的时间，让设计人员将注意力集中在分析模拟结果、比较方案等创造性的工作中。所有的模拟结果以纯文本的格式存储，用户可以很方便地使用其它数据处理工具（如EXCEL）进行整理和分析。通过一定的实践，当用户能够熟练地使用DeST的界面后，准备数据以及运行程序所消耗的时间大概是分析所需要的时候的1/4或者

26、更少。          6 DeST能够解决的问题下面用两个实例来演示如何采用DeST 设计的不同阶段进行分析。其一是ANNEX30案例1中1个9层的办公楼，该楼位于德国，在模拟时采用了比利时的室外气象参数。设计要求房间温度全年控制在2226，相对湿度必须满足40%60%的范围。该案例的详细数据参阅文献10。图8该建筑的标准层平面图。另案例是位于天津的一个商业建筑。 图8 ANNEX案例1的建筑平面图 6.1 初步设计阶段 在本阶段，通过DeST计算出不同朝向下各房间逐时的基础室温，对该建筑物的不同朝向进行了比较

27、。图9、10分别是不同朝向下（正面朝南：朝向一；正面朝东：朝向二）两房间全年的基础室温分布。图9 休息室1的基础室温分布 图10 会议室1的基础室温分布 从上述结果可以看出，朝向对于类似会议室1的房间没有太大的影响，而对于类似于休息室1的房间，则有显着的影响。朝东时该类房间要比朝南时温度偏高许多，这说明太阳辐射对于此类房间是一个很重要的影响因素，通过此比较可对建筑的最初设计提供参 考。 6.2 方案设计：水系统类型比较 图11是位于天津的一个宾馆的示意图。该楼计划 4个朝向都采用风机盘管系统。为满足客户的要求，房间温度必须能够在2228之间任意调整。因此，不同的控制精工要求对于风机盘管制供水需

28、求有很大的影响，尤其在过渡季。控制精度要求越高，同时需要提供冷水和热水的时间越长。 图11 天津某商业建筑物的标准层示意图 在模拟时，考虑了3种不同的控制精度，表2是3种不同的精度及其对应的设定值范围。 在每一种控制精度下，DeST对该楼进行了全年的模拟，计算出各月份需要冷热源同时提供冷水和热水的小时数，表3是过滤季时结果。 表2 3种不同的温度精度要求 精度房间的设定值±121±1 23±1 25±1±222±2 24±2 26±2±323±3 24±3 25±3 表3需要

29、同时提供冷热水的运行时间/h 精度3月4月5月热冷热冷热冷±1744304604450345709±2744175579204272525±366513315479288精度9月10月11月热冷热冷热冷±1256684589544720149±216956957132372075±353663665366811注：“热”，“冷”分别表示需要供热水和冷水。 从结果可以看出，当要求的精度较高时，风机盘管需要设计成四管制以满足各个时刻同时的冷热水要求，否则在过渡季中将有上千h不满足。如果要求的精度不高（如±3），则两管制的系统基本可以满足