close

Вход

Забыли?

вход по аккаунту

?

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

код для вставкиСкачать
УДК 004.021
МОДЕЛИРОВАНИЕ ПРОЦЕССОВ РАЗВИТИЯ ГОРЕНИЯ ПОЖАРНОЙ НАГРУЗКИ
С ПОМОЩЬЮ КОНЕЧНЫХ ЦЕПЕЙ МАРКОВА
А.И. Подрезова, В.А. Ловчиков, Ю.Д. Моторыгин, И.А. Пешков
В статье рассматриваются особенности различных видов классических математических моделей для расчета
процесса развития горения. Отмечены их достоинства и недостатки. Показано, что для исследования процесса зажигания можно использовать стохастические или вероятностные методы описания процессов. Приведен пример, иллюстрирующий возможности использования конечных Марковских цепей для анализа процесса возникновения горения
Ключевые слова: моделирование, конечные цепи Маркова, динамика развития пожара
Современная система оценки пожарной опасности включает в себя различные методы стандартных испытаний для определения показателей пожарной опасности [1, 2] и методы исследования динамики горения при полномасштабных испытаниях.
В последние годы получило развитие математическое моделирование процессов горения.
Наибольшее применение находят методы
стандартных испытаний. Получить полную картину
пожарной опасности таким способом не представляется возможным, так как испытания проводятся для
определения одного показателя пожарной опасности, при этом используются строго фиксированные
значения термических воздействий на материалы.
Другие пожароопасные свойства оцениваются уже
при иных условиях.
Полномасштабные испытания являются дорогостоящими и трудоемкими. Количество проведенных, на сегодняшний день, натурных испытаний
единично. Провести натурные испытания даже основных пожароопасных ситуаций с применением
различных видов пожарной нагрузки также не представляется возможным.
Интегральная модель пожара является наиболее простой и позволяет получить информацию, или
сделать прогноз о средних значениях параметров
состояния среды в помещении для любого момента
времени развития пожара. При этом для того, чтобы
сопоставлять среднеобъемные характеристики среды их предельным значениям в зоне пожара используются формулы, полученные на основе экспериментальных данных пространственного распределения температур, концентраций продуктов горения,
оптической плотности дыма и так далее. Такой подход не всегда корректен, а в ряде случаев, в частноПодрезова Анна Игоревна – СПбУ ГПС МЧС России,
соискатель, тел. (473) 264-56-86, e-mail: cosaug@mail.ru
Ловчиков Владимир Александрович – СПбУ ГПС МЧС
России, д-р хим. наук, профессор, тел. (812) 388-86-28
Моторыгин Юрий Дмитриевич – СПбУ ГПС МЧС
России, канд. техн. наук, доцент, тел. (812) 388-86-28,
e-mail udm@gema-peter.ru
Пешков Игорь Александрович - Главное управление МЧС
России по Белгородской области, начальник управления
материально-технического обеспечения, канд. техн. наук,
тел. (473) 294-03-40
176
сти при неоднородности газовой среды, может давать большую погрешность. Интегральная модель
пожара так же не позволяет описывать процессы горения, не достигшие стадии объемного распространения пламени или локальные пожары не получившие развития. Кроме того, в реальных условиях существенное влияние на развитие горения может оказывать состав и распределение пожарной нагрузки,
конвективные потоки и особенности конструкций
помещения.
Зонная модель позволяет получить информацию о размерах характерных пространственных зон,
возникающи при пожаре в помещении и средних
параметров состояния среды. В качестве характерных пространственных зон обычно выделяются
припотолочная область пространства, область восходящего над очагом горения потока нагретых газов
и область не задымленной холодной части пространства в начальный период пожара. При разработке таких моделей приходится делать большое
количество допущений и упрощении основанных на
начальных предположениях об изменениях границ
выбранных зон. Модель позволяет исследовать пожары в помещениях с ярко выраженными границами конвективных потоков и с достаточно четкой
однородностью среды, что оказывается часто недостаточным для решения задач обеспечения пожарной
безопасности в зданиях и сооружениях.
В полевых моделях вместо одной или нескольких зон пространство разбивается на большое
количество незначительных по величине объемов.
Для каждого из этих объемов с помощью численных
методов решается система уравнений в частных
производных выражающих принципы локального
сохранения энергетических и массовых характеристик процесса. Однако система уравнений, описывающая изменение во времени указанных параметров газовой среды внутри помещения чрезвычайно
громоздка и под множеством цифр теряется сама
физика происходящих процессов. Малейшая ошибка на каком-либо шаге вычислений может привести
к большим погрешностям в конечном результате.
Это не позволяет применять данное модельное описание достаточно широко.
Этим традиционно применяемым в пожарном
деле детерминированным методам моделирования,
основанным на уравнениях состояния, законах со-
хранения и сложных системах дифференциальных
уравнений описывающих теплообмен при пожаре
существует альтернативный подход – стохастические или вероятностные методы описания процессов. На сегодня они успешно применяются в самых
разнообразных областях науки и техники. Среди
них простотой и ясностью физического смысла выделяются модели, основанные на теории конечных
цепей Маркова /4, 5, 6/. Конечной цепью Маркова
называется процесс, который переходит из состояния в состояние с определенной вероятностью, так
называемой вероятностью перехода. Число состояний конечно, а значение переходной вероятности
полностью определяется тем, в каком состоянии находится процесс, то есть она является условной. Вероятности перехода образуют стохастическую матрицу Р номер строки которой указывает из какого
состояния происходит переход, а номер столбца в
какое состояние попадает процесс в результате перехода. Все возможные пути процесса описываются
степенями матрицы переходных вероятностей – Pn
[4]. Для переходной матрицы Р существует собственный вектор α, такой что
αР = α,
где α = (α1 α2…..αn)
Здесь n – число состояний моделируемого процесса. Собственный вектор - строка α содержит такое же количество компонент, как и строка переходных вероятностей в матрице Р, а их физический
смысл – среднее время нахождения процесса в состоянии n.
Параметры Марковской модели могут быть
определены экспериментально или с помощью каких либо методов оценки переходных вероятностей.
Для этого могут применяться экспертные методы
или расчеты с помощью традиционных моделей.
В данной работе предлагается одна из возможных экспериментальных методик оценки параметров процесса горения.
Поскольку Марковская цепь событий представляет собой процесс перехода из состояния в состояние, то при ее разработке первой задачей является отождествление модельного состояния с состоянием процесса горения пожарной нагрузки. Так
как горение в большой степени зависит от способа
распределения горючего материала в пространстве,
что по умолчанию предполагается во всех способах
оценки пожарной опасности, то в качестве модельного объекта был выбран пустотелый куб. Используя кубы одинакового размера можно имитировать
простейшее структурированное распределение горючего в пространстве. С технической точки зрения
достоинством кубического модуля является постоянство удельной пожарной нагрузки для любого
размера стороны куба. Это позволяет сохранять постоянным значение важного параметра при моделирования процессов горения.
Например, при исследовании простейшего
случая – горения двух соприкасающихся модулей,
следует различать время горения первого модуля τ1
(рис.1.), модуля который поджигается извне, время
горения первого и второго модуля одновременно τ2
и время горения только второго модуля τ3.
τ1
τ2
τ3
1
2
Рис. 1. Временные характеристики процесса горения
соприкасающихся модулей
Здесь: τ1- время от начала горения модуля №1
до воспламенения модуля №2; τ2 -время совместного горения модуля №1 и модуля №2; τ3 -время горения модуля №2.
Для такого процесса существует матрица переходных вероятностей P:
P ⎞
⎛P
P = ⎜ 11 12 ⎟
⎜P
P ⎟
⎝ 21 22 ⎠
и собственный вектор матрицы α:
α=(α1 α2),
такой что:
αΡ=α.
Здесь α1 вероятность нахождения процесса в
первом модуле или состоянии, α2 во втором состоянии.
Для рассматриваемого процесса справедливо:
α 1Ρ11 + α1 Ρ12 +α 2Ρ21 + α2 Ρ22 = 1,
1
1
τ1 τ2 τ2 τ3
2
2
Слагаемые α 1Ρ11 ~ τ1, α2 Ρ22 ~ τ3 и поскольку из
свойств Марковской цепи известно, что
1
1
α1 Ρ12 = α 2Ρ21 , то α1 Ρ12 ~ τ2 , а α 2Ρ21 ~ τ2 .
2
2
Простые алгебраические преобразования приводят к следующему выражению:
1
τ1 + τ 2
α1=
2
,
τ1 + τ 2 + τ 3
а для вероятности нахождения процесса во
втором модуле справедливо:
1
α2=
2
τ2 + τ3
τ1 + τ 2 + τ 3
.
Поскольку
α1 Ρ11 =
τ1
τ1 + τ 2 + τ 3
то
Ρ11 =
τ1
1
τ1 + τ 2
2
.
,
Аналогично для других вероятностей:
1
2
Ρ12 =
τ2
1
;
τ1 + τ 2
2
1
Ρ21 =
2
1
2
Ρ22 =
τ2
τ 2 +τ3
τ3
1
;
.
τ 2 +τ3
2
Следующей задачей является выбор оптимального размера модуля пожарной нагрузки. Для этого
необходимо экспериментально исследовать зависимость времени горения от изменения размера куба.
В таблице приведены соответствующие экспериментальные данные, показывающие зависимость от размера ребра модельного куба а времени
горения t. В таблице также приведены данные о
массе модельного кубика M, площади его основания
S, массовой скорости выгорания U и удельной массовой скорости выгорания V.
На рис.2 представлены табличные данные
зависимости
удельной
массовой
скорости
выгорания от размера ребра модельного модуля
пожарной нагрузки.
0,00300
V г/(см.кв.)
0,00250
0,00200
Ряд1
0,00150
0,00100
0,00050
0,00000
0
5
10
15
20
25
а см
Рис. 2. Зависимость удельной массовой скорости
выгорания от размера ребра модельного модуля
Из рисунка видно, существуют три области:
- для малых размеров характерно резкое падение удельной массовой скорости выгорания,
- для размеров следующей области наблюдаются практически стабильные значения скорости,
- при дальнейшем увеличении размера ребра
удельная массовая скорость медленно падает.
Это позволило выбрать в качестве оптимального модельного куба модуль с размером ребра 7,5
см.
Рисунок 3 демонстрирует эксперимент горения
двух совмещенных модулей. Два модуля с длиной
грани 7,5 см устанавливаются вплотную друг к другу на теплоизоляционную платформу в металлическом вытяжном шкафу. На модуль №1 на высоте две
№ а
Mг
tс
S см. кв.
U г/с
1
2
3
4
2
2
2
2
0,56
0,55
0,59
0,56
61
52
71
57
4
4
4
4
0,00918
0,01058
0,00831
0,00982
V г/(с·см.
кв)
0,00230
0,00264
0,00208
0,00246
5
2
0,54
50
4
0,01080
0,00270
6
7
8
4
4
4
2,23
2,25
2,22
109
102
109
16
16
16
0,02046
0,02206
0,02037
0,00128
0,00138
0,00127
9
1
0
1
1
1
2
1
3
1
4
1
5
1
6
1
7
1
8
1
9
2
0
2
1
2
2
2
3
2
4
2
5
2
6
2
7
2
8
2
9
4
2,23
104
16
0,02144
0,00134
6
4,9
100
36
0,04900
0,00136
6
4,85
100
36
0,04850
0,00135
6
4,73
98
36
0,04827
0,00134
6
4,8
90
36
0,05333
0,00148
9
10,26
102
81
0,10059
0,00124
9
10,27
103
81
0,09971
0,00123
9
10,59
110
81
0,09627
0,00119
9
1
2
1
2
1
2
1
5
1
5
1
5
1
5
1
5
2
0
2
0
2
0
2
0
10,39
97
81
0,10711
0,00132
17,76
105
144
0,16914
0,00117
17,87
110
144
0,16245
0,00113
17,92
107
144
0,16748
0,00116
27,92
110
225
0,25382
0,00113
27,9
117
225
0,23846
0,00106
27,63
125
225
0,22104
0,00098
27,85
118
225
0,23602
0,00105
27,61
113
225
0,24434
0,00109
48,62
165
400
0,29467
0,00074
48,52
149
400
0,32564
0,00081
48,65
157
400
0,30987
0,00077
48,61
162
400
0,30006
0,00075
третьих длины грани с помощью шприца наносится
в виде пятна инициатор горения - автомобильное
дизельное топливо, объемом 40 мкл. Инициатор горения поджигается при помощи открытого пламени.
Наблюдается потемнение верхней грани модуля и
последующее его горение. Развитие горения происходит во всех направлениях от места поджога. Через
определенное время пламя переходит с модуля №1
на модуль №2, начинается совместное горение обоих модулей (рис.1.). Далее горит только модуль №2.
Прекращение горения отмечается по обрыву языка
пламени.
Таким образом в данной работе:
Предложены модели развития горения на
основе конечных цепей Маркова.
2. Показано, что для этой цели могут быть
использованы структурные модули в виде
куба.
3. Установлено, что их удельная массовая
скорость сложным образом зависит от
размера ребра куба.
4. Проведена оценка значений марковских
параметров модельного процесса горения
пожарной нагрузки составленной из двух
соприкасающихся модулей
1.
Литература
Рис. 3. Горение двух совмещенных модулей
По полученным временам τ1, τ2, τ3 были рассчитаны ориентировочные значения переходных вероятностей. Их усредненные значения приведены
ниже
⎛ 0.65 0.35 ⎞
⎟,
⎝ 0.43 0.57 ⎠
Ρ =⎜
α =(0.56 0.44).
1. ГОСТ 12.1.044–89 Пожаровзрывоопасность веществ и материалов. Номенклатура показателей и методы
их определения.
2. ГОСТ 12.1.004-91 Пожарная безопасность. Общие требования. М.издательство стандартов 1992.-78с.
3. Кошмаров Ю. А. Прогнозирование опасных факторов пожара в помещении. Учебное пособие – М. Академия ГПС МВД РФ, 2000. -118с.
4. Кемени Д., Снелл Д. Конечные цепи Маркова М.
Наука 1970. – 271 с.
5. Баруга-Рид. А.Т Элементы теории Марковских
процессов их приложение М. 1969.-568с.
6. Феллер В. Введение в теорию вероятностей ее
приложению. М. -1984.-528с.
Санкт-Петербургский университет Государственной противопожарной службы МЧС России
Главное управление МЧС России по Белгородской области
MODELLING OF DEVELOPMENTS OF BURNING OF FIRE LOADING
BY THE MARKOV CHAINS WITH A FINITE STATE SPACE
A.I. Podrezova, Y.D. Motorygin, V.A. Lovchikov, I.A. Pechkov
In article features of various kinds of classical mathematical models for calculation of process of ignition are considered.
Their merits and demerits are noted. In work it is shown that for research of process of ignition it is possible to use stochastic or likelihood methods of the description of processes. The example illustrating possibilities of use of final Markov chains for the analysis of
process of occurrence of burning is resulted
Key words: environment, research of ignition process, finite Markov chains
Документ
Категория
Без категории
Просмотров
44
Размер файла
391 Кб
Теги
процессов, моделирование, конечный, пожарной, горение, помощь, маркова, развития, нагрузки, цепей
1/--страниц
Пожаловаться на содержимое документа