Уравнение интегральной математической модели пожара называется. Математическое моделирование пожара

Александренко М.В. 1 , Акулова М.В. 2 , Ибрагимов А.М. 3

1 Студент,

Ивановский государственный политехнический университет

МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ ПОЖАРА

Аннотация

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

Ключевые слова: математическая модель, пожар.

Alexandrenko M.V. 1 , Akulova M.V. 2 , Ibragimov A.M. 3

Ivanovo State Polytechnic University

MATHEMATICAL MODELLING OF THE FIRE

Abstract

The article considers types of mathematical models of the fire and their scope. Mathematical modeling allows to predict dynamics of the fire in rooms of buildings of different function and consequently allows to bring research of fire danger of objects to qualitatively new stage of development, to provide transition from comparative methods to expected, considering object service conditions.

Keywords : mathematical model, fire.

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

Моделирование пожара в помещениях основано на представлении пожара как физического явления передачи тепла и массы в соответствующих условиях его развития. Условия развития пожара характеризуются видом пожарной нагрузки и конструктивно-планировочными характеристиками здания (помещения).

По типу математического аппарата различают следующие модели: детерминированные; вероятностные; смешанные (детерминированные – вероятностные); имитационные.

Наиболее эффективным инструментом прогноза и изучения пожаров являются детерминированные математические модели.

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

Приведём краткую характеристику каждой из моделей.

  1. Детерминированные математические модели

Все многообразие детерминированных математических моделей развития пожара в помещениях (внутренние пожары) можно разделить на три группы:

–интегральные (модели первого поколения);

–зонные (модели второго поколения);

–полевые (CFD) (модели третьего поколения).

1.1. Интегральные математические модели

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

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

1.2. Зонные математические модели

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

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

1.3. Полевые математические модели

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

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

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

  1. Вероятностные математические модели

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

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

  1. Смешанные (детерминированные – вероятностные) математические модели

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

  1. Имитационные математические модели

Имитационное моделирование представляет интерес в исследовании сложных систем при априорной неопределенности. В модели может быть задано вероятное протекание пожара, вероятные законы распределения и распространения тепловых потоков, имитируется процесс работы конструкций.

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

Анализ конструктивно-планировочных характеристик помещения;

Определение вида, количества и размещения пожарной нагрузки;

Определение вида возможного пожара и его базовых параметров;

Выбор метода расчета и проведение расчета, оценка вероятностных характеристик пожара;

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

Заключение

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

Литература

  1. Клуб студентов «Технарь». Конспекты по математическим моделям [Электронный курс] URL: http://www.c-stud.ru (дата обращения 10.03.2015)
  2. Расчет необходимого времени эвакуации людей из помещений при пожаре: Рекомендации. – М.: ВНИИПО МВД СССР, 1989. – 22 с.
  3. Методические рекомендации «Применение полевого метода математического моделирования пожара в помещениях.
  4. ГОСТ 12.1.004-91* Пожарная безопасность. Общие требования.
  5. СНиП 21-01-97* Пожарная безопасность зданий и сооружений.

References

  1. Club of students “Technician”. Abstracts on mathematical models of URL: http://www.c-stud.ru (date of the address 10.03.2015)
  2. Calculation of necessary time of evacuation of people from rooms at the fire: Recommendations. – M.: VNIIPO MVD USSR, 1989. – 22 s.
  3. Methodical recommendations “Application of a field method of mathematical modeling of the fire in rooms.
  4. GOST 12.1.004-91 * Fire safety. General requirements.
  5. SNiP 21-01-97 * Fire safety of buildings and constructions.

Расчет динамики опасных факторов пожара в помещении с использованием интегральной математической модели пожара

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

Прогнозирование обстановки на пожаре к моменту прибытия первых

Подразделений на тушение

Расчет огнестойкости ограждающих строительных конструкций

С учетом параметров реального пожара

Расчет динамики опасных факторов пожара в помещении с использованием зонной математической модели пожара

Заключение

Литература

Введение

Для разработки экономически оптимальных и эффективных проти­вопожарных мероприятий необходим научно-обоснованный прогноз ди­намики опасных факторов пожара. Прогнозирование динамики опасных факторов пожара необходимо:

-при создании и совершенствовании систем сигнализации и автома­тических систем пожаротушения;

-при разработке оперативных планов тушения пожаров;

-при оценке фактических пределов огнестойкости;

И для многих других целей.

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

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

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



Исходные данные

Краткая характеристика объекта

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

Размеры цеха в плане:

- ширина =36 м;

- длина = 18 м;

- высота = 6м.

План цеха показан на рис.п.1.1

Рис. п.1.1. План отделочного цеха мебельного комбината

В наружных стенах помещения цеха имеется 3 одинаковых оконных проема, один из которых открытый. Расстояние от пола до нижнего края каждого оконного проема = 0,8 м. Высота оконных проемов = 2,4 м. Ширина каждого оконного проема = 6,0 м. Остекление оконных проемов выполнено из обычного стекла. Остекление разрушается при среднеобъемной температуре газовой среды в помещении, равной 300 0 C.

В противопожарной стене, отделяющей отделочный цех от других помещений, имеется технологический проем шириной 3 м и высотой 3 м. При пожаре этот проем открыт.

Отделочный цех имеет два одинаковых дверных проема, соединяющих цех с наружной средой. Их ширина равна 0,9 м и высота 2 м. При пожаре дверные проемы открыты.

Полы цеха бетонные, с асфальтовым покрытием.

Горючий материал представляет собой деревянные детали мебели, покрытые лаком. Горючий материал расположен на полу. Размер площадки, занятой горючим материалом: длина – 20 м, ширина – 10 м. Количество горючего материала составляет 10 тонн.

Сбор исходных данных

Геометрические характеристики объекта.

Выбирается положение центра ортогональной системы координат в левом нижнем углу помещения на плане (рис. п.1.1). Координатная ось x направлена вдоль длины помещения, ось y - вдоль его ширины, ось z - вертикально вдоль высоты помещения.

Геометрические характеристики:

помещение: длина L =36 м; ширина В = 18 м; высота Н = 6 м.

двери(количество дверей N д o =2): высота h д1,2 = 2,0 м; ширина b д1,2 = 0,9 м; координаты левого нижнего угла двери: у д1 = 10 м; х д1 = 0,0 м; у д2 = 7 м; х д2 = 36,0 м;

открытые окна (количество открытых окон N о o = 1): высота h о o 1 = 2,4 м; ширина b о o 1 = 6,0 м; координаты одного нижнего угла окна: x о o 1 = 3,0 м; у о o 1 = 0 м; z о o 1 = 0,8 м;

закрытые окна (количество закрытых окон N з o =2): высота h з o 1,2 = 2,4 м; ширина b з o 1,2 = 6,0 м; координаты одного нижнего угла окна: x з o 1 = 15 м; y з o 1 = 0,0 м; z T кр = 300 о С; x з o 2 = 27 м; y з o 1 = 0,0 м; z зо1 = 0,8 м; температура разрушения остекления T кр = 300 о С;

технологический проем(количество проемов N п o =1): высота h п1 = 3,0 м; ширина b п1 = 3,0 м; координаты левого нижнего угла проема: у п1 = 18 м; х п1 = 20,0 м.

V - свободный объем помещения, ;

а - коэффициент отражения предметов на путях эвакуации;

Е - начальная освещенность, лк;

Предельная дальность видимости в дыму, м;

Дымообразующая способность горящего материала, ;

L - удельный выход токсичных газов при сгорании 1 кг материала, кг/кг;

X - предельно допустимое содержание токсичного газа в помещении, (; ; );

Удельный расход кислорода, кг/кг.

Если под знаком логарифма получается отрицательное число, то данный ОФП не представляет опасности.

Параметр z вычисляют по формуле:

Высота площадки, на которой находятся люди, над полом помещения, м;

Разность высот пола, равная нулю при горизонтальном его расположении, м.

Следует иметь в виду, что наибольшей опасности при пожаре подвергаются люди, находящиеся на более высокой отметке. Поэтому, например, при определении необходимого времени эвакуации людей из партера зрительного зала с наклонным полом значение h следует находить, ориентируясь на наиболее высоко расположенные ряды кресел. Параметры А и n вычисляют так:

для случая горения жидкости с установившейся скоростью:

,

Удельная массовая скорость выгорания жидкости, ;

для случая горения жидкости с неустановившейся скоростью:

для кругового распространения пожара:

,

V - линейная скорость распространения пламени, м/с;

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

,

b - перпендикулярный к направлению движения пламени размер зоны горения, м.

При отсутствии специальных требований значения а и Е принимаются равными 0,3 и 50 лк соответственно, а значение м.

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

При решении задач с использованием двухзонной модели пожар в здании характеризуется усредненными по массе и объему значениями параметров задымленной зоны:

Т - температура среды в задымленной зоне, К;

Оптическая плотность дыма, Нп/м;

Массовая концентрация i-того токсичного продукта горения в задымленной зоне, кг/кг;

Массовая концентрация кислорода, кг/кг;

Z - высота нижней границы слоя дыма, м.

В свою очередь перечисленные параметры выражаются через основные интегральные параметры задымленной зоны с помощью следующих формул:

, (П6.26)

, (П6.27)

, (П6.29)

где m, - общая масса дыма и соответственно i-го токсичного продукта горения в задымленной зоне, кг;

Масса кислорода в задымленной зоне, кг;

Энтальпия продуктов горения в задымленной зоне, кДж;

S - оптическое количество дыма, ;

Плотность дыма при температуре Т, ;

Объем задымленной зоны, ;

Н, А - высота и площадь помещения, м;

Удельная теплоемкость дыма, .

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

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

, (П6.30)

где t - текущее время, с;

Массовый расход дыма соответственно через конвективную колонку и открытые проемы в помещении, кг/с;

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

, (П6.31)

где , , - тепловая мощность, соответственно, вносимая в задымленную зону конвективной колонкой, удаляемая с дымом через открытые проемы и теряемая в конструкции, кВт;

Введение

В современных условиях разработка экономически оптимальных и эффективных противопожарных мероприятий немыслима без научно-обоснованного прогноза динамики опасных факторов пожара (ОФП).

Прогнозирование ОФП необходимо:

· при создании и совершенствовании систем сигнализации и автоматических систем пожаротушения;

· при разработке оперативных планов тушения (планировании действий боевых подразделений на пожаре);

· при оценке фактических пределов огнестойкости;

· для расчета пожарного риска и многих других целей.

Современные методы прогнозирования ОФП позволяют не только спрогнозировать вероятные пожары, но и смоделировать уже произошедшие пожары для их анализа и оценки действия РТП.

Опасными факторами пожара, воздействующими на людей и материальные ценности (согласно Федеральному закону Российской Федерации от 22 июля 2008 г. №123-ФЗ «Технический регламент о требованиях пожарной безопасности»), являются:

· пламя и искры;

· повышенная температура окружающей среды;

· пониженная концентрация кислорода;

· токсичные продукты горения и термического разложения;

· снижение видимости в дыму;

· тепловой поток.

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

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

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

Методы прогнозирования ОФП различают в зависимости от вида математической модели пожара. Математические модели пожара в помещении условно делятся на три вида: интегральные, зонные и полевые (дифференциальные).

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

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

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

1. Тема и задачи курсовой работы

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

Курсовая работа ставит перед слушателями следующие задачи:

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

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

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

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

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

2. Требования к содержанию и оформлению курсовой работы

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

Необходимый справочный материал дан в приложениях к указаниям и в рекомендуемой литературе.

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

Работа должна быть выполнена аккуратно, чернилами черного цвета или напечатана черным шрифтом на печатных листах формата А4. Текст в пояснительной записке следует писать разборчиво, без сокращений слов (за исключением общепринятых сокращений), на одной стороне листа. Компьютерный вариант работы набирается в текстовом процессоре Word, шрифт Times New Roman с 1-1,5 межстрочным интервалом. Размер шрифта для текста - 12 или 14, для формул - 16, для таблиц - 10, 12 или 14. Размеры полей на листе - 2 см со всех сторон. Абзацный отступ не менее 1 см.

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

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

Страницы курсовой работы должны быть пронумерованы арабскими цифрами. Первой страницей является титульный лист, второй - задание на выполнение курсовой работы, третьей - содержание и т.д. На первой странице курсовой работы номер не ставится. Страницы курсовой работы, кроме титульного листа, и задания на курсовую работу должны быть пронумерованы. Бланк задания на выполнение курсовой работы приведен в приложении 1.

На титульном листе должны быть указаны:

наименование министерства, учебного заведения и кафедры, на которой выполняется курсовая работа;

тема курсовой работы и вариант задания;

Ф.И.О. слушателя, выполнившего курсовую работу;

звание, должность, Ф.И.О. научного руководителя;

город и год выполнения курсовой работы.

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

Если работа удовлетворяет требованиям, предъявляемым к ней, то руководитель допускает ее к защите. Работа, признанная не отвечающим предъявленным требованиям, возвращается обучаемому на доработку.

Защита курсовых работ слушателями факультета заочного обучения может проводиться во время сессии. Результаты защиты оцениваются по четырехбалльной системе: «отлично», «хорошо», «удовлетворительно», «неудовлетворительно». Руководитель проекта проставляет оценку на титульном листе работы, в ведомости, зачетной книжке обучаемого и заверяет подписью. Проставляются только положительные оценки.

При получении неудовлетворительной оценки слушатель обязан повторно выполнить работу по новой теме или переработать прежнюю.

3. Выбор варианта задания и исходные данные

Вариант задания на выполнение курсовой работы определяется по номеру в списке учебной группы (по номеру в журнале группы). Номер варианта указывается на титульном листе курсовой работы. В зависимости от года поступления слушателей на обучение (набор 2010 г., 2011 г. и т.д.) исходные данные для расчетов (температура атмосферного воздуха и внутри помещения, размеры помещения и проемов, параметры горючей нагрузки и т.д.) приведены в таблицах 1-5 (Приложение 2).

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

Дополнительные данные для всех вариантов:

критическая температура для остекления - 300°С;

число проемов - 2 (окна и дверь);

противодымная механическая вентиляция - отсутствует;

автоматическая установка пожаротушения (АУП) - отсутствует;

все остальные не указанные параметры принять по умолчанию.

Сокращения , принятые при изложении курса «Прогнозирование опасных факторов пожара»:

ОФП - опасные факторы пожара;

ПДЗ - предельно-допустимое значение опасного фактора пожара;

ПРД - плоскость равных давлений (нейтральная плоскость);

ГМ - горючий материал.

1. В соответствии с вариантом задания в 1 главе курсовой работы произвести расчет исходных параметров горючей нагрузки в рассматриваемом помещении.

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

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

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

5. По табличным данным построить соответствующие графические зависимости среднеобъемных параметров от времени развития пожара: m (t);

µ m (t); l вид (t); (t); (t); (t); с m (t); Y*(t); S пож (t); G в (t); G г (t); ДP(t).

6. Сделать описание и сравнительные выводы по полученным графикам, объяснить скачки на графиках (при их наличии).

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

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

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

b) по методике определения времени от начала пожара до блокирования эвакуационных путей в результате распространения на них опасных факторов пожара согласно приложению №5 к приказу МЧС России от 10.07.2009 №404 к пункту 33 (Методики определения расчетных величин пожарного риска на производственных объектах).

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

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

10. Провести расчеты параметров ОФП для уровня рабочей зоны (ОФП л) при свободном развитии пожара в момент времени 11 минут, согласно формуле:

(ОФП л - ОФП 0) = (ОФП m - ОФП 0)·Z,

где ОФП л - локальное значение ОФП;

ОФП 0 - начальное значение ОФП;

ОФП m - среднеобъемное значение опасного фактора пожара;- безразмерный параметр, вычисляемый по формуле:

, при H £ 6 м,

где h - высота рабочей зоны, м;

Н - высота помещения, м.

11. Результаты расчетов ОФП для уровня рабочей зоны внести в таблицу в 5 главе курсовой работы.

12. На основании полученных расчетов для времени 11 минут:

а) привести схему газообмена в помещении для времени развития пожара 11 минут при свободном развитии пожара;

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

13. Сделать общий вывод по курсовой работе. Вывод должен включать:

а) краткое описание объекта;

b) анализ ОФП, достигших своего предельно допустимого значения на 11 минуте при свободном развитии пожара;

c) сравнение критического времени наступления ПДЗ по опасным факторам пожара согласно расчетам компьютерной программы INTMODEL и методики определения времени от начала пожара до блокирования эвакуационных путей в результате распространения на них опасных факторов пожара согласно приложению №5 к приказу МЧС России от 10.07.2009 №404;

d) анализ своевременности срабатывания установленных в помещении пожарных извещателей при необходимости предложения по их замене;

e) описание действий персонала объекта при возникновении пожара, исходя из данных, полученных при расчетах;

f) описание действий пожарных подразделений, исходя из положения, что время их прибытия - 10 минута от начала развития пожара;

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

h) вывод о целесообразности и перспективах использования компьютерных программ для расчета динамики ОФП при пожаре.

14. В конце курсовой работы привести список использованной литературы.

5. Образец выполнения курсовой работы

МЧС РОССИИ

Федеральное Государственное бюджетное образовательное

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

«Уральский институт Государственной противопожарной службы

Министерства Российской Федерации по делам гражданской обороны,

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

Кафедра физики и теплообмена

КУРСОВАЯ РАБОТА

Тема: Прогнозирование опасных факторов пожара в складском помещении

Вариант №35

Выполнил:

слушатель учебной группы З-461

старший лейтенант внутренней службы Иванов И.И.

Проверил:

старший преподаватель кафедры

физики и теплообмена, к.п.н., капитан внутренней службы

Субачева А.А.

Екатеринбург

на выполнение курсовой работы

по дисциплине «Прогнозирование опасных факторов пожара»

Слушатель Иванов Иван Иванович

Вариант №35 Курс 4 Группа З-461

Наименование объекта: склад хлопка в тюках

Исходные данные

Блок атмосфера

давление, мм. рт. ст.

температура, 0 С

Блок помещение

высота, м

ширина, м

температура, 0 С

проем 1 - штатный (дверь)

нижний срез, м

∑ ширина, м

верхний срез, м

вскрытие, 0 С

проем 2 - штатный (окна)

∑ ширина, м

нижний срез, м

вскрытие, 0 С

верхний срез, м

вид горючего материала

хлопок в тюках

дымовыделение Нп*м 2 /кг

выделение СО, кг/кг

ширина, м

выделение СО 2 , кг/кг

количество ГН, кг

удельная скорость выгорания, кг/м 2 *с

выделение тепла МДж/кг

скорость распространения пламени, м/с



потребление кислорода кг/кг

Срок сдачи: «____»__________

Слушатель____________________ Руководитель_______________

1. Исходные данные

Помещение пожара расположено в одноэтажном здании. Здание построено из сборных железобетонных конструкций и кирпича. В здании наряду с помещением склада находятся два рабочих кабинета. Оба помещения отделены от склада противопожарной стеной. План объекта приведен на рисунке 1.

(Требуется проставить на схеме размеры помещения и расчетную массу горючей нагрузки согласно своему варианту!)

Рис. 1. План здания

Размеры склада:

длина l 1 = 60 м;

ширина l 2 = 24 м;

высота 2h = 6 м.

В наружных стенах помещения склада имеется 10 одинаковых оконных проемов. Расстояние от пола до нижнего края каждого оконного проема Y H = 1,2 м. Расстояние от пола до верхнего края проема Y B = 2,4 м. Суммарная ширина оконных проемов = 24 м. Остекление оконных проемов выполнено из обычного стекла. Остекление разрушается при среднеобъемной температуре газовой среды в помещении, равной 300°С.

Помещение склада отделено от рабочих кабинетов противопожарными дверьми, ширина и высота которых 3 м. При пожаре эти проемы закрыты. Помещение склада имеет один дверной проем, соединяющий его с наружной средой. Ширина проема равна 3,6 м. Расстояние от пола до верхнего края дверного проема Y в = 3, Y н =0. При пожаре этот дверной проем открыт, т.е. температура вскрытия 20 0 C.

Полы бетонные, с асфальтовым покрытием.

Горючий материал представляет собой хлопок в тюках. Доля площади, занятая горючей нагрузкой (ГН) = 30%.

Площадь пола, занятая ГН, находится по формуле:

=;

где − площадь пола.

Количество горючего материала на 1 Р 0 = 10. Общая масса горючего материала .

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

Свойства ГН характеризуются следующими величинами:

теплота сгорания Q = 16,7 ;

удельная скорость выгорания = 0,0167 ;

скорость распространения пламени по поверхности ГМ ;

дымообразующая способность D = 0,6 ;

потребление кислорода = 1,15 ;

выделение диоксида углерода = 0,578 ;

выделение оксида углерода = 0,0052 .

Механическая вентиляция в помещениях отсутствует. Естественная вентиляция осуществляется через дверные и оконные проемы.

Отопление центральное водяное.

Внешние атмосферные условия:

ветер отсутствует, температура наружного воздуха 20 0 C = 293 К (согласно выбранному варианту);

давление (на уровне Y=h) Р а = 760 мм. рт. ст., т.е. = 101300 Па.

Параметры состояния газовой среды внутри помещения перед пожаром :

Т = 293 К (согласно выбранному варианту);

Р = 101300 Па;

Другие параметры:

критическая температура для остекления − 300 о С;

материал ограждающих конструкций - железобетон и кирпич;

температура воздуха в помещении - 20 о С;

автоматическая система пожаротушения − отсутствует;

противодымная механическая вентиляция − отсутствует.

2. Описание интегральной математической модели свободного развития пожара в складском помещении

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

уравнение материального баланса газовой среды в помещении :

V(dс m /dф) = G B + ш - G r , (1)

где V - объем помещения, м 3 ; с m - среднеобъемная плотность газовой среды кг/м 3 ; ф - время, с; G B и G r - массовые расходы поступающего в помещение воздуха и уходящих из помещения газов, кг/с; ш - массовая скорость выгорания горючей нагрузки, кг/с;

уравнение баланса кислорода :

Vd(p 1)/dф = x 1в G B - x 1 n 1 G r - ш L 1 Ю, (2)

где x 1 - среднеобъемная массовая концентрация кислорода в помещении; х 1в - концентрация кислорода в уходящих газах; n 1 - коэффициент, учитывающий отличие концентрации кислорода в уходящих газах х 1г от среднеобъёмного значения x 1 , n 1 = х 1г /x 1 ; L 1 - скорость потребления кислорода при горении, p 1 - парциальная плотность кислорода в помещении;

уравнение баланса продуктов горения :

Vd(p 2)/dф = ш L 2 Ю - x 2 n 2 G r , (3)

где X i - среднеобъемная концентрация i-гo продукта горения; L i - скорость выделения i-гo продукта горения (СО, СО2); n i - коэффициент, учитывающий отличие концентрации i-гo продукта в уходящих газах x iг от среднеобъёмного значения x i , n i = x iг /х i ; р 2 - парциальная плотность продуктов горения в помещении;

уравнение баланса оптического количества дыма в помещении :

Vd ()/d =Dш - n 4 G r / р m - к c S w , (4)

где - среднеобъемная оптическая плотность дыма; D - дымообразующая способность ГМ; n 4 - коэффициент, учитывающий отличие концентрации дыма в уходящих из помещения нагретых газах от среднеобъемной оптической концентрации дыма, n4= м mг /м m ;

уравнение баланса энергии U:

dU/dф = hQ p н ш + i г ш + С рв Т в G в - С р Т m m G r - Q w , (5)

где P m - среднеобъемное давление в помещении, Па; С р m , Т m - среднеобъемные значения изобарной теплоемкости и температуры в помещении; Q p н - низшая рабочая теплота сгорания ГН, Дж/кг; С рв, Т в - изобарная теплоемкость и температура поступающего воздуха, К; i г - энтальпия газификации продуктов горения ГН, Дж/кг; m - коэффициент, учитывающий отличие температуры Т и изобарной теплоемкости С рг уходящих газов от среднеобъемной температуры Т m и среднеобъемной изобарной теплоемкости С р m ,

m = С рг Т г /С р m Т m ;

Ю - коэффициент полноты сгорания ГН; Q w - тепловой поток в ограждение, Вт.

Среднеобъемная температура Т m связана со среднеобъёмным давлением Р m и плотностью р m уравнением состояния газовой среды в помещении:

P m = с m R m T m . (6)

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

VdP m / dф = ш + G B - G r + G пр - G выт + G ов, (7)

Вышеуказанная система уравнений решается численными методами с помощью компьютерной программы. Примером может служить программа INTMODEL.

. Расчет динамики ОФП с помощью компьютерной программы INTMODEL

Результаты компьютерного моделирования

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

Таблица 1. Динамика развития параметров газовой среды в помещении и координат ПРД

Вpемя, мин Температура t m , 0 С Оптическая плотность дыма µ m , Нп/м Дальность видимости l m , м ,

масс.%,

масс.%, масс.%с m , кг/м 3 Нейтральная плоскость - ПРД Y*, мG в, кг/сG г, кг/сДP, ПаS пож, м 2









Изменение среднеобъемных параметров газовой среды во времени

Рис. 2. Изменение среднеобъемной температуры газовой среды во времени

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

Выводы по графику: На графике температуры можно условно выделить 3 стадии развития пожара. Первая стадия - нарастание температуры (приблизительно до 22 мин.), вторая - квазистационарная стадия (с 23 мин. до 50 мин.), и третья - стадия затухания (с 50 мин. до полного выгорания горючей нагрузки).

Рис. 3. Изменение оптической плотности дыма во времени

Описание графика: В начальной стадии пожара выделение дыма незначительно, полнота сгорания максимальна. В основном дым начинает выделяться после 22 минуты от начала возгорания, а превышение ПДЗ по среднеобъемному значению плотности дыма произойдет примерно на 34 минуте. Начиная с 52 минуты, с переходом в режим затухания, задымление уменьшается.

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

Рис. 4. Изменение дальности видимости в помещении во времени

Описание графика: На протяжении 26 минут развития пожара дальность видимости в горящем помещении остается удовлетворительной. С переходом в режим ПРВ видимость в горящем помещении быстро ухудшается.

Выводы по графику: Дальность видимости связана с оптической плотностью дыма соотношением . То есть дальность видимости обратно пропорциональна оптической плотности дыма, поэтому при увеличении задымления дальность видимости уменьшается и наоборот.


Рис. 5. Изменение среднеобъемной концентрации кислорода во времени

Описание графика: В первые 9 минут развития пожара (начальная стадия) среднеобъемная концентрация кислорода почти не изменяется, т.е. потребление кислорода пламенем низкое, что может быть объяснено малыми размерами очага горения в это время. По мере увеличения площади горения содержание кислорода в помещении снижается. Примерно с 25 минуты от начала горения содержание кислорода стабилизируется на уровне 10-12 масс.% и остается почти неизменным примерно до 49-й минуты пожара. Таким образом, с 25-й по 49-ю минуту в помещении реализуется режим ПРВ, т.е. горение в условиях недостатка кислорода. Начиная с 50-й минуты содержание кислорода увеличивается, что соответствует стадии затухания, при которой поступающий воздух снова постепенно заполняет помещение.

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

Рис. 6. Изменение среднеобъемной концентрации СО во времени развития пожара

Описание графика: сделать описание и выводы по графикам по аналогии с вышеприведенными.

Выводы по графику:

Рис. 7. Изменение среднеобъемной концентрации СО 2 во времени

Описание графика:

Выводы по графику:

Рис. 8. Изменение среднеобъемной плотности газовой среды во времени

Описание графика:

Выводы по графику:

Рис. 9. Изменение положения плоскости равных давлений во времени

Описание графика:

Выводы по графику:

Рис. 10. Изменение притока свежего воздуха в помещение от времени развития пожара

Описание графика:

Выводы по графику:

Рис. 11. Изменение оттока нагретых газов из помещения от времени развития пожара

Описание графика:

Выводы по графику:

Рис. 12. Изменение разности давлений во времени

Описание графика:

Выводы по графику:

Рис. 13. Изменение площади горения при пожаре во времени

Описание графика:

Выводы по графику:

Описание обстановки на пожаре в момент времени 11 минут

Согласно п. 1 ст. 76 ФЗ-123 «Технический регламент о требованиях пожарной безопасности», время прибытия первого подразделения пожарной охраны к месту вызова в городских поселениях и городских округах не должно превышать 10 минут. Таким образом, описание обстановки на пожаре проводится на 11 минуту от начала пожара.

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

− достигается температура 97°С (переходит пороговое значение 70°C);

− дальность видимости практически не изменилась и составляет 64,62 м, т.е. еще не переходит пороговое значение в 20 м;

− парциальная плотность газов составляет:

с= 0,208 кг/м 3 , что меньше предельной парциальной плотности по кислороду;

с= 0,005 кг/м 3 , что меньше предельной парциальной плотности по углекислому газу;

с= 0,4*10 -4 кг/м 3 , что меньше предельной парциальной плотности по угарному газу;

ПРД будет находиться на уровне 0,91 м;

площадь горения составит 24,17 м 2 .

Таким образом, расчеты показали, что на 11 минуту свободного развития пожара, следующие ОФП достигнут своего предельно допустимого значения: среднеобъемная температура газовой среды (на 10 минуте).

. Время достижения пороговых и критических значений ОФП

Согласно ФЗ-123 «Технический регламент о требованиях пожарной безопасности», необходимым временем эвакуации считается минимальное время достижения одним из опасных факторов пожара своего критического значения.

Необходимое время эвакуации из помещения по данным математического моделирования

Таблица 2. Время достижения пороговых значений

Пороговые значения

Время достижения, мин

Предельная температура газовой среды t = 70°C

Критическая дальность видимости 1 кр = 20 м

Предельно допустимая парциальная плотность кислорода с = 0,226 кг/м 3 10


Предельно допустимая парциальная плотность двуокиси углерода (с) пред = (с) пред = 0,11 кг/м 3 не достигается


Предельно допустимая парциальная плотность оксида углерода (с) пред = (с) пред = 1,16*10 -3 кг/м 3 не достигается


Максимальная среднеобъемная температура газовой среды Т m = 237 + 273 = 510 К

Критическая температура для остекления t = 300°C

не достигается

Пороговая температура для тепловых извещателей ИП-101-1А t п opor = 70°C


В данном случае минимальным временем для эвакуации из помещения склада является время достижения предельной температуры газовой среды, равное 10 мин.

Вывод:

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

b ) сделать вывод о своевременности срабатывания пожарных извещателей, установленных в помещении (см. п. 8 таблица 2). В случае неэффективной работы пожарных извещателей предложить им альтернативу (приложение 3).

Определение времени от начала пожара до блокирования
эвакуационных путей опасными факторами пожара

Рассчитаем необходимое время эвакуации для помещения с размерами 60·24·6, пожарной нагрузкой в котором является хлопок в тюках. Начальная температура в помещении 20°С.

Исходные данные:

помещение

свободный объем

безразмерный параметр

;

температура t 0 = 20 0 С;

вид горючего материала - хлопок в тюках - ТГМ, n=3;

теплота сгорания Q = 16,7 ;

удельная скорость выгорания = 0,0167

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

Критическая продолжительность пожара:

t кр = miníý = í746; 772; ý = 746 с.

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

Необходимое время эвакуации людей из складского помещения:

t нв = 0,8*t кр /60 = 0,8*746/60 = 9,94 мин.

Сделать заключение о достаточности / недостаточности времени на эвакуацию по данным расчета.

Вывод: сравнить необходимое время эвакуации, полученное различными методами, и, при необходимости, объяснить различия в результатах.

. Расчет динамики ОФП для уровня рабочей зоны. Анализ обстановки на пожаре на момент времени 11 минут

Уровень рабочей зоны согласно ГОСТ 12.1.004-91 «Пожарная безопасность. Общие требования» принимается равным 1,7 метра.
















На 11 минуте горения газообмен протекает со следующими показателями: приток холодного воздуха составляет 3,26 кг/с, а отток нагретых газов из помещения - 10,051 кг/с.

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

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

Общий вывод по работе

Сделать общий вывод по работе, включающий:

а) краткое описание объекта;

b ) общая характеристика динамики ОФП при свободном развитии пожара;

c ) сравнение критического времени наступления ПДЗ по опасным факторам пожара согласно расчетам компьютерной программы INTMODEL и методики определения времени от начала пожара до блокирования эвакуационных путей в результате распространения на них опасных факторов пожара согласно приложению №5 к приказу МЧС России от 10.07.2009 №404;

d ) анализ срабатывания установленных в помещении пожарных извещателей при необходимости предложения по их замене;

e ) характеристика оперативной обстановки на момент прибытия пожарных подразделений, предложения по проведению безопасной эвакуации людей;

f ) вывод о целесообразности и перспективах использования компьютерных программ для расчета динамики ОФП при пожаре.

Литература

1. Терентьев Д.И. Прогнозирование опасных факторов пожара. Курс лекций / Д.И. Терентьев, А.А. Субачева, Н.А. Третьякова, Н.М. Барбин // ФГБОУ ВПО «Уральский институт ГПС МЧС России». - Екатеринбург, 2012. - 182 с.

2. Кошмаров Ю.А. Прогнозирование ОФП в помещении: Учебное пособие / Ю.А. Кошмаров/ - М.: Академия ГПС МВД России, 2000. -118 с.

Федеральный закон Российской Федерации от 22 июля 2008 г. №123-ФЗ «Технический регламент о требованиях пожарной безопасности».

Приказ МЧС РФ от 10.07.2009 №404 (с изменениями от 14 декабря 2010 г.) «Об утверждении методики определения расчетных величин пожарного риска на производственных объектах». - Пожаровзрывобезопасность. - №8. - 2009. - Стр. 7-12.

Приказ МЧС РФ от 30.06.2009 №382 (с изменениями от 11 апреля 2011) «Об утверждении методики определения расчетных величин пожарного риска в зданиях, сооружениях и строениях различных классов функциональной пожарной опасности». - Пожарная безопасность №3. - 2009. - Стр. 7-13.

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

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

Методы прогнозирования ОФП различают в зависимости от вида математической модели пожара. Математические модели пожара в помещении условно делятся на три вида: интегральные, зонные и полевые (дифференциальные).

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

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

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

Основное преимущество интегральной модели: быстрый и низкотрудоемкий инженерный расчет динамики опасных факторов пожара.

Основные недостатки:

Область корректного применения интегральной модели (по объемам и геометрии помещений, расположению горючего материала и т.д.) является нерешенной проблемой;

Необходимость использования дополнительной экспериментальной информации или моделей более высокого уровня (зонных или полевых) для получения распределения параметров тепломассообмена по объему помещения;

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

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

Основные преимущества:

Быстрый и низкотрудоемкий инженерный расчет динамики опасных факторов пожара;

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

Основные недостатки:

Область корректного применения зонной модели (по объемам и геометрии помещений, расположению горючего материала и т.д.) является нерешенной проблемой;

Необходимость использования дополнительной экспериментальной информации или модели более высокого уровня (полевой) для получения распределения параметров тепломассообмена по объемам зон помещения;

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

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

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

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

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


Характеристика объекта