Реферат: Методика определения расчетных величин пожарного риска в объектах защиты на основе полевого метода моделирования пожаров Общие положения

МЕТОДИКА

определения расчетных величин пожарного риска в объектах защиты на основе полевого метода моделирования пожаров


Общие положения


1.1 Настоящая методика определяет порядок получения расчетных величин пожарного риска на основе методов компьютерного моделирования для зданий, сооружений, производственных объектов, строений (далее – объект защиты), которые соответствует Ф1 – Ф5 классам функциональной пожарной опасности. Данная методика может применяться для определения расчетных величин пожарного риска при проектировании, реконструкции, реставрации, капитальном ремонте и эксплуатации объекта защиты. Использование данной методики рекомендуется совместно со специализированным пакетом программ для компьютерного моделирования пожаров «σFire».

Применение данной методики не предполагается использовать для проведения расчетов пожарного риска объектов защиты специального назначения, в том числе объектов военного назначения, объектов производства, переработки, хранения радиоактивных и взрывчатых веществ и материалов, объектов уничтожения и хранения химического оружия и средств взрывания, наземных космических объектов и стартовых комплексов, горных выработок, объектов, расположенных в лесах.

1.2 В соответствии с постановлением Правительства РФ от 31 марта 2009 г. №272 «О порядке проведения расчетов по оценке пожарного риска» расчеты по оценке пожарного риска проводятся путем сопоставления расчетных величин пожарного риска с соответствующими нормативными значениями пожарных рисков, установленными Федеральным законом "Технический регламент о требованиях пожарной безопасности" (далее – Технический регламент).

Определение расчетных величин пожарного риска осуществляется на основании:

а) анализа пожарной опасности зданий;

б) определения частоты реализации пожароопасных ситуаций;

в) построения полей опасных факторов пожара для различных сценариев его развития;

г) оценки последствий воздействия опасных факторов пожара на людей для различных сценариев его развития;

д) наличия систем обеспечения пожарной безопасности зданий.

1.3. Расчетные величины пожарного риска являются количественной мерой возможности реализации пожарной опасности объекта защиты и ее последствий для людей.

Определение расчетных величин пожарного риска основывается на расчете индивидуального пожарного риска для людей находящихся внутри объекта защиты. Индивидуальный пожарный риск определяется как частота воздействия опасных факторов пожара (ОФП) на человека. Для определения частоты воздействия ОФП выбирается сценарий развития пожара наиболее опасного для жизни и здоровья людей.

1.4. Расчеты пожарного риска оформляются в виде отчета, в который включаются:

а) наименование использованной методики;

б) описание объекта защиты, в отношении которого проведен расчет по оценке пожарного риска;

в) результаты проведения расчетов по оценке пожарного риска;

г) перечень исходных данных и используемых справочных источников информации;

д) вывод об условиях соответствия (несоответствия) объекта защиты требованиям пожарной безопасности.

1.5. В данной методике используются основные понятия, установленные статьей 2 Технического регламента о требованиях пожарной безопасности.

Математическая модель развития пожара

На данный момент в мире все большее распространение получает подход оценки пожарной опасности строительных и технологических объектов при помощи математического моделирования.

Строительные сооружения или технологические объекты могут быть уникальными, поэтому типовые решения по обеспечению их противопожарной защиты не всегда могут быть приемлемыми. Выбор того или иного способа защиты от пожара зависит от наличия и размещения исполнительных элементов системы противопожарной защиты, а также от возможных вариантов возникновения и развития опасных факторов пожара.

Полевые модели являются наиболее мощным и универсальным инструментом компьютерного моделирования. В полевых моделях выделяется расчетная область делится на большое количество контрольных объемов. Для каждого из этих объемов с помощью численных методов решается система уравнений в частных производных, выражающих принципы локального сохранения массы, импульса, энергии и масс компонентов. С его помощью можно рассчитать температуры, скорости, концентрации компонентов смеси, тепловые потоки и т.д. в каждой точке расчетной области. Использую полевые модели можно провести расчет пожара на объекте практически любой геометрической формы с учетом основных физико-химических процессов.

Для использования полевых моделей требуются значительно большие вычислительные ресурсы по сравнению с другими менее универсальными подходами (интегральным или зональным). Однако, в настоящее время появились достаточно мощные и доступные вычислительные системы, расширились познания в области физики и химии процессов, происходящих при горении, созданы эффективные методы расчета турбулентной аэродинамики, сложного теплообмена и химического реагирования, созданы компьютерные программы, позволяющих с достаточной для практической деятельности точностью моделировать процессы, происходящие при пожарах. Проведение вычислительного эксперимента с помощью этих программ характеризуется низкими затратами и возможностью получения максимально полной информации об исследуемом объекте. В своей основе полевой метод не содержит никаких априорных допущений о структуре течения, и в связи с этим принципиально применим для рассмотрения любого сценария развития пожара.

Представленные ниже математические модели реализованы в специализированном программном комплексе σFire, созданном для компьютерного моделирования развития пожара. В программном комплексе σFire реализована поддержка многоядерной архитектуры ЭВМ, что позволяет эффективно использовать суперкомпьютеры (кластер СФУ).


2.1. Математическая модель процессов переноса

2.1.1. Уравнения для ламинарного режима течения.

Уравнение сохранения массы или уравнение неразрывности имеет вид:

(2.1.1.1)

Уравнения сохранения количества движения:

(2.1.1.2)

где -тензор вязких напряжений, - вектор объемных сил.

Составляющие тензора вязких напряжений определяются как:

(2.1.1.3)

где - динамическая (молекулярная) вязкость, - компоненты вектора скорости.


2.1.2. Уравнение энергии

В пакете “Fire” уравнение сохранения энергии рассматривается в следующем виде:

(2.1.2.1)

где - коэффициент теплопроводности, – источниковый член, отвечающий за приток (отток) энергии в процессе химического реагирования, излучения, или каких-либо других процессах.

Энтальпия многокомпонентной среды определяется по правилу смеси

(2.1.2.2)

где энтальпия компонент hm(T) вычисляется как

(2.1.2.3)

Удельная теплоемкость компонент задается в виде полинома 4-ой степени от температуры

(2.1.2.4)

Температура смеси Т в каждой точке рассчитывается по вычисленному из уравнения (2.1.3.1) значению энтальпии h и составу смеси Ym.


2.1.3. Уравнение переноса компонент среды

Уравнение переноса концентрации компонент:

(2.1.3.1)

где - массовая доля m-той компоненты,- коэффициент диффузии компонента m, – источниковый член, отвечающий за изменение компоненты в процессах химического реагирования или в каких-либо других процессах.


2.1.4. Термодинамические свойства среды

Уравнения состояния

1). Несжимаемая жидкость

(2.1.4.1)

где -константа.

2). Модель термического расширения

(2.1.4.2)

где - коэффициент термического расширения, -опорная температура, -опорная плотность.

3). Модель идеального газа

(2.1.4.3)

где P статическое давление, M - молярная масса газа, R - универсальная газовая постоянная, T - температура газа.


Свойства смеси

В программе “Fire” свойства смеси определяются через массовые доли ее компонент .

Теплопроводность смеси

(2.1.4.4)

Плотность смеси

(2.1.4.5)

Удельная теплоемкость смеси

(2.1.4.6)

Энтальпия смеси

(2.1.4.7)

Динамическая вязкость смеси

(2.1.4.8)

Молярная масса смеси

(2.1.4.9)

2.2. Моделирование турбулентности

Считается, что при помощи уравнения Навье-Стокса можно описать как ламинарное так и турбулентное течение (Ландау Л.Д., Лифшиц Е.М., 1986), практически же реальные достижения в расчетах турбулентных потоков связаны с использованием полуэмпирических моделей турбулентности использующих подход Рейнольдса. Обоснованием необходимости применения такого подхода может служить приведенная в ряде работ оценка, согласно которой отношение характерных размеров энергосодержащих (крупных) вихрей и мелкомасштабных вихрей, определяющих «турбулентный фон», имеет порядок Re3/4. Естественно, что, ориентируясь на возможности современных ЭВМ, решить поставленные задачи в практически важном диапазоне Re  103 в настоящее время маловероятно.

Применим операцию осреднения по времени к уравнениям (2.1.1.1) и (2.1.1.2) получим уравнения движения в форме Рейнольдса (RANS).

(2.2.1)

где - усредненное по времени поле скоростей, - тензор рейнольдсовых напряжений, - символ Кронекера.

Для определения компонент тензора рейнольдсовых напряжений используем гипотезу Буссинеска

(2.2.2)

где - турбулентная вязкость, - кинетическая энергия турбулентных пульсаций.

Система уравнении (2.2.1)-( 2.2.2) является не замкнутой поэтому для ее замыкания применяются различные полуэмпирические модели турбулентности которые описаны ниже.


2.2.1. Модели турбулентности

Стандартная k- модель (Launder and Spalding, 1972)

Уравнения стандартной высоко-рейнольдсовой k- модели:

(2.2.1.1)

(2.2.1.2)

где P - скорость генерации турбулентности:

(2.2.1.3)

турбулентная вязкость:

(2.2.1.4)

Тензор Рейнольдсовых напряжений принимает форму:

. (2.2.1.5)

Константы замыкания стандартной модели (Launder and Spalding, 1972):

C=0,09, k=1,0, =1,3, С1=1,44, C2=1,92.


Chen k- модель

Уравнения для модификации k- модели Chen, которые определяют кинетическую энергию турбулентности и скорость ее диссипации для k- модели, имеют вид (Chen and Kim, 1987):

(2.2.1.6)

(2.2.1.7)

Эмпирические константы C=0,09, k=0,8, =1,15, С1=1,15, C2=1,9, C3=0,25 взяты из работы (Chen and Kim, 1987).


M-SST модель турбулентности

Модель Ментера записывается путем суперпозиции моделей k- и k-, основанной на том, что модели типа k- лучше описывают свойства свободных сдвиговых течений, в модели - имеют преимущество при моделировании пристеночных течений. Плавный переход от k- модели в пристеночной области к k- модели вдали от твердых стенок обеспечивается введением весовой эмпирической функцией F1 .

Вторая важная деталь модели состоит в изменении стандартной связи между k,  и турбулентной вязкостью. Модификация этой связи состоит во введении перехода к формуле Бредшоу в пристеночной области. Согласно предположению Бредшоу напряжение сдвига в пограничном слое пропорционально энергии турбулентных пульсаций.

Уравнения M-SST модели [M-SST]:

Транспортные уравнения для k и 

(2.2.1.8)

(2.2.1.9)


В члене генерации турбулентной энергии введен ограничитель:

(2.2.1.10)

Весовая функция и ее аргумент

, (2.2.1.11)

где положительная часть перекрестных диффузионных членов в уравнении переноса 

. (2.2.1.12)


Выражения для вихревой вязкости с учетом гипотезы Бредшоу:

, (2.2.1.13)

где величина завихренности

. (2.2.1.14)

Функция переключения F2 определяется подобно F1 :

. (2.2.1.15)

Константы в уравнениях переноса записываются путем суперпозиции констант для модели k- модели (Wilcox) и констант стандартной k- модели:

константы:

(2.2.1.16)


Набор констант для пристеночного слоя SST модели:

. (2.2.1.17)

Набор констант для свободных сдвиговых слоев:


. (2.2.1.18)


Другие константы, использующиеся в модели:



2.2.2. Метод пристеночных функций

Вблизи стенок турбулентный поток имеет сложную структуру со значительными градиентами. Во внутренней области, занимающей порядка 20% толщины пограничного слоя, генерируется около 80% всей энергии турбулентности. Для моделирования пристеночных течений часто применяется метод пристеночных функций, который позволяет экономить вычислительные ресурсы и учитывать влияния различных факторов (шероховатости стенки, градиента давления вблизи стенки, вдув-отсос с поверхности тела) за счет введения эмпирической информации. Метод пристеночных функций, связывающий параметры течения с расстоянием до стенки, представляет собой аналитическое решение упрощенных решений турбулентного течения.

В основе упрощения лежат следующие предположения:

1) локальное равновесие энергии турбулентных пульсаций (локальная скорость производства турбулентности уравновешивается скоростью вязкой диссипации);

2) постоянство касательного напряжения вдоль слоя;

3) пренебрежимо малый градиент давления вдоль потока;

4) локальная изотропность диссипирующих вихрей;

5) универсальность профиля скорости (слабая зависимость от числа Рейнольдса).

Следует отметить, что этот метод является в основном эмпирическим. Метод пристеночных функций, основанный на обобщении экспериментальных данных для развитых турбулентных течений, достаточно хорошо описывает пристеночные течения при значительных турбулентных числах Рейнольдса с пренебрежимо малой областью вязкого слоя. В случае не выполнения гипотез, лежащих в основе пристеночных функций, требуется модификация соотношений метода или использование низкорейнольдсовых диссипативных моделей.

2.3. Дискретизация уравнений переноса

Для решения уравнения переноса энергии используется широко известный метод контрольного объема, суть которого заключается в разбиении расчетной области на контрольные объемы (рис. 2.3.1) и интегрировании исходных уравнений сохранения по каждому контрольному объему для получения конечно разностных соотношении. Метод и его свойства достаточно подробно описаны в литературе (Белов и др., 1989), (Ferziger and Peric, 2002).

В пакете «Fire» используются многоблочные структурированные криволинейные сетки.




Рис. 2.3.1. К вычислению потоков через грань контрольного объема.


В расчетном алгоритме применяется подход стыковки отдельных сеток с различной плотностью расчетных узлов вдоль выделенной поверхности раздела, при этом сетки не пересекаются. В этом случае граничный контрольный объем может соприкасаться с несколькими контрольными объемами из соседнего блока. Соответственно граничная грань ячейки разбивается на несколько подобластей. Для каждой области пересечения строится полный поток. Рассчитанные потоки разбиваются на явную и неявную части и заносятся в правую часть и матрицу коэффициентов системы уравнений соответственно. Использование значения величины из соседнего блока с предыдущей итерации для вычисления потока через общую грань несколько снижает скорость сходимости по сравнению с одноблочной техникой. Но многоблочная технология позволяет существенно снизить вычислительные затраты для построения сетки и проведения расчетов тепловых и гидродинамических процессов в многосвязных областях сложной геометрии.


2.4. Граничные условия.


Условия на входе потока

В пакете «Fire» существует возможность задавать на входе либо фиксированное значение скорости, либо массовый расход. При этом если задан расход то компоненты скорости u, v, w в каждом входе рассчитываются по заданному массовому расходу, углу наклона оси входного потока к осям координат и по доли потока, заданное для данного входного «окна». Например, компонента u на входе определяется следующим образом:

(2.4.1)

где G - массовый расход жидкости на входе,  - суммарная плотность, А-площадь сечения входного окна,  - угол между осью OX и осью входного потока.

Кинетическая энергия турбулентности k и скорость диссипации турбулентности  на входе оцениваются по соотношениям вида:

(2.4.2)

где константы koin, oin подбираются эмпирически.

Обычно для развитых турбулентных потоков на входе уровень турбулентных пульсации составляет порядка 1%, поэтому koin=0,01, а oin =C/L, где L - характерный размер энергосодержащих вихрей.

Энтальпия на входе вычисляется по заданным температуре входа и составу газов или жидкости.

Граничные условия для Ym задаются в виде массовых долей компонентов на входе.


Условия на выходе потока

На выходной границе при решении уравнений для u, v, w, k, , h, Ym ставятся условия отсутствия градиентов (так называемые «мягкие условия»):

, (2.4.3)

где n - вектор внешней нормали к расчетной области.


Условие симметрии (скольжение)


На плоскости симметрии ставятся условия равенства нулю: производной по нормали к плоскости симметрии всех скалярных величин и тангенциальной составляющей скорости , нормальной к плоскости составляющая скорости .

, , (2.4.4)

где n- вектор нормали к плоскости симметрии.


Твердая стенка

Нормальную и тангенциальную компоненты скорости на стенках полагаем равными нулю, что моделирует соответственно непротекание и прилипание.

, (2.4.5)

Для определения турбулентных характеристик вблизи стенки используется метод пристеночных функций.

Для массовых долей задаются условия отсутствия диффузионных потоков через стенку (конвективные потоки отсутствуют в силу (2.4.5))

(2.4.6)

где n- вектор нормали к стенки.


Теплоотдача на стенках


Граничные условия для уравнения на энтальпию можно задать двумя способами:

1) Задание постоянной температуры на стенке

(2.4.7)

2) Задание теплового потока и теплоотдачи на стенке

(2.4.8)

где - температура на стенке, - заданная температура, - температура внешней среды, -заданный тепловой поток через стенку, -коэффициент теплоотдачи.


2.5. Методы решения систем линейных алгебраических уравнений


После дискретизации исходных уравнении по времени и пространству получаем системы линейных алгебраических уравнений:



(2.5.1)

Для решения систем линейных алгебраических уравнений в пакете реализовано несколько методов, описание которых можно найти в соответствующей литературе:

полинейный,

переменных направлений,

неполной факторизации Булеева (Булев Н.И., 1989),

метод сопряженных невязок с факторизацией по Булееву (Булев Н.И., 1989)


2.6. Радиационный перенос энергии.


Моделирование процесса радиационного теплообмена является очень сложной и ресурсоемкой задачей, т.к. в отличие от остальных процессов тепломассопереноса каждый элементарный объем среды находится в непосредственном взаимодействии со всеми другими элементарными объемами и решение интегро-дифференциальных уравнений, описывающих это явление, очень трудоемко. Вычислительные затраты на решение уравнения радиационного теплопереноса (УРТ) (ур. 2.6.1) может существенно превосходить затраты на решение всех остальных процессов.

(2.6.1)

Граничные условия на УРТ для случая диффузионного отражения

(2.6.2)

Радиационный источниковый член в уравнении энергии

(2.6.3)

где

где I – интенсивность (Вт/(ср*м2)); Ib – интенсивность абсолютно черного тела, s коэффициент рассеивания, м-1; к – коэффициент поглощения, м-1; b - коэффициент затухания (b = к + s), м-1; W - телесный угол, рад; r – радиус-вектор, м; s – угловое направление, м; Ф – функция рассеивания.

Коэффициенты поглощения газа к вычисляются по модели суммы серых газов (WSGG модели), коэффициенты поглощения рассчитываются через оптическую плотность дыма.

Для задачи радиационного обмена в разработан ряд подходов. Это методы имитационного моделирования, в которых прослеживается движение отдельных лучей и фиксируется изменение их энергии вдоль траектории. К ним относятся различные варианты метода Монте-Карло и методы дискретных направлений.

Ряд других методов основан на представлении исходного интегро-дифференциального уравнения переноса излучения в форме дифференциальных уравнений второго порядка. Наиболее простой из этих методов - это диффузионный метод, который получается интегрированием уравнения переноса излучения по всему телесному углу. Диффузионный метод плохо работает в оптически тонких средах, а метод Монте-Карло достаточно требователен к вычислительным затратам, поэтому здесь представлены два метода дискретных направлений: конечно-объемный метод (КОМ) и дискретно-ординтный метод (ДОМ)

Конечно-объемный метод

КОМ – лучевой метод. Лучевые методы не имеют ограничений по применению, их точность зависит, прежде всего, от дискретизации углового пространства, однако, требуют существенных вычислительных ресурсов. В КОМ переход к разностному аналогу ур. 2.6.1 происходит интегрированием по контрольному объему и угловому пространству (ур. 2.6.1).

(2.6.4)

(2.6.5)

где l - l–ый контрольный телесный угол для КОМ.

Разностный аналог УРТ для КОМ

(2.6.6)

(2.6.7)

(2.6.8)

где i – грань контрольной ячейки, а l – контрольный телесный угол

Дискретно-ординтный метод

В ДОМ переход к разностному аналогу ур. 2.6.1 происходит интегрированием по контрольному объему (ур. 2.6.1).

(2.6.9) где где μ, ,  - направляющие косинусы; l - l–ое угловое направление для ДОМ.

Интегралы по телесному пространству в (6.1) заменяются квадратурами, суммируемыми по ординатным направлениям

(2.6.10)




где на квадратуры накладывается следующее условие:

Разностный аналог УРТ для ДОМ




(2.6.11)

где:

(2.6.12)

(2.6.13)

2.7. Химическое реагирование и горение.

2.7.1. Модели горения газообразного топлива

Расчет химической кинетики горения газового топлива основан на использовании глобальных необратимых реакцией между реагентами (горючими и окислителем) с образованием продуктов реакции.





сохранение массы:



изменение концентрации за счет реагирования:

(2.7.1.1)

сток, (компонент является реагентом):

(2.7.1.2)

источник, (компонент является продуктом реакции):



(2.7.1.3)

где скорость реагирования

изменение энтальпии:



изменение термической энтальпии:



(2.7.1.4)




2.7.2. Модели скоростей реакции.

Кинетическая модель.

Скорость горения i-го реагента в этой модели определяется соотношением Аррениусовского типа:


, (2.7.2.1)

где - молярная концентрация и скорость реакции i-го реагента, - предэкспоненциальный множитель и энергия активации реакции соответственно, эмпирические константы, XOX – молярная концентрация окислителя (кислорода).


Модель "обрыва вихря" (eddy break up model).

Скорость горения в этой модели ограничивается скоростью турбулентного перемешивания топлива и окислителя (Magnussen, B.F., and Hjertager, B.W., 1981):


(2.7.2.2)

где A и В - эмпирические константы равные соответственно 4,0 и 0,5;

Yi, YOX, YPR - концентрации (массовые доли) реагента, окислителя и продуктов сгорания;

SOX, SPR - стехиометрические коэффициенты в реакции горения.


Гибридная модель (кинетика/обрыв вихря)

Для описания реагирования в турбулентных течениях с большим временем перемешивания применяется гибридная модель, использующая оба вышеприведенных механизма определения скорости реакции. По этой модели в качестве результирующей скорости выбирается наименьшая из скоростей (2.7.2.1, 2.7.2.2):


. (2.7.2.3)

Модель комбинированного временного масштаба

Модель "обрыва вихря" дает завышенную скорость горения в пристеночных областях. Эту проблему удается избежать путем добавления времени, определяемого кинетическим механизмом горения, к времени турбулентного перемешивания:



(2.7.2.4)


^ Основные расчетные величины и зависимости, используемые для оценки индивидуального пожарного риска


3.1. Показателем оценки уровня обеспечения пожарной безопасности людей на объектах является вероятность предотвращения воздействия (PВ) ОФП.

Вероятность предотвращения воздействия ОФП определяют для пожароопасной ситуации, при которой место возникновения пожара находится на первом этаже вблизи одного из эвакуационных выходов объекта защиты.

Вероятность предотвращения воздействия ОФП (PВ) на людей в объекте вычисляют по формуле

, (3.1.1)

где QB — расчетная вероятность воздействия ОФП на отдельного человека в год.

Уровень обеспечения безопасности людей при пожарах отвечает требуемому, если

, (3.1.2)

где QBH — допустимая вероятность воздействия ОФП на отдельного человека в год (10-6 год-1).

3.2. Вероятность (QB) вычисляют для людей в каждом здании (помещении) по формуле

, (3.2.1)

где Qп – частота возникновения пожара в здании в течение года, определяется расчетом по методике, изложенной в ГОСТ 12.1.004-91, или на основании статистических данных (Таблица 3.2.1). Если известны данные о количестве людей в объекте защиты, то используется уточненная оценка, если нет – оценка в расчете на одно учреждение. При отсутствии статистической информации допускается принимать Qп = 410-2 для каждого объекта защиты. Кроме этого оценку частотных характеристик возникновения пожара допускается выполнять по статистическим данным, опубликованным в научно-техническом журнале «Пожарная безопасность».


Таблица 3.2.1

Статистические данные о частоте возникновения пожара в зданиях

№ п/п

Наименование здания

Частота возникновения пожара в течение года




В расчете на одно учреждение

Уточненная оценка

1.

Детские дошкольные учреждения (детский сад, ясли, дом ребенка)

7,34*10-3


9,72*10-5

(в расчете на одного ребенка)

2.

Общеобразовательные учреждения (школа, школа-интернат, детский дом, лицей, гимназия, колледж)

1,16*10-2


4,16*10-5

(в расчете на одного учащегося)

3.

Начального профессионального образования (профессиональное техническое училище)

1,98*10-2


4,59*10-5

(в расчете на одного учащегося)

4.

Среднего профессионального образования (среднее специальное учебное заведение)

2,69*10-2


2,94*10-5

(в расчете на одного учащегося)

5.

Высшего профессионального образования (высшее учебное заведение)

1,398*10-1


2,43*10-5

(в расчете на одного учащегося)

6.

Прочие внешкольные и детские учреждения

1,52*10-2


2,38*10-5

(в расчете на одного учащегося)

7.

Детский оздоровительный лагерь, летняя детская дача

1,26*10-3


3,23*10-5

(в расчете на одного отдыхающего)

8.

Больница, госпиталь, клиника, родильный дом, психоневрологический интернат и др. стационары

3,66*10-2


2,358*10-4

(в расчете на одно койко-место)

9.

Санаторий, дом отдыха, профилакторий, дом престарелых и инвалидов

2,99*10-2


1,767*10-4

(в расчете на одно койко-место)

10.

Амбулатория, поликлиника, диспансер, медпункт, консультация

8,88*10-3


5,37*10-5

(в расчете на одно посещение)

11.

Предприятия розничной торговли:

универмаг, промтоварный магазин;

универсам, продовольственный магазин;

магазин смешанных товаров;

аптека, аптечный ларек;

прочие здания торговли


2,03*10-2


1,579*10-3

(в расчете на одного работающего)

12

Предприятия рыночной торговли:

крытый, оптовый рынок (из зданий стационарной постройки)

торговый павильон, киоск, ларек, палатка, контейнер

1,13*10-2


1,678*10-3

(в расчете на одного работающего)

13

Предприятия общественного питания

3,88*10-2


2,063*10-3

(в расчете на одного работающего)

14

Гостиницы, мотели

2,81*10-2


3,255*10-4

(в расчете на одно место)

15

Спортивные сооружения

1,83*10-3




16

Клубные и культурно-зрелищные учреждения

6,90*10-3




17

Библиотеки

1,16*10-3




18

Музеи

1,38*10-2




19

Зрелищные учреждения (театры, цирки)

9,66*10-2

4,03*10-7

(в расчете на одно посещение)



Рпр – вероятность присутствия людей в объекте защиты, определяемая из соотношения Рпр= tфункц/24, где tфункц – время нахождения людей в объекте защиты в часах;

^ PЭ — вероятность эвакуации людей;

PП.З, — вероятность эффективной работы технических решений противопожарной защиты.

3.3. Вероятность эвакуации (PЭ) вычисляют по формуле

(3.3.1)

где РЭ.П — вероятность эвакуации по эвакуационным путям;

PД.В — вероятность эвакуации по наружным эвакуационным лестницам, переходам в смежные секции здания.

3.4. Вероятность (РЭ.П) вычисляют по зависимости

(3.4.1)

где tбл — время от начала пожара до блокирования эвакуационных путей в результате распространения на них ОФП, имеющих предельно допустимые для людей значения, мин;

tр — расчетное время эвакуации людей, мин;

tн.э — интервал времени от возникновения пожара до начала эвакуации людей, мин.

3.5. В объектах защиты с массовым пребыванием людей в одно время (50 и более человек) вероятность эвакуации при пожаре рассчитывается по (3.6):

(3.5.1)

где tск – время существования скоплений людей на участках пути (плотность людского потока на путях эвакуации превышает значение 0,5)

3.6. Расчетное время эвакуации tР определяют при максимально возможной расчетной численности людей в объекте защиты на основе модели движения людей при эвакуации (пункт 3.7);

3.7 Аналитическая модель движения людского потока (ГОСТ 12.1.004-91)

Расчетное время эвакуации людей tр из помещений и зданий устанавливают по расчету времени движения одного или нескольких людских потоков через эвакуационные выходы от наиболее удаленных мест размещения людей.

При расчете весь путь движения людского потока подразделяют на участки (проход, коридор, дверной проем, лестничный марш, тамбур) длиной li и шириной δi. Начальными участками являются проходы между рабочими местами, оборудованием, рядами кресел и т. п.

При определении расчетного времени длину и ширину каждого участка пути эвакуации для проектируемых зданий и сооружений принимают по проекту, а для построенных – по факту. Длину пути по лестничным маршам, а также по пандусам измеряют по длине марша. Длину пути в дверном проеме принимают равной нулю. Проем, расположенный в стене толщиной более 0,7 м, а также тамбур следует считать самостоятельными участками горизонтального пути, имеющими конечную длину li.

Расчетное время эвакуации людей tр следует определять как сумму времени движения людского потока по отдельным участкам пути ti по формуле

tp = t1 + t2 + t3 +,...- + ti, (3.7.1)

где t1 — время движения людского потока на первом (начальном) участке, мин;

t1, t2, t3 ,..., ti — время движения людского потока на каждом из следующих после первого участка пути, мин.

Время движения людского потока по первому участку пути ti, мин, рассчитывают по формуле

(3.7.2)

где l1 — длина первого участка пути, м;

v1— скорость движения людского потока по горизонтальному пути на первом участке, м/мин (определяют по таблице 1 в зависимости от плотности D).

Плотность однородного людского потока на первом участке пути D1 рассчитывают по формуле

(3.7.3)

где ^ N1 – число людей на первом участке, чел;

f – средняя площадь горизонтальной проекции человека, м2/чел.;

δ1 – ширина первого участка пути, м.

Скорость v1 движения людского потока на участках пути, следующих после первого, принимают по таблице 1 в зависимости от интенсивности движения людского потока по каждому из этих участков пути, которую вычисляют для всех участков пути, в том числе и для дверных проемов, по формуле

(3.7.4)

где δi, δi-1— ширина рассматриваемого i-го и предшествующего ему участка пути, м;

qi, qi-1— интенсивности движения людского потока по рассматриваемому i-му и предшествующему участкам пути, м/мин [интенсивность движения людского потока на первом участке пути q =qi-1 определяют по таблице 1 по значению D1, установленному по формуле (3.7.3)].

Если значение qi определяемое по формуле (3.7.4), меньше или равно qmax, то время движения по участку пути ti, мин, равно:

(3.7.5)

при этом значения qmax, м/мин, следует принимать равными:

16,5 — для горизонтальных путей;

19,6 — для дверных проемов;

16,0 — для лестницы вниз;

11,0— для лестницы вверх.

Если значение qi, определенное по формуле (3.7.4), больше qmax то ширину δi данного участка пути следует увеличивать на такое значение, при котором соблюдается условие

qi ≤ qmax (3.7.6)

При невозможности выполнения условия (3.7.6) интенсивность и скорость движения людского потока по участку i определяют по таблице 3.7.1 при значении D = 0,9 и более.

Таблица 3.7.1

Интенсивность и скорость движения людского потока при различной на разных участках путей эвакуации в зависимости от плотности

Плотность потока D, м2/м2

Горизонтальный путь

Дверной проем, интенсивность q, м/мин

Лестница вниз

Лестница вверх

Скорость v, м/мин

Интенсивность q, м/мин

Скорость v, м/мин

Интенсивность q, м/мин

Скорость v, м/мин

Интенсивность q, м/мин

0,01

100

1,0

1,0

100

1,0

60

0,6

0,05

100

5,0

5,0

100

5,0

60

3,0

0,10

80

8,0

8,7

95

9,5

53

5,3

0,20

60

12,0

13,4

68

13,6

40

8,0

0,30

47

14,1

16,5

52

15.6

32

9,6

0,40

40

16,0

18,4

40

16,0

26

10,4

0,50

33

16,5

19,6

31

15,6

22

11,0

0,60

28

16,3

19,05

24,5

14,1

18,5

10,75

0,70

23

16,1

18,5

18

12,6

15

10,5

0,80

19

15,2

17,3

13

10,4

13

10,4

0,90 и более

15

13,5

8,5

8

7,2

11

9,9

Примечание — Интенсивность движения в дверном проеме при плотности потока 0,9 и более, равная 8,5 м/мин, установлена для дверного проема шириной 1,6 м и более, а при дверном проеме меньшей ширины интенсивность движения следует определять по формуле q = 2,5 + 3,75 δ



1 — начало участка i

Рис. 3.7.1— Слияние людских потоков


При слиянии в начале участка i двух и более людских потоков (рис.  3.7.1) интенсивность движения qi, м/мин, рассчитывают по формуле

, (3.7.7)

где qi-1 — интенсивность движения людских потоков, сливающихся в начале участка i, м/мин;

δi-1— ширина участков пути слияния, м;

δi—ширина рассматриваемого участка пути, м.

Если значение qi определенное по формуле (3.7.7), больше qmax то ширину δi, данного участка пути следует увеличивать на такое значение, чтобы соблюдалось условие (3.7.6). В этом случае время движения по участку i определяют по формуле (3.7.5).

3.8. Значение времени начала эвакуации tн.э для объектов защиты без систем оповещения вычисляют по результатам исследования поведения людей при пожарах в объектах защиты конкретного назначения.

При наличии в объекте защиты системы оповещения о пожаре значение н.э принимают равным времени срабатывания системы с учетом ее инерционности. При отсутствии необходимых исходных данных для определения времени начала эвакуации в объекте защиты без систем оповещения величину н.э следует принимать равной 0,5 мин — для этажа пожара и 2 мин — для вышележащих этажей.

Если местом возникнов