Материал: Комплексные проблемы техносферной безопасности. материалы Международной научно-практической конференции. Колодяжный С.А., Небольсин В.А
3.Гареев А.М., Галимова Р.Г. Методические указания по расчету тепло- и влагообеспеченности территории. Часть 2 / А.М. Гареев, Р.Г. Галимова. – Уфа: Изд-во «РИЦ БашГУ», 2012. –108 с.
4.Переведенцев Ю.П. Изменения климатических условий и ресурсов Среднего Поволжья: учебное пособие по региональной климатологии / Ю.П. Переведенцев, М.А. Верещагин, К.М. Шанталинский и др. – Казань: Изд-во «Центр Интех», 2011. – 296 с.
ФГБОУ ВО «Башкирский государственный университет», г. Уфа
R.G. Galimova, K.D. Silantiev
ON SOME FEATURES OF THE TEMPORAL AND SPATIAL ANALYSIS OF THE DANGEROUS METEOROLOGICAL PHENOMENA «STRONG HEAT»
IN THE TERRITORY OF THE REPUBLIC OF BASHKORTOSTAN
The article considers options for identifying such a dangerous meteorological phenomenon as «strong heat» by the criteria of Roshydromet using software products in the long-term series of observations of daily intervals. The analyzed period from 1961 to 2015. Studies on the spatio-temporal analysis of the dangerous meteorological phenomenon under consideration were carried out for the territory of the Republic of Bashkortostan
Key words: dangerous meteorological phenomenon «strong heat», Republic of Bashkortostan
Federal State Budgetary Educational Institution of Higher Education
«Bashkir State University», Ufa
УДК 614.841:69
С.А. Сазонова, С.Д. Николенко
ЧИСЛЕННАЯ РЕАЛИЗАЦИЯ ЗОННОЙ МАТЕМАТИЧЕСКОЙ МОДЕЛИ ПОЖАРА
Рассмотрена трехзонная математическая модель пожара в помещении для первой фазы начальной стадии пожара. В среде MathCAD выполнен расчет изменения координаты границы припотолочного слоя в помещении по уравнениям зонной математической модели пожара в случае распространении пожара по твердым горючим материалам
Ключевые слова: трехзонная математическая модель пожара, помещения, граница припотолочного слоя, пожарная безопасность, MathCAD
Зонные математические модели в основном используются для исследования динамики опасных факторов пожара в начальной стадии пожара [1]. В начальной стадии распределение параметров состояния газовой среды по объему помещения очень неравномерно. В этот период времени пространство внутри помещения можно условно разделить на ряд характерных зон с существенно различающимися температурами и составами газовых сред. В объеме помещения можно выделить три характерные зоны: конвективную колонку над очагом пожара, припотолочный слой нагретого газа и воздушную зону с практически неизменными параметрами состояния, равными своим начальным значениям. Математическая модель пожара, базирующаяся на разбиении пространства на характерные области, получила название трехзонной модели. Схема этой модели показана на рис. 1 [1].
На рис. 1 используются следующие обозначения: ук ─ координата нижней грани припотолочного слоя, отсчитываемая от поверхности горения; yдв ─ высота дверного проема; d ─ эквивалентный диаметр очага горения; 2h ─ высота помещения; Gк ─ поток газа, поступающего в припотолочный слой из конвективной колонки, кг/с; GB ─ поток воздуха, поступающий колонку из зоны III, кг/с; Gг ─ поток газа, вытесняемого из помещения, кг/с; ψ ─ скорость выгорания, кг/с; δ ─ расстояние от пола до поверхности горения, м.
41
Рис. 1. Схема трехзонной модели пожара: I – зона конвективной струи (конвективная колонка); II ─ зона припотолочного нагретого газа; III – зона холодного воздуха;
IV – зона наружного воздуха (наружная атмосфера) (оригинальный)
Произведем расчеты изменения координаты границы припотолочного слоя с использованием зонной математической модели пожара в помещении при распространении пожара по твердым горючим материалам (древесина и оргстекло). При расчете используем встроенную функцию системы MathCAD rkfixed. Данные для расчетов следующие:
1.Коэффициент полноты сгорания для твердых сгораемых материалов η = 0,95.
2.φ = 0,6.
3.ср= 1006 Дж/кг/град.
4.π0 = 1,293 кг/м3.
5.Т0 = 293 К.
В табл. 1 приведены численные значения параметров при горении.
Таблица 1
Численные значения параметров при горении
Параметр и единица его измерения |
Численные значения параметров при горении |
||
древесины |
оргстекла |
||
|
|||
Теплота сгорания, , Дж/кг |
17∙ 106 |
26 ∙ 106 |
|
|
|
||
Удельная скорость выгорания , |
0,0236 |
0,011 |
|
|
|
|
|
Приведем пример выполнения расчетов в среде MathCAD. Введем исходные данные: η: = 0,95 φ: = 0,6
ср: = 1006 Дж/кг/град Т0: = 293 К
l1: = 12м π0: = 1,293 кг/м
ψуд: = 0.0236 кг/м/с l2: = 6м
h2: = 6м
QHP : = 17 ∙ 106 Дж/кг g: = 9,8 м/с
χ: = 0,7
δ: = 0
Вычислим площадь потолка Fпот, Qпож, константы С1 и C2 и y0.
42
Fпот: = l1·l2 Fпот: = 72 н
Qпож: = η·ψуд·Qнр·Fr
Qпож = 8,385 ·106Вт
C1: = 
·Qпож С1=0,122
С2: = 
·
С2 = 0,012
У0: = 1,5
У0 = 7,036 м.
Далее решаем уравнение с разделяющимися переменными (1), используя встроенную функцию системы MathCAD rkfixed.
ук = h2-δ D(τ,ук): = -С1-С2·(ук+у0)
Z: = rkfixed(ук,0,29.9,20,D) n: = 0..20
На рис. 2 и в табл. 2 приведены результаты расчетов.
Рис. 2. Результаты расчетов (оригинальный)
Решение задач в среде MathCAD на основе моделирования динамики опасных факторов пожара в помещениях подробно рассмотрены в работах [2, 3, 4]. Время эвакуации людей из помещения может быть сокращено на вредных и опасных производствах, где производственные процессы сопровождаются выделением вредных горючих веществ. Для таких производств дополнительно рассматриваются задачи экологические, безопасности труда, а так же комплексной безопасности [5] при решении проблемы пожаровзрывобезопасности. Так как пожаровзрывобезопасность необходимо обеспечить в производственных помещениях, то требуется учитывать наличие в них технических трубопроводов, транспортирующих взрывоопасные целевые продукты [6, 7]. Вопросам предотвращения утечек на трубопроводных системах посвящены работы [8, 9]. Такие задачи решаются с привлечением информационных технологий [10] и с обеспечением информационной безопасности [11].
43
Таблица 2
Результаты расчетов Z
Z= |
0 |
1 |
0 |
|
|
1 |
1,49 |
5,59 |
2 |
2,9 |
5,19 |
3 |
4,48 |
4,79 |
4 |
5,9 |
4,41 |
5 |
7,47 |
4,03 |
6 |
8,9 |
3,6 |
7 |
10,46 |
3,29 |
8 |
11,9 |
2,93 |
9 |
13,45 |
2,57 |
10 |
14,9 |
2,23 |
11 |
16,44 |
1,88 |
12 |
17,9 |
1,55 |
13 |
19,43 |
1,22 |
14 |
20,9 |
0,89 |
15 |
22,42 |
0,57 |
Комплексное решение поставленных задач обеспечит требуемый уровень пожарной безопасности и позволит учесть возможные неблагоприятные сценарии развития пожара, что позволит точнее определить необходимое время эвакуации людей при пожаре из рассматриваемых помещений.
Литература
1.Кошмаров Ю.А. Прогнозирование опасных факторов пожара в помещении [Текст]: учеб. пособие / Ю.А. Кошмаров. – М.: 2000. – 118 с.
2.Сазонова С.А. Расчет коэффициента теплопотерь на начальной стадии пожара с применением информационных технологий [Текст] / С.А. Сазонова, С.Д. Николенко // Моделирование систем и процессов. – 2016. – Т. 9. – № 4. – С. 63-68.
3.Сазонова С.А. Численное решение задач в сфере пожарной безопасности [Текст] / С.А. Сазонова, С.Д. Николенко // Моделирование систем и процессов. – 2016. – Т. 9. – № 4.
–С. 68-71.
4.Николенко С.Д. Автоматизация расчетов по интегральной математической модели времени эвакуации людей при пожаре [Текст] / С.Д. Николенко, С.А. Сазонова // Моделирование систем и процессов. – 2017. – Т. 10. – № 1. – С. 43-49.
5.Николенко С.Д. Обеспечение безопасности земляных работ с применением расчетов прикладной механики [Текст] / С.Д. Николенко, С.А. Сазонова // Моделирование систем и процессов. – 2016. – Т. 9. – № 4. – С. 47-51.
6.Квасов И.С. Оценивание параметров трубопроводных систем [Текст]. В сборнике: Информационные технологии и системы. Материалы III Всероссийской научно-технической конференции. / И.С. Квасов, С.А. Сазонова, В.Е. Столяров. – Воронеж, 1999. – С. 112-113.
7.Квасов И.С. Оценивание параметров трубопроводных систем на основе функционального эквивалентирования [Текст]. В книге: Понтрягинские чтения - Х / И.С. Квасов, С.А. Сазонова, В.Е. Столяров. Воронежский государственный университет, 1999. – С. 219.
8.Квасов И.С. Информационные системы технической диагностики трубопроводных сетей [Текст]. В книге: Математическое моделирование в естественных и гуманитарных науках. Тезисы докладов / И.С. Квасов, С.А. Сазонова, В.Е. Столяров. Воронежский государственный университет. – 2000. – С. 105.
9.Квасов И.С. Синтез систем сбора данных для распределительных гидравлических
44
сетей [Текст]. В сборнике: Информационные технологии и системы Материалы III Всероссийской научно-технической конференции / И.С. Квасов, С.А. Сазонова, В.Е. Столяров. – Воронеж, 1999. – С. 113-115.
10.Жидко Е.А. Информационные риски в экологии ХХI века: концепция управления [Текст] / Е.А. Жидко, Л.Г. Попова // Информация и безопасность. – 2010. – Т.13. – №2. – С.
175-184.
11.Жидко Е.А. Информационная безопасность модернизируемой России: постановка задачи [Текст] / Е.А. Жидко, Л.Г. Попова // Информация и безопасность. – 2011. – Т. 14. –
№2. – С. 181-190.
ФГБОУ ВО «Воронежский государственный технический университет»
S.A. Sazonova, S.D. Nikolenko
NUMERICAL IMPLEMENTATION OF THE ZONE MATHEMATICAL FIELD MODEL
A three-zone mathematical model of a fire in a room for the first phase of the initial stage of a fire is considered. In the MathCAD environment, the changes in the coordinate of the boundary of the ceiling layer in the room are calculated according to the equations of the zone mathematical model of the fire in the case of fire spreading over solid combustible materials
Key words: three-zone mathematical model of fire, premises, boundary of the ceiling layer, fire safety, MathCAD
Federal State Budgetary Educational Institution of Higher Education «Voronezh State Technical University»
УДК 614.841:69
С.А. Сазонова, С.Д. Николенко, Е.М. Локтев
ЧИСЛЕННАЯ РЕАЛИЗАЦИЯ ИНТЕГРАЛЬНОЙ МАТЕМАТИЧЕСКОЙ МОДЕЛИ ПРИ АНАЛИЗЕ РАЗВИТИЯ ПОЖАРА
Статья посвящена численной реализации математической модели пожара для помещений, записанной в безразмерном виде. В принятой за основу для вычислений математической модели пожара, усредненные термодинамические параметры состояния в совокупности входят в состав уравнений состояния, материального баланса пожара, кислородного баланса, баланса продуктов горения, баланса инертного газа, энергии и начальных условий. Приведены результаты вычислений. Выполнен анализ полученных результатов
Ключевые слова: математическая интегральная модель пожара, численная реализация, пожарная безопасность, стадии пожара, начальная стадия пожара, развитие пожара
При рассмотрении интегральной математической модели пожара в помещении среднеобъемными параметрами состояния газовой среды в произвольный промежуток времени t. К таким параметрам относятся плотность πm, давление рm, температура Тm, концентрация хi компонентов газовой среды [1].
В соответствии с [1] в состав интегральной математической модели пожара входят: усредненное уравнение состояния среды; дифференциальное уравнение материального баланса пожара в помещении; уравнение кислородного баланса пожара; уравнением баланса продуктов сгорания; уравнение материального баланса для инертного газа; уравнение энергии пожара. Перечисленные уравнения решаются совместно с начальными условиями, учитывающими данные о параметрах состояния среды в помещении перед пожаром.
Совокупность перечисленных уравнений представляет собой математическое описание пожара. В этих уравнениях присутствуют шесть переменных, их число соответствует количественно числу уравнений.
45