第四节 软件选取

一、火灾模拟

(一)概述

火灾数值模拟是火灾研究的重要内容之一，但由于火灾现象的复杂性，近几十年来才建立起描述火灾现象的实用数学模型。火灾模型主要分为确定性模型和随机性模型。

火灾数值模型主要有专家系统(Expert System)、区域模型(Zone Model)、场模型(Field Model)、网络模型(Network Model)和混合模型(Hybrid Model)。场模型也即CFD(计算流体动力学)模型，主要是指利用计算流体动力学技术对火灾进行模拟的模型，由于CFD模型可以得到比较详细的物理量时空分布，能精细地体现火灾现象，加之高速、大容量计算机的发展，使得CFD模型得到了越来越广泛的应用。

目前用于火灾模拟的CFD模型主要有FDS、PHOENICS、FLUENT等。FDS是专门针对火灾模拟而开发的CFD软件，简单易用。因此，它在火灾模拟中的应用最为广泛。而PHOENICS和FLUENT是计算流体力学的通用软件，将其用于火灾模拟需要有较强的流体力学背景，因此，应用较少。目前，国内外对FDS的研究比较多，而对于PHOENICS和FLUENT在火灾模拟方面的应用研究较少，对各个软件的对比研究则更少。

在火灾模拟中，怎样才能使模拟结果更加准确、可信是一个急需解决的问题。验证与确认是评价数值解精度和可信度的主要手段。长期以来，CFD工作者对CFD软件的验证与确认工作一直没有给予足够的重视。因此，对于计算结果的可信度，CFD研究人员并不能给出明确的回答。这使得CFD软件的使用者对CFD也持一种矛盾的心态，既想利用CFD这种快捷经济的设计工具，又对CFD的计算结果心存疑虑。如果有条件，可以结合数值计算和模拟实体火灾的方式，进一步验证模型的可靠性。

(二)选取

从软件易用性来看，火灾专用模拟软件相对简单，在应用中不需要作复杂设置，使用者只需掌握火灾基本知识即可得到合理的结果，而通用CFD软件对使用者要求较高，使用者需要对流体力学知识有深入了解，才能得到合理结果，因此，一般火灾模拟选择专用软件为宜。

利用火灾模型进行数值分析前，应着重考虑该模型对所模拟问题的适用性及预测能力，一般情况下，需要事先利用相关试验(已有其他人员进行的试验或自己进行相关试验)对模型进行确认研究。从模拟结果的准确性来看，火灾专用模拟软件由于是专门针对火灾开发的，在概念模型层面相对于通用软件更接近于真实模型，其数学模型更能反映火灾过程，因此，一般情况下，建议选择火灾专用软件，除非在专用软件无法模拟的情况下才选择通用软件。

使用火灾专用软件时，应着重考虑网格独立性、边界条件设置对模拟结果的影响，使用通用软件时，还应考虑湍流模型、燃烧模型、辐射模型的选择。

火灾模型的验证和确认应包含其对各类火灾参数的预测能力研究，如火场温度、热辐射通量、反应产物的浓度变化(着重研究C0、C0，、烟密度等)、火场能见度等。

对于通用的CFD软件，如PHOENICS、FLUENT、CFX等，由于其发展比较成熟，其程序一般能够比较准确地反映其所确立的概念模型，因此，对这类模型可以着重于确认研究;对于专用火灾模拟软件，如FDS等，已经进行了较多的确认和验证工作，对于比较常见的火灾场景，如建筑室内火灾等，可以直接用来模拟分析，而对于一些特殊的场景，如火灾在狭长双层玻璃幕墙内的蔓延模拟，还需进行进一步确认研究;对于自行编制的火灾模拟程序，模型的验证工作是至关重要的，应确保程序能够准确反映概念模型。

利用Steckler房间火灾试验对PHOENICS和FLUENT进行确认研究，就该类实例来说，FLUENT的准确度要高于PHOENICS，但就工程应用来说，在选择合理的湍流模型、辐射模型，并经过网格独立性检验后，二者的模拟结果一般均可满足工程需要。

火灾发展具有确定性和随机性的特点，火灾试验的影响因素较多，在选择确认试验时，应尽量选择可重复性强的试验，并应注重采用不同火灾场景下的火灾试验对其进行确认研究，以便更好地检验模型的可信度。

二、疏散模拟

(一)概述

(二)疏散模型分类

人员疏散计算方法主要有两种：水力疏散模型和人员行为模型。

1.水力疏散模型

水力疏散模型通常对人员疏散过程作如下保守假设：

1)疏散人员具有相同的特征，并且都具有足够的身体条件疏散到安全地点。

2)疏散人员是清醒的，在疏散开始的时刻一起井然有序地进行疏散，且人员在疏散过程中不会中途返回选择其他疏散路径。

3)在疏散过程中，人流的流量与疏散通道的宽度成正比分配，即从某一出口疏散的人数按其宽度占出口总宽度的比例进行分配。

4)人员从各个疏散门疏散且所有人的疏散速度一致，保持不变。

2.人员行为模型

下面将通过分析这些人员疏散模型的功能与特点，对这些软件进行适当分类。

(1)一般分类。疏散模型在处理疏散的一般问题时，均采用了三种不同基本方法：优化法、模拟法和风险评估法。

优化法假定人员以最有效的方式进行疏散，而不考虑外部环境的影响及非疏散行为。通常，模型认为人员选择的疏散路线是最佳的。这一类模型适用于大量的人群或将所有人员当做一个有共同特性的群体来考虑的情况，而不考虑个体行为。

模拟法试图表现实际的疏散行为与运动，不仅要得到准确的结果，而且要反映疏散时选择的疏散路线及人员所做的决定。由于各个模型在考虑人员行为时的详细程度不同，因此结果的准确度也不相同。

风险评估模型能识别出火灾时与疏散有关的危险或相关事故，并能对最后的风险进行量化。通过多次重复运算，可以估算出与不同防烟分区设计或防火保护措施有关的各种重要变量的统计数据。有关模型的类别与名称见表4-3-1。

(2)建筑空间的表示。各种疏散模型都必须对建筑空间进行描述，以模拟人员在建筑内部的疏散过程。在模型中，空间被划分为许多小的区域，每个区域都与相邻的区域相连。根据对空间划分的精细程度，常将模型中的空间划分分为两种方法：精细网络法和粗糙网络法。

对于精细网络法，整个封闭空间用覆盖一些瓦片状的网格来表示，各个模型中节点的网格大小与形状都不同。例如，Exodus采用0.5m×0.5m的正方形网格节点，Simulex则采用0.2m×0.2m的正方形网格节点，而Egress则采用六边形的网格节点，每个网格的大小足以容纳一个人。这些网格节点之间的连接也不相同，在Exodus中每个网格与相邻的八个网格节点相连，而Egress中每个网格与相邻的

六个网格节点相连。因此，Egress的方法可以准确地表示封闭空问的几何形状及内部障碍物的位置，并在疏散的任意时刻都能将每个人置于准确的位置。

对于粗糙网络法，空间的描述是按照实际建筑结构的划分来确定的，每个网格节点表示一个房间或走廊，然后根据它们之间的实际连接关系构建其网络模型。在这类模型中，仅能表示人员从一个建筑单元移动到另一个建筑单元，而无法描述人员在一个建筑单元内的运动，它也无法处理一些局部的现象，例如超越前方人员、避开障碍物等。有关不同空间划分的模型分类见表4-3-2。

(3)人群分析。各类疏散模型在对人员进行分析时，采用了两种方法：个体分析法和群体分析法。个体分析法允许用户设定或由随机方式确定个体特性，人员决策与运动由这些个体特性决定。需要注意的是，不能将个体的独立决定与不能执行群体行为混为一谈，定义个体时并不排斥他具有群体行为，而是先考虑每个人的个体特性，然后再为他指定一个行为，而这个行为也许就是群体行为。群体分析法将人群视为一个具有共同特性的群体。在描述疏散过程时，不针对逃生的个体，而针对大量的人群。这种方法难以模拟事件对个体的影响(例如：火灾烟气毒性的影响)，而只能对整个人群的普遍影响进行模拟。例如，它不能表示老年人或残疾人等特殊人群的生存率，而只能表示受影响的人的比例。它的好处是模型的运算速度相对较快。有关基于人群分析的模型分类见表4-3-3。

(4)行为分析。人员在逃生时的决策过程是复杂的，疏散模型根据模拟人员决策过程时所采用的分析方法可分为以下几类：无行为准则模型、函数模拟行为模型、复杂行为模型、基于行为准则的模型以及基于人工智能的模型。

无行为准则模型完全依赖于人群的物理运动和几何形状的物理表达，来影响人员的疏散，并对其进行预测判断。

函数模拟行为模型把人员的行为用一个方程或一个方程组来描述，以此达到控制人的响应的目的。

这类模型可以将人定义为个体，但由于所有个体均受到同一函数相同的影响，且会以一定的方式对这种影响产生反作用，因此实际上削弱了个体行为。该函数或者按照现实生活中人员的行为来建立，或者引用其他从事人体行为模拟研究领域的成果(例如，磁模型的方程来源于物理学)。复杂行为模型通过复杂的物理方法来含蓄表示行为决策准则。此类模型一般基于第二手数据的应用，包括心理的或社会的影响，因而它依赖于第二手数据的准确性与有效性。基于行为准则的模型预先规定了一套人员的行为准则，然后再根据这些准则来确定疏散过程中人员的行为。例如，“假如人在一个充满烟气的房间里，他会通过最近的出口离开”等类似准则。但是，这种行为决策方式会导致人员在相同的环境下以某种确定的方式进行反应，从而与实际中的人员反应有所差异。

人工智能模型将个体人员设计成能对周围环境进行智能分析的模拟人或与之相近的智能人，因此可以准确地表现其决策过程，但这会使用户对人员行为的控制权被计算机所代替。疏散模型按人员行为分类见表4-3-4。

3.人员行为特性

火灾是具有突发性的意外事件，伴有火焰、浓烟、强烈的热辐射、噪声和有毒气体，常在短时间内给人以毁灭性的伤害。身处火场的人们往往需要承受巨大的心理压力，从而表现出各种各样的异常行为。研究发现，不同的心理素质、阅历和经验，会导致人在遭遇火灾时，呈现不同的心理反应和行为。但是，如果在遭遇火灾时，能保持良好的心理状态，及时采取自救行动，往往能够化险为夷，成功疏散，避免伤亡。

对于人员特性的考虑可以分为两方面，即单个人员独立考虑和全局考虑。大多数模型可以根据用户要求或计算机自动设定每个人员的移动属性(如步行速度)，并记录每个人员在任何时刻的移动历史轨迹。这类模型也不排除群集行为特性，但它是按单个人员检查和分配各自的移动特性的，它需要较多的计算机容量，程序处理的难度也稍大。另一类是按群集方式来考虑人员的移动特性，它将一群或一组人群按同一特性考虑，即将一群人按同一移动速度考虑，认为他们同时到达或离开建筑物的某个网格节点，它具体计算建筑内人员疏散成功率，其操作简单，使用方便，运行速度也较快。

4.软件介绍

(1)STEPS。该模型可用于模拟在正常或紧急情况下，人员在不同类型建筑物中的疏散情况。

模型是一个由一系列的网格单元组成的网格系统，在网格系统中，一个人只能占有一个单元。网格单元的默认尺寸是0.Sm×0.5m。另一个细网格选型可适用于多人占有一个网格单元，但仍处于测试阶段。

(2)Simulex。该模型是一个能够模拟人群从复杂建筑物中疏散的模型。该模型采用一个连续的空间体系，各层的平面图和楼梯都划分成一个个0.2m×0.2m的块或网格。该模型包含一个算法，它能够计算出每个网格到最近安全出口的距离，并将这些信息标注在一个距离图表上。

(3)SGEM。该空间网格疏散模型可利用CAD平面图生成复杂建筑的疏散图案，比较后得出最佳疏散设计路线。这个模型已经用于一些咨询项目。

模型的结构基本上是细网格模型。最初，此模型将建筑物分成一些节点，这些节点代表建筑物的空间或区域(不受保护、部分保护和全部保护)，其中区域之间至少有一个弧形口连接，由此形成一个粗的网格。然后，将每个粗的网格单元再分成有限的网格，每个网格的大小是0.4m×0.4m，一个人占据一个网格。此外，同一时间内一个人只能占据一个网格单元。

(4)buildingExODus。该模型可用于模拟疏散大量被很多障碍围困的人。模型由airEXODUS、

buildingExODus、maritimeEXODUS、railEXODUS、vrEXODUS六个部分组成。buildingEx0Dus试图考虑人与人、人与火灾以及人与建筑物之间的相互作用。该模型包括6个在模拟疏散方面相互联系、相互传递信息的子模型，即人员、运动、行为、毒性、危险性和几何学子模型。

该模型是一种行为模式模型，是一个细网格系统。其利用二维空间网格绘制出几何结构、位置、障碍物等。这种网格由“节点”和“弧”组成。每个节点都代表了建筑平面图上的小空间，而弧在建筑平面图上把这些节点连接在一起。人们通过利用这些弧从建筑物的一个节点到另一个节点。这些信息存储在几何子模型中。同时，在整个模拟过程中，每个节点都有毒气等级、烟气温度和浓度等与其相关的动态环境。

5.软件选用

疏散软件能够人为设置出口障碍，通过对建筑平面信息的识别，可以在某些安全出口受阻的情况下创建另一个替代的有效距离地图来引导人员疏散，从而得出建筑物的最优化疏散设计方案。人员的行为特性具有不确定性，无论是采用函数型行为模式，还是采用规则决定方法，都无法绝对准确地模拟出火灾时人员疏散的真实情况，在今后的软件后续开发中可以增加更多的实际因素参数，以修正模拟疏散结果。

疏散软件基本可以定义个体人员的基本行为特性，模拟出一定的心理反应，能够有条件地寻找出有效的疏散路径。国际上人员疏散软件的开发已经积累了一定的经验，并且经过了相当数量的工程实践检验，但毕竟国外的人员疏散行为、习惯，心理反应等与我国居民有很大差异，一些软件还不能简单地直接拿来应用，需结合我国的建筑特点和疏散演习等进行校正和应用。一般情况下，当建筑的结构简单、布局规则、疏散路径容易辨别、建筑的功能较为单一且人员密度较大时，适合采用水力疏散模型来进行人员疏散的计算，而其他情况则适于采用人员行为模型。原因是：前一种情况下发生火灾，人员可以较为快速准确地获取关键信息并采取疏散行动，人员疏散通常很快由个体行为转化为群体行为，符合水力疏散模型的适用条件。而其他情况下，人员获取信息的渠道相对较少，从接到报警信号到识别火灾以及疏散过程中寻找合适的疏散路径都在很大程度上依赖于个体的判断。因此，人员疏散由个体行为转化为群体行为的时间非常长，有时候可能自始至终都是个体行为，这种情况下选用人员行为模型显然更合适。如果有条件，可以结合数值计算和建筑疏散演习等方式，进一步验证模型的可靠性。

三、模型评价

下面针对建筑消防性能化设计计算方法的确定性模型需要重点评价的几个方法进行论述。

(一)计算模型的适用性

以火灾动力学软件FDS为例，它可用来模拟火灾热和燃烧产物的运输、气体和固体表面之间的辐射和对流传热、热解、火灾蔓延与增长、喷淋等。对于开放空间或燃料控制的火灾，FDS能相对准确地进行模拟。但FDS的局限性在于其限于低速流动模拟;通过分解压力项，处理状态方程，从而滤除声波的影响。针对相对封闭房间内氧控制的火灾场景，有可能会发生爆燃现象，在此过程中压力波对火焰的传播起着较大的影响。在模拟此类火灾场景时，尽管FDS能模拟并有可能获得看似正确的计算结果，但从模型基本理论上已不再适用。因此，针对计算模型的适用性问题，不仅要从计算结果进行考虑，还要从模型的自身假设进行分析。计算软件为了能够模拟更多的问题，往往采用普适性的算法，对于有些根本不满足计算模型理论的场景，计算结果也可能会与实验结果偏差不大，这样的结果不能轻易相信，且也不能说明类似这样的场景就可以采用这样的方法来计算。在计算模型理论都无法满足的情况下，根本就不允许采用这样的模型来进行计算。

(二)计算的收敛性

在数值方法中，需要对连续性的数学模型进行离散化后再求解，也就是用一个离散的数值模型来近似模拟。时间和空间都需要离散化。一个连续性的数学模型有很多不同的离散方法，形成很多不同的离散模型。为了获得一个好的近似解，要求离散模型能够模拟连续模型的性质和行为。这就要求离散方法采用高阶精度的格式，同时要保证其不会带来计算结果的非物理振荡，能更好地收敛于真实解。

对于定常模拟来说，只需要求最终的计算结果逼近真实解。但对于非定常模拟来说，则要求每一计算时间步内的结果也要收敛，且要达到能接受的计算精度。如果模型没有发生时间步的截断而且能保持长的时间步，那表明该模型没有收敛性问题，反之如果经常发生时间步截断，那模型计算将很慢，收敛性差。时间步的大小主要取决于非线性迭代次数。如果模型只进行一次非线性迭代计算就可以收敛，那么表明模型很容易收敛;如果需要2～3次非线性迭代计算，则模型较易收敛，如果需要4～9次非线性迭代计算，则模型不易收敛，大于10次的模型可能有问题。

影响计算收敛性的因素很多，如网格尺度、计算格式精度、初始流场参数、化学反应的刚度、计算模型等。

(三)网格尺度的合理性

对于建筑火灾场模拟计算，首先应该考虑网格尺度的合理性问题，而这一问题也是场模拟计算中非常重要的问题。网格尺度的合理性问题直接影响计算结果的误差，甚至影响计算结果是否定性合理。

网格尺度的合理性一方面是计算结果不依赖于网格尺度的变化，即网格的独立性;另一方面，在保证网格独立性的同时，应考虑计算资源的能力，尽可能减少计算量，提高计算网格的经济性。在场模拟计算中，如何做到这两点呢?

1.网格独立性

通常的做法是，下一次要考虑的网格尺度一般为前一次网格尺度的1/2，即网格加密一倍。如果加密一倍的计算结果与该次加密前的计算结果之间的误差在可接受的范围内，网格不再加密，即可采用该次加密前的网格尺度的计算结果作为最终结果来进行分析评判。如果加密一倍的计算结果与该次加密前的网格尺度的计算结果之间的误差不在可接受的范围内，应进一步进行加密。当然，加密的起点也应有一定的基础，可以基于计算者的经验、基于模型分析、基于计算问题的分析、基于前人或公开发表类似问题的经验等。基于这样的基础，可以加密，也可以加粗网格。总之，针对具体的问题，也不一定遵循前述加密原则，可适当增大加密强度。

2.网格经济性

一般，加密一倍网格，计算量增大8倍，计算时间可能增大几十倍，甚至上百倍。一方面要保证一定的计算精度，另一方面要考虑合适的计算量。因此，采用能满足该精度的最粗网格，也可以采用局部加密技术，在高密度梯度区(如火源)、壁面附近等加密网格，在低密度梯度区或影响相对较小的区域加粗网格。

(四)时间步长的合理性

在求解微分方程时，必须注意时间步长的选择。首先应考虑系统的稳定性。在分析和求解瞬态算法时，为了解的收敛，必须考虑稳定性。对时间步长进行限制的算法，称为有条件稳定。未对时间步长进行限制的算法称为无条件稳定。在求解连续性问题0DE的解析解时，稳定积分能给出衰减解。对于某些时间步长，不稳定方法会产生无界或快速振荡的数值解。要意识到即使是稳定连续性模型，数值模型也有可能不稳定。因此，原连续性模型不稳定时，任何数值模型都得不到精确解。相反，在条件不稳定的情况下，无条件稳定的算法也能够得到稳定的数值模型。这意味着无条件稳定的算法不能考虑快速增长的现象。

在建筑性能化设计计算的火灾场景模拟中，时间步通常是条件稳定。时间步过大，会出现数值振荡，进而导致不收敛，计算不能进展下去。时间步一般满足流动的CFL条件，如FDS中的时间步dt=5