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

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

Рубрика Математика
Вид курсовая работа
Язык русский
Дата добавления 21.03.2015
Размер файла 2,0 M

Отправить свою хорошую работу в базу знаний просто. Используйте форму, расположенную ниже

Студенты, аспиранты, молодые ученые, использующие базу знаний в своей учебе и работе, будут вам очень благодарны.

Размещено на http://www.allbest.ru/

Содержание

Введение

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

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

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

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

5. Прогнозирование обстановки на пожаре к моменту прибытия первых подразделений на тушение

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

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

Заключение

Список литературы

Введение

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

Научные методы прогнозирования ОФП основываются на математическом моделировании пожара. Эти методы не только позволяют предсказать развитие пожара, но и восстановить картину уже происшедшего пожара.

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

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

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

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

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

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

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

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

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

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

- ширина =9 м;

- длина = 12 м;

- высота = 3,6 м.

План склада показан на рис.1.1

Рис. 1.1 План склада для хранения

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

Склад для хранения имеет два одинаковых дверных проема. Их ширина равна 0,8 м и высота 1,9 м. При пожаре дверные проемы открыты.

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

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

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

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

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

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

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

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

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

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

Свойства горючей нагрузки.

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

низшая рабочая теплота сгорания Qрн = 17 МДж/кг;

скорость распространения пламени wлс = 0,005 м/с;

удельная скорость выгорания о = 0,008 кг/(м2с);

удельное дымовыделение Dопг = 1096 Нпм2/кг;

удельное потребление кислорода при горении LО2 = -0,968 кг/кг;

выделение окиси углерода LСО=0,163 кг/кг;

выделение двуокиси углерода LСО2=0,812 кг/кг;

Остальные характеристики горючей нагрузки:

суммарная масса горючей нагрузки Мо= 1300 кг;

длина открытой поверхности lпн = 10 м;

ширина открытой поверхности bпн = 5 м;

высота открытой поверхности от уровня пола hпн = 0 м;

Начальные и граничные условия.

Задаемся начальными и граничными условиями:

температура газовой среды помещения равна Tm0 =19 оС;

температура наружного воздуха составляет Ta = -23 оС;

давления в газовой среде помещения и наружном воздухе на уровне пола равны pa=101300 Па.

Выбор сценария развития пожара.

Место возникновения горения расположено в центре площадки, занятой ГМ.

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

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

Согласно исходным данным в базовой системе дифференциальных уравнений следует положить, что

Gпр = 0; Gвыт = 0;

Gов = 0; Qo = 0,

где Gпр и Gвыт - расходы приточного и вытяжного вентиляторов;

Gов - расход газообразного огнетушащего вещества; Qo - тепловой поток, излучаемый системой отопления.

Для пожара при заданных условиях можно принять в уравнении энергии, что

= 0.

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

С учетом сказанного система основных уравнений ИММП имеет вид

;

;

;

;

где V - объем помещения, мз;

сm, Tm, pm - соответственно среднеобъемные плотность, температура и давление;

мm- среднео6ъемная оптическая плотность дыма, Нп/м;

ж = Xm/L - приведенная среднеобъемная концентрация продукта горения;

X02 - среднео6ъемная концентрация кислорода.

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

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

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

На этой схеме использованы следующие обозначения:

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

удв - высота дверного проема;

dэ - эквивалентный диаметр очага горения;

2h - высота помещения;

GK - поток газа, поступающего в припотолочный слой из конвективной колонки, кг-с";

GB- поток воздуха, поступающий в колонку, из зоны III, кг-с-1;

Gr - поток вытесняемого газа из помещения, кг-с';

- Размещено на http://www.allbest.ru/

скорость выгорания, кг-с-1;

- расстояние от пола до поверхности горения, м.

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

3. Расчет динамики опасных факторов пожара в помещении

Для прогнозирования динамики ОФП использована интегральная математическая модель пожара, которую реализует программа INТМОDЕL. В этой программе для численного решения системы дифференциальных уравнений использован метод Рунге-Кутта-Фельберга 4-5 порядка точности с переменным шагом.

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

Полностью исходные данные для расчетов представлены в табл. 3.1. Результаты расчета динамики опасных факторов пожара на складе для хранения представлены в таблицах 3.2 - 3.5.

Таблица 3.1

Исходные данные для расчета динамики опасных факторов пожара в помещении

Атмосфера: Давление, мм рт. ст

Температура, °С

760

-23

Помещение:

Длина, м

Ширина, м

Высота, м

12

9

3,6

Температура, °С

19

Количество проемов, шт

5

Координаты первого проема:

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

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

ширина, м

вскрытие, °С

0

1,9

0,8

19

Координаты второго проема:

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

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

ширина, м

вскрытие, °С

0

1,9

0,8

19

Координаты третьего проема:

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

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

ширина, м

вскрытие, °С

0,8

2,6

2

19

Координаты четвертого проема:

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

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

ширина, м

вскрытие, °С

0,8

2,6

2

300

Координаты пятого проема:

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

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

ширина, м

вскрытие, °С

0,8

2,6

2

250

Горючая нагрузка: Вид горючей нагрузки: крупы (рис, мука, гречиха, пшеница)

Длина; м

Ширина, м

Количество, кг

Выделение тепла, МДж·кг-1

Потребление О2, кг·кг-1

Дымовыделение, Нп·м2·кг-1

Выделение СО, кг·кг-1

Выделение С02, кг·кг-1

Скорость выгорания, кг/м2 ·час

Линейная скорость пламени, мм·с

10

5

1300

17

0,968

1096

0,163

0,812

28,8

5

Таблица 3.2

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

Время, мин

Т-ра, оС

Конц. О2, масс.%

Задымл. Нп/м

Дальн. вид., м

Конц. СО, масс. %

Конц. СО2 масс. %

Конц. ОВ

масс. %

0

19

23

0

15

0

0

77

1

19

22,999

0,001

15

0

0

77

2

20

22,995

0,007

15

0,001

0,004

76,996

3

22

22,981

0,024

15

0,003

0,014

76,987

4

25

22,955

0,057

15

0,007

0,034

76,968

5

30

22,91

0,109

15

0,013

0,066

76,937

6

38

22,842

0,184

12,96

0,023

0,117

76,889

7

47

22,745

0,284

8,39

0,038

0,188

76,821

8

58

22,612

0,411

5,79

0,058

0,287

76,728

9

71

22,439

0,567

4,20

0,083

0,415

76,607

10

85

22,227

0,743

3,20

0,115

0,571

76,458

11

99

21,982

0,938

2,54

0,151

0,753

76,286

12

113

21,705

1,155

2,06

0,192

0,959

76,091

13

126

21,398

1,398

1,70

0,238

1,187

75,874

14

139

21,063

1,678

1,42

0,289

1,438

75,637

15

151

20,701

2,011

1,18

0,343

1,71

75,379

16

163

20,315

2,415

0,99

0,402

2,003

75,101

17

174

19,908

2,918

0,82

0,465

2,315

74,805

18

183

19,502

3,537

0,67

0,528

2,631

74,505

18,30

185

19,384

3,745

0,64

0,547

2,724

74,417

19

189

19,117

4,274

0,56

0,589

2,936

74,215

20

193

18,754

5,133

0,46

0,648

3,229

73,938

21

196

18,414

6,125

0,39

0,705

3,511

73,671

22

198

18,092

7,258

0,33

0,759

3,783

73,412

23

200

17,789

8,538

0,28

0,812

4,047

73,162

24

202

17,503

9,968

0,24

0,864

4,303

72,92

25

203

17,234

11,549

0,21

0,914

4,551

72,684

26

204

16,98

13,279

0,18

0,962

4,792

72,455

27

205

16,741

15,152

0,16

1,009

5,027

72,233

28

206

16,517

17,161

0,14

1,055

5,255

72,017

29

207

16,307

19,295

0,12

1,099

5,475

71,808

30

208

16,11

21,542

0,11

1,142

5,69

71,605

31

208

15,926

23,886

0,1

1,184

5,897

71,408

32

209

15,755

26,31

0,09

1,224

6,098

71,218

33

209

15,596

28,796

0,08

1,263

6,291

71,034

34

209

15,455

31,308

0,08

1,299

6,471

70,864

35

207

15,407

33,461

0,07

1,319

6,571

70,769

36

207

15,364

35,354

0,07

1,338

6,663

70,681

37

207

15,326

37,029

0,06

1,355

6,748

70,601

38

207

15,294

38,515

0,06

1,37

6,825

70,528

39

207

15,266

39,832

0,06

1,384

6,894

70,463

40

207

15,242

40,997

0,06

1,396

6,956

70,404

41

207

15,223

42,025

0,06

1,408

7,012

70,351

42

207

15,207

42,928

0,06

1,418

7,062

70,304

43

207

15,194

43,72

0,05

1,426

7,106

70,262

44

207

15,183

44,41

0,05

1,434

7,145

70,225

45

207

15,175

45,01

0,05

1,441

7,18

70,192

46

207

15,168

45,529

0,05

1,447

7,211

70,162

47

208

15,163

45,977

0,05

1,453

7,238

70,137

48

208

15,159

46,361

0,05

1,458

7,262

70,114

49

208

15,156

46,69

0,05

1,462

7,283

70,094

50

208

15,154

46,971

0,05

1,466

7,302

70,076

51

208

15,152

47,209

0,05

1,469

7,319

70,061

52

211

15,105

47,485

0,05

1,481

7,376

70,006

53

219

15,026

48,278

0,05

1,5

7,472

69,915

54

216

15,08

48,743

0,05

1,493

7,437

69,948

55

213

15,114

48,96

0,05

1,489

7,417

69,967

56

212

15,135

49,029

0,05

1,487

7,406

69,977

57

203

15,206

49,008

0,05

1,474

7,344

70,036

58

187

15,389

48,305

0,05

1,439

7,17

70,201

59

166

15,705

46,422

0,05

1,378

6,864

70,491

60

145

16,087

43,615

0,05

1,303

6,492

70,844

61

125

16,483

40,296

0,06

1,225

6,105

71,211

62

108

16,865

36,807

0,06

1,151

5,733

71,563

63

94

17,218

33,368

0,07

1,082

5,391

71,888

64

82

17,538

30,101

0,08

1,02

5,081

72,181

65

72

17,827

27,065

0,09

0,964

4,803

72,445

66

64

18,088

24,282

0,1

0,914

4,554

72,682

67

57

18,323

21,745

0,11

0,869

4,329

72,895

68

52

18,535

19,442

0,12

0,829

4,127

73,086

69

48

18,729

17,369

0,14

0,791

3,943

73,261

70

44

18,91

15,508

0,15

0,757

3,772

73,423

71

41

19,078

13,842

0,17

0,725

3,613

73,573

72

39

19,235

12,352

0,19

0,696

3,465

73,714

73

37

19,384

11,021

0,22

0,668

3,326

73,846

74

35

19,524

9,833

0,24

0,641

3,195

73,97

75

33

19,657

8,772

0,27

0,616

3,071

74,088

76

32

19,783

7,825

0,3

0,593

2,953

74,199

77

30

19,903

6,981

0,34

0,57

2,842

74,305

78

29

20,017

6,228

0,38

0,549

2,736

74,406

79

28

20,126

5,556

0,43

0,529

2,635

74,502

80

27

20,23

4,957

0,48

0,509

2,538

74,593

81

26

20,33

4,422

0,54

0,491

2,446

74,681

82

26

20,425

3,946

0,6

0,473

2,358

74,764

82,5

25

20,472

3,727

0,64

0,465

2,315

74,805

Таблица 3.3

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

Время,

Плотн. газа,

Изб. давл.,

Высота ПРД,

Пpиток

воздуха

Истечение

газа

СкВыг

мин

кг/м3

Па

м

м3/с

кг/с

м3/с

кг/с

г/с

0

1,21

33,03

1,42

1,881

2,657

2,197

2,657

0

1

1,21

33,1

1,41

1,85

2,614

2,235

2,703

1,60

2

1,21

33,62

1,38

1,772

2,503

2,357

2,843

6,50

3

1,20

34,92

1,34

1,671

2,36

2,565

3,076

14,70

4

1,1853

37,25

1,29

1,573

2,222

2,853

3,382

26,10

5

1,1649

40,7

1,25

1,497

2,114

3,204

3,732

40,9

6

1,1377

45,26

1,21

1,453

2,052

3,597

4,092

58,9

7

1,1048

50,78

1,18

1,444

2,039

4,01

4,43

80,2

8

1,0673

57,03

1,16

1,465

2,069

4,428

4,726

104,9

9

1,0275

63,61

1,16

1,535

2,168

4,786

4,917

129,8

10

0,9877

70,15

1,17

1,653

2,335

5,057

4,995

151,6

11

0,9504

76,23

1,19

1,788

2,526

5,266

5,005

172,7

12

0,9164

81,76

1,20

1,919

2,71

5,443

4,987

193,7

13

0,8858

86,74

1,21

2,036

2,876

5,603

4,963

215

14

0,8583

91,21

1,22

2,136

3,017

5,757

4,941

237

15

0,8334

95,26

1,23

2,22

3,136

5,908

4,924

259,7

16

0,8108

98,94

1,24

2,29

3,235

6,059

4,912

283,4

17

0,7905

102,21

1,25

2,388

3,373

6,137

4,851

305,2

18

0,7752

104,58

1,27

2,539

3,586

6,066

4,702

319,4

18,30

0,7717

105,04

1,28

2,574

3,636

6,043

4,666

322,3

19

0,7653

106,08

1,29

2,658

3,755

5,986

4,581

329,7

20

0,7588

107,08

1,30

2,722

3,845

5,959

4,522

339,7

21

0,7542

107,81

1,31

2,755

3,892

5,96

4,495

349,6

22

0,7506

108,38

1,31

2,774

3,918

5,973

4,483

359,4

23

0,7476

108,86

1,31

2,785

3,934

5,989

4,478

369,1

24

0,7451

109,27

1,31

2,793

3,945

6,007

4,476

378,7

25

0,7429

109,64

1,31

2,799

3,954

6,023

4,475

388,2

26

0,7409

109,95

1,31

2,804

3,961

6,039

4,474

397,6

27

0,7392

110,23

1,31

2,808

3,967

6,052

4,474

406,7

28

0,7378

110,47

1,31

2,812

3,972

6,063

4,473

415,6

29

0,7365

110,67

1,31

2,816

3,977

6,072

4,473

424,3

30

0,7355

110,84

1,31

2,818

3,981

6,08

4,472

432,6

31

0,7347

110,97

1,31

2,821

3,984

6,087

4,472

440,5

32

0,734

111,08

1,31

2,825

3,99

6,088

4,469

447,9

33

0,7335

111,17

1,31

2,816

3,977

6,109

4,481

455,9

34

0,7338

111,84

1,20

2,251

3,179

7,101

5,211

510

35

0,7359

114,31

0,78

0,718

1,015

11,451

8,427

694

36

0,7368

114,28

0,76

0,693

0,979

11,614

8,557

701,4

37

0,7372

114,33

0,74

0,669

0,945

11,782

8,686

709,3

38

0,7374

114,42

0,72

0,645

0,911

11,947

8,81

717,3

39

0,7373

114,54

0,71

0,623

0,881

12,101

8,922

724,9

40

0,7372

114,66

0,69

0,604

0,853

12,243

9,025

731,9

41

0,737

114,78

0,68

0,587

0,829

12,367

9,114

738,1

42

0,7367

114,9

0,67

0,572

0,808

12,476

9,191

743,6

43

0,7365

115

0,66

0,56

0,791

12,569

9,257

748,2

44

0,7363

115,09

0,65

0,549

0,776

12,647

9,312

752,2

45

0,7361

115,17

0,64

0,541

0,764

12,711

9,357

755,5

46

0,7359

115,24

0,64

0,534

0,754

12,764

9,394

758,2

47

0,7358

115,29

0,63

0,528

0,746

12,807

9,423

760,3

48

0,7356

115,34

0,63

0,524

0,74

12,842

9,447

762,1

49

0,7355

115,38

0,63

0,52

0,735

12,869

9,466

763,5

50

0,7355

115,41

0,63

0,517

0,731

12,89

9,48

764,6

51

0,7354

115,43

0,62

0,515

0,728

12,908

9,492

765,5

52

0,7303

108,18

1,83

6,203

8,762

2,252

1,645

247,8

53

0,7195

119,84

0,37

0,242

0,342

15,344

11,039

890,7

54

0,7236

118,43

0,48

0,348

0,491

14,397

10,418

844,5

55

0,7269

117,45

0,54

0,414

0,585

13,821

10,047

815,3

56

0,7294

116,67

0,59

0,477

0,673

13,311

9,709

790,9

57

0,7431

108,95

1,41

3,357

4,742

5,173

3,844

409,9

58

0,7694

104,24

1,48

3,685

5,205

4,454

3,426

331,2

59

0,8054

98,14

1,52

3,846

5,432

3,907

3,146

255

60

0,846

91,43

1,55

3,871

5,467

3,506

2,966

191,5

61

0,8874

84,68

1,56

3,794

5,359

3,223

2,86

143,2

62

0,9271

78,28

1,56

3,655

5,163

3,024

2,804

108

63

0,9635

72,43

1,56

3,485

4,923

2,884

2,779

82,9

64

0,9961

67,21

1,55

3,309

4,674

2,78

2,769

64,7

65

1,0247

62,63

1,54

3,139

4,433

2,7

2,767

51,4

66

1,0497

58,64

1,53

2,983

4,214

2,635

2,766

41,4

67

1,0707

55,31

1,51

2,804

3,961

2,626

2,812

33,7

68

1,0877

52,6

1,49

2,681

3,787

2,595

2,823

27,80

69

1,102

50,32

1,48

2,58

3,645

2,565

2,826

23,20

70

1,1141

48,37

1,47

2,496

3,526

2,535

2,825

19,50

71

1,1246

46,69

1,47

2,426

3,426

2,507

2,82

16,40

72

1,1337

45,24

1,46

2,365

3,341

2,481

2,813

14

73

1,1416

43,97

1,46

2,313

3,268

2,457

2,804

11,90

74

1,1485

42,85

1,45

2,268

3,204

2,434

2,796

10,20

75

1,1547

41,86

1,45

2,229

3,148

2,413

2,786

8,70

76

1,1601

40,99

1,45

2,194

3,099

2,394

2,778

7,50

77

1,165

40,21

1,44

2,163

3,055

2,377

2,769

6,40

78

1,1693

39,51

1,44

2,136

3,017

2,361

2,76

5,50

79

1,1731

38,89

1,44

2,111

2,982

2,346

2,753

4,80

80

1,1766

38,34

1,44

2,089

2,951

2,333

2,745

4,10

81

1,1797

37,84

1,44

2,07

2,923

2,321

2,738

3,60

82

1,1825

37,39

1,43

2,052

2,899

2,31

2,732

3,10

82,5

1,1838

37,18

1,43

2,044

2,887

2,305

2,728

2,90

Таблица 3.4

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

Вpемя,

КонцОВ,

Т-pа,

КонцО2

ПолнСг

УделСкВыг

ВыгМасса

СкВыг,

ПлощГоp,

мин

масс.%

гp.С

масс.%

масс.%

кг/м2/ч

кг

г/с

м2

0

77

19

23

89,805

20,8

0

0

0

1

77

19

22,999

89,805

20,8

0,026

1,60

0,28

2

76,996

20

22,995

89,805

20,8

0,233

6,50

1,13

3

76,987

22

22,981

89,805

20,801

0,817

14,70

2,54

4

76,968

25

22,955

89,805

20,803

1,975

26,10

4,52

5

76,937

30

22,91

89,802

20,808

3,902

40,9

7,07

6

76,889

38

22,842

89,796

20,816

6,807

58,9

10,18

7

76,821

47

22,745

89,781

20,828

10,897

80,2

13,86

8

76,728

58

22,612

89,75

20,848

16,369

104,9

18,11

9

76,607

71

22,439

89,688

20,882

23,38

129,8

22,38

10

76,458

85

22,227

89,581

20,94

31,783

151,6

26,07

11

76,286

99

21,982

89,414

21,029

41,471

172,7

29,56

12

76,091

113

21,705

89,168

21,154

52,424

193,7

32,96

13

75,874

126

21,398

88,823

21,317

64,648

215

36,31

14

75,637

139

21,063

88,358

21,52

78,17

237

39,64

15

75,379

151

20,701

87,751

21,76

93,03

259,7

42,96

16

75,101

163

20,315

86,981

22,038

109,279

283,4

46,29

17

74,805

174

19,908

86,03

22,413

126,951

305,2

49,03

18

74,505

183

19,502

84,939

23,018

145,706

319,4

49,96

18,30

74,417

185

19,384

84,644

23,206

151,469

322,3

50

19

74,215

189

19,117

83,772

23,741

165,177

329,7

50

20

73,938

193

18,754

82,558

24,458

185,241

339,7

50

21

73,671

196

18,414

81,311

25,17

205,893

349,6

50

22

73,412

198

18,092

80,046

25,875

227,135

359,4

50

23

73,162

200

17,789

78,772

26,573

248,961

369,1

50

24

72,92

202

17,503

77,501

27,267

271,367

378,7

50

25

72,684

203

17,234

76,24

27,951

294,344

388,2

50

26

72,455

204

16,98

75

28,626

317,887

397,6

50

27

72,233

205

16,741

73,787

29,285

341,984

406,7

50

28

72,017

206

16,517

72,608

29,927

366,62

415,6

50

29

71,808

207

16,307

71,468

30,548

391,783

424,3

50

30

71,605

208

16,11

70,373

31,146

417,456

432,6

50

31

71,408

208

15,926

69,325

31,718

443,617

440,5

50

32

71,218

209

15,755

68,326

32,251

470,245

447,9

50

33

71,034

209

15,596

67,381

32,824

497,315

455,9

50

34

70,864

209

15,455

66,531

36,719

524,71

510

50

35

70,769

207

15,407

66,24

49,971

552,596

694

50

36

70,681

207

15,364

65,976

50,499

580,696

701,4

50

37

70,601

207

15,326

65,743

51,071

608,93

709,3

50

38

70,528

207

15,294

65,541

51,647

637,308

717,3

50

39

70,463

207

15,266

65,368

52,192

665,827

724,9

50

40

70,404

207

15,242

65,222

52,696

694,479

731,9

50

41

70,351

207

15,223

65,101

53,145

723,254

738,1

50

42

70,304

207

15,207

65,002

53,537

752,137

743,6

50

43

70,262

207

15,194

64,921

53,874

781,115

748,2

50

44

70,225

207

15,183

64,855

54,158

810,174

752,2

50

45

70,192

207

15,175

64,803

54,393

839,3

755,5

50

46

70,162

207

15,168

64,762

54,588

868,482

758,2

50

47

70,137

208

15,163

64,73

54,745

897,711

760,3

50

48

70,114

208

15,159

64,705

54,871

926,977

762,1

50

49

70,094

208

15,156

64,685

54,973

956,273

763,5

50

50

70,076

208

15,154

64,671

55,049

985,589

764,6

50

51

70,061

208

15,152

64,659

55,113

1014,911

765,5

50

52

70,006

211

15,105

64,367

17,845

1045,043

247,8

50

53

69,915

219

15,026

63,864

64,131

1078,523

890,7

50

54

69,948

216

15,08

64,209

60,805

1109,906

844,5

50

55

69,967

213

15,114

64,421

58,698

1140,565

815,3

50

56

69,977

212

15,135

64,551

56,947

1170,787

790,9

50

57

70,036

203

15,206

64,996

29,511

1197,837

409,9

50

58

70,201

187

15,389

66,13

23,849

1220,266

331,2

50

59

70,491

166

15,705

68,033

18,356

1237,931

255

50

60

70,844

145

16,087

70,241

13,786

1251,338

191,5

50

61

71,211

125

16,483

72,427

10,31

1261,375

143,2

50

62

71,563

108

16,865

74,42

7,779

1268,91

108

50

63

71,888

94

17,218

76,163

5,966

1274,642

82,9

50

64

72,181

82

17,538

77,658

4,655

1279,078

64,7

50

65

72,445

72

17,827

78,935

3,699

1282,575

51,4

50

66

72,682

64

18,088

80,027

2,979

1285,374

41,4

50

67

72,895

57

18,323

80,962

2,426

1287,646

33,7

50

68

73,086

52

18,535

81,765

2,002

1289,509

27,80

50

69

73,261

48

18,729

82,469

1,669

1291,053

23,20

50

70

73,423

44

18,91

83,092

1,401

1292,345

19,50

50

71

73,573

41

19,078

83,647

1,184

1293,433

16,40

50

72

73,714

39

19,235

84,145

1,004

1294,353

14

50

73

73,846

37

19,384

84,595

0,856

1295,136

11,90

50

74

73,97

35

19,524

85,002

0,731

1295,803

10,20

50

75

74,088

33

19,657

85,371

0,626

1296,374

8,70

50

76

74,199

32

19,783

85,708

0,538

1296,864

7,50

50

77

74,305

30

19,903

86,016

0,463

1297,284

6,40

50

78

74,406

29

20,017

86,298

0,399

1297,647

5,50

50

79

74,502

28

20,126

86,556

0,344

1297,959

4,80

50

80

74,593

27

20,23

86,794

0,297

1298,229

4,10

50

81

74,681

26

20,33

87,012

0,257

1298,462

3,60

50

82

74,764

26

20,425

87,213

0,223

1298,663

3,10

50

82,5

74,805

25

20,472

87,308

0,207

1298,754

2,10

50

Таблица 3.5

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

Вpемя,

Т-pа,

ТемПов,

КоэфТепОб

ПлТеплПот

ТеплПот,

мин

гp.С

гp.С

Вт/м2/К

Вт/м.2

кВт

0

19

19

0

0

0

1

19

19

2,133

0,3

0,53

2

20

19

3,588

2,50

4,23

3

22

20

5,244

11,20

19,32

4

25

20

6,88

33,2

57,25

5

30

21

8,454

75,8

130,49

6

38

23

9,952

145,5

250,58

7

47

25

11,364

247,4

426,07

8

58

28

12,686

384,2

661,72

9

71

31

13,691

545

938,7

10

85

35

14,136

704,8

1213,91

11

99

39

14,601

873

1503,62

12

113

43

15,073

1044,3

1798,69

13

126

48

15,543

1215,8

2094,02

14

139

52

16,008

1386,1

2387,28

15

151

57

16,468

1554,8

2677,89

16

163

61

16,922

1722

2965,96

17

174

66

17,365

1885

3246,69

18

183

69

17,722

2016,9

3473,76

18,30

185

70

17,795

2043,7

3520,05

19

189

72

17,965

2106,4

3627,97

20

193

73

18,129

2166,7

3731,89

21

196

75

18,248

2210,8

3807,86

22

198

76

18,343

2245,7

3867,88

23

200

76

18,422

2274,9

3918,22

24

202

77

18,49

2300,3

3961,88

25

203

78

18,551

2322,6

4000,33

26

204

78

18,604

2342,3

4034,27

27

205

79

18,651

2359,6

4064,07

28

206

79

18,692

2374,6

4089,95

29

207

80

18,727

2387,5

4112,06

30

208

80

18,756

2398,2

4130,59

31

208

80

18,78

2407

4145,72

32

209

80

18,799

2414

4157,79

33

209

81

18,812

2419,1

4166,57

34

209

80

18,803

2415,6

4160,5

35

207

80

18,743

2393,4

4122,28

36

207

80

18,718

2384,3

4106,63

37

207

79

18,707

2380

4099,24

38

207

79

18,703

2378,7

4096,99

39

207

79

18,704

2379,2

4097,87

40

207

79

18,709

2380,8

4100,56

41

207

79

18,714

2382,9

4104,19

42

207

80

18,721

2385,2

4108,21

43

207

80

18,727

2387,6

4112,23

44

207

80

18,733

2389,8

4116,04

45

207

80

18,738

2391,8

4119,52

46

207

80

18,743

2393,6

4122,6

47

208

80

18,747

2395,1

4125,27

48

208

80

18,751

2396,5

4127,55

49

208

80

18,754

2397,6

4129,46

50

208

80

18,756

2398,4

4130,91

51

208

80

18,758

2399,2

4132,27

52

211

81

18,904

2452,9

4224,86

53

219

85

19,224

2570,9

4428,01

54

216

83

19,099

2525

4348,91

55

213

82

19,002

2488,9

4286,87

56

212

82

18,929

2462,2

4240,78

57

203

78

18,546

2320,6

3996,99

58

187

71

17,863

2068,8

3563,22

59

166

62

17,038

1764,5

3039,03

60

145

54

16,23

1467,5

2527,58

61

125

48

15,517

1206,2

2077,56

62

108

42

14,919

988,4

1702,39

63

94

38

14,432

811,6

1397,84

64

82

34

14,038

669,5

1153,15

65

72

31

13,721

555,5

956,71

66

64

29

13,464

463,6

798,44

67

57

28

12,575

370,9

638,77

68

52

26

11,995

307,1

528,92

69

48

25

11,471

256,9

442,4

70

44

24

10,993

216,6

373,09

71

41

24

10,553

183,9

316,79

72

39

23

10,144

157,1

270,51

73

37

23

9,763

134,7

232,08

74

35

22

9,405

116,1

199,91

75

33

22

9,068

100,3

172,77

76

32

22

8,75

86,9

149,76

77

30

21

8,448

75,6

130,13

78

29

21

8,161

65,8

113,33

79

28

21

7,888

57,4

98,88

80

27

21

7,627

50,2

86,43

81

26

21

7,377

43,9

75,66

82

26

20

7,138

38,5

66,32

82,5

25

20

7,022

36,1

62,12

Примечание:

1. При ф = 53,0 мин площадь ГМ охвачена огнем полностью.

2. При ф = 82,5 мин - полное выгорание горючей нагрузки.

Графики зависимостей Тm (ф), µm (ф), XО2 (ф), XCO2 (ф), XCO (ф), Sпож (ф), Y* (ф), lвид (ф) представлены на рис. п.3.1.-п.3.8.

Рис. 3.1 График зависимости среднеобъемной температуры от времени Тm (ф)

Рис. 3.2 График зависимости среднеобъемной концентрации оксида углерода от времени СО (ф)

Рис. 3.3 График зависимости среднеобъемной концентрации диоксида углерода от времени СО2 (ф)

Рис. 3.4 График зависимости среднеобъемной концентрации кислорода от времени О2 (ф)

Рис. 3.5 График зависимости среднеобъемной оптической плотности дыма от времени µm (ф)

Рис. 3.6 График зависимости дальности видимости от времени lвид (ф)

Рис. 3.7 График изменения плоскости равных давлений во времени Y* (ф)

Рис. 3.8 График изменения площади пожара во времени Sпож. (ф)

Вывод:

1. Максимальная зависимость среднеобъемной температуры от времени Тm (ф)=219 С0.

2. Максимальная зависимость среднеобъемной концентрации оксида углерода от времени СО (ф)= 1,493.

3. Максимальная зависимость среднеобъемной концентрации диоксида углерода от времени СО2 (ф)= 7,472.

4. Максимальная зависимость зависимости среднеобъемной концентрации кислорода от времени О2 (ф)=23.

5. Максимальная зависимость среднеобъемной оптической плотности дыма от времени µm (ф)=49,029 Нп\м.

6. Максимальная зависимость дальности видимости от времени lвид (ф)=15 м.

7. Максимальная изменения плоскости равных давлений во времени Y* (ф)=1,56 м.

8. Максимальная изменения площади пожара во времени Sпож. (ф)=50 м2.

интегральный зонный модель пожар

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

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

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

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

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

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

К опасным факторам пожара, воздействующим на людей и имущество, относятся:

1) пламя и искры;

2) тепловой поток;

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

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

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

6) снижение видимости в дыму.

Критические значения ОФП принимаем по [2,3] (табл. 4.1).

Таблица 4.1

Предельно допустимые значения ОФП

ОФП, обозначение, размерность

ПДЗ

Температура t, єC

70

Парциальная плотность, кг м-1:

- кислорода, с1

- оксида углерода, с2

- диоксида углерода, с2

0,226

0,00116

0,11

Оптическая плотность дыма, м, Непер?м-1

2,38

Таким образом, критическое значение температуры на уровне рабочей зоны равно 70С. Для определения времени достижения температурой этого значения рассчитаем, какова же будет среднеобъемная температура, если на уровне рабочей зоны температура будет критической. Связь между локальными и среднеобъемными значениями ОФП по высоте помещения имеет следующий вид [11]:

(ОФП - ОФПо) = (ОФПm - ОФПо)Z(4.1)

где ОФП - локальное (предельно допустимое) значение ОФП;

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

ОФПm - среднеобъемное значение опасного фактора;

Z - параметр, вычисляемый по формуле:

(4.2)

где H - высота помещения, м;

h - уровень рабочей зоны, м.

Высоту рабочей зоны h определяем по формуле

h = hпл + 1,7(4.3)

где hпл. - высота площадки, на которой находятся люди, над полом помещения, м.

Наибольшей опасности при пожаре подвергаются люди, находящиеся на более высокой отметке [2]. В нашем случае принимаем hпл. = 0. Тогда

H = 0+1,7 = 1,7 м

Значение параметра Z на уровне рабочей зоны будет равно:

Тогда при достижении на уровне рабочей зоны температуры 70С среднеобъемная температура будет равна:

Этого значения среднеобъемная температура достигает, примерно, через 10 минут после начала пожара (табл. 3.2).

Для успешной эвакуации людей дальность видимости при задымлении помещения при пожаре должна быть не меньше расстояния от наиболее удаленного рабочего места до эвакуационного выхода. Дальность видимости на путях эвакуации должна быть не менее 20 м [ 2 ]. Дальность видимости связана с оптической плотностью дыма следующим соотношением [11]:

lпр = 2,38/(4.4)

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

2,38/20 = 0,119 Нп/м

При этом среднеобъемный уровень задымленности будет равен:

Нп/м

По табл. 3.2 получаем = 5 минут.

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

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

кгм3

Для определения времени достижения концентрацией кислорода этого значения строим график зависимости среднеобъемной плотности кислорода от времени пожара (рис. 3.4).

(4.5)

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

Предельная парциальная плотность оксида углерода на путях эвакуации составляет 1,1610-3 кг/м3. При достижении на уровне рабочей зоны парциальной плотностью СО этого значения, среднеобъемная плотность оксида углерода составит:

кгм3

Такого значения среднеобъемная парциальная плотность СО за время расчета не достигает (табл. 3.2-3.3).

Предельное значение парциальной плотности СО2 на уровне рабочей зоны равно 0,11 кг/м3. При этом среднеобъемное значение плотности диоксида углерода будет равно:

кгм3

Такого значения парциальная плотность СО2 за время расчета не достигает (табл. 3.2-3.3).

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

(4.6)

Рис. 4.1 Зависимость парциальной плотности кислорода в помещении от времени пожара

Как видим, на уровне рабочей зоны, среднеобъемная плотность кислорода=0,226 кгм3 не достигает парциальной плотности.

Рис. 4.2 Зависимость парциальной плотности оксида углерода в помещении от времени пожара

Как видим, при достижении на уровне рабочей зоны среднеобъемной плотности оксида углерода=0,0013 кгм3, то как видно на рис. 4.2 время составляет 5 мин.

Как видим, при достижении на уровне рабочей зоны среднеобъемной плотности диоксида углерода=0,12 кгм3 , то как видно на рис. 4.3 время составляет 6 мин.

Рис. 4.3 Зависимость парциальной плотности диоксида углерода в помещении от времени пожара

Вывод:

1. Максимальная зависимость парциальной плотности кислорода в помещении от времени пожара = 26.

2. Максимальная зависимость парциальной плотности оксида углерода в помещении от времени пожара =1,1.

3. Максимальная зависимость парциальной плотности диоксида углерода в помещении от времени пожара =5,45.

5. Прогнозирование обстановки на пожаре к моменту прибытия первых подразделений на тушение

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

Первые подразделения прибывают на пожар через 10 мин после его начала. В это время площадь пожара составляет 26,07 м2, среднеобъемная температура в помещении составляет 85оС, тогда температура на уровне рабочей зоны для личного состава (принимаем 1,7 м) будет составлять (формула 4.1)

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

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

Среднеобъемная оптическая плотность дыма в помещении 10 минуте пожара составляет 0,743 Нп/м.

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

Тогда дальность видимости на уровне рабочей зоны составит

lвид=2,38/0,676=3,5 м.

Среднеобъемное значение парциальной плотности кислорода в помещении на 10 минуте пожара составляет 0,9877 кг/м3 (рис. 3.3).

Содержание кислорода на рабочем уровне составит:

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

Рис. 5.1а Склад для хранения

Рис. 5.1б Разрез 1-1 склада

Среднеобъемное значение парциальной плотности оксида углерода в помещении на 4 минуте пожара определим по формуле

(5.1)

Рис. 5.1в Разрез 2-2 склада

Рис. 5.1г Разрез 3-3 склада

Рис. 5.1 д Разрез 4-4 склада

Тогда

Содержание оксида углерода на рабочем уровне составит

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

Среднеобъемное значение парциальной плотности диоксида углерода в помещении на 10 минуте пожара определим по формуле

(5.2)

Тогда

Содержание диоксида углерода на рабочем уровне составит

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

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

Рассчитать температурное поле в перекрытии через 34 мин после начала нагревания и установить время достижения на рабочей арматуре температуры 209°С.

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

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

Начальная температура перекрытия 19°С, такую же температуру имеет воздух над перекрытием.

Задачу решаем численным методом конечных разностей.

Решение

Принимаем среднюю температуру плиты за весь период нагревания равной 160 °С. Тогда усредненные за весь период нагревания коэффициенты теплопроводности и температуропроводности будут равные л = 1,2 Вт/(м*К), a=6,3*10-7 м2/c.

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

(п. 6.1)

где Тт - среднеобъемная температура; х0 - половина расстояния от очага горения до места выхода газов из помещения; х - координата, отсчитываемая от очага горения по горизонтали; у - координата, отсчитываемая от поверхности пола по вертикали.

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

tm = 165,5°С.

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

=11,63ехр(0,0023-165,5) = 17 Вт/(м2К).

Максимальная расчетная толщина слоя:

Минимальное число слоев:

Разобьем плиту на 5 слоев. При этом толщина слоя будет равна:

Расчетный интервал времени:

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

В начальный момент времени температура во всех слоях равна начальной:

t0,0t1,0t2,0t3,0t4,0t5,019°С.

Температурное поле через 1 = 25 мин, температура греющей среды (определяется аналогично максимальной температуры среды):

tlml = 133,5°С;

коэффициент теплоотдачи с обогреваемой стороны:

= 11,63*ехр(0,0023*133,5) = 15,8 Вт/м2К;

толщина пограничного слоя:

температура в 1-м слое через 1:

т.е. температуры в слоях с 1-го по 25-й через 1 останутся равными начальным.

Температура в 0-м слое через 1:

температура на обогреваемой поверхности

Температурное поле через 2 = 50 мин, температура греющей среды:

tlm,2 = 143,5°С;

коэффициент теплоотдачи с обогреваемой стороны:

=11,63*ехр(0,0023*143,5) =16,2 Вт/м2К;

толщина пограничного слоя:

температура в 1-м слое через 2:

температура во 2-м слое через 2 :

т.е. температуры в слоях со 2-го по 5-й через 2 останутся равными 160 °С;

температура в 0-м слое через 2:

температура обогреваемой поверхности:

Результаты расчета заносим в таблицу 6.1.

Таблица 6.1

Результаты расчетов

tlm

to

t1c

t1

t2

t3

t4

t5

0

0

19

-

-

19

19

19

19

19

19

19

1

25

135,5

15,8

0,075

71

45

19

19

19

19

19

2

50

145,5

16,2

0,074

90

67,5

45

19

19

19

19

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

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

Длина помещения, м - 12.0

Шиpина помещения, м - 9.0

Высота помещения, м - 3.6

Площадь ГН, м2 - 50.00

Гоpючая нагpузка - лен

Результаты расчета:

Вpемя, мин

Т-ра, гр.С

ВысГр, м

ПлГор, м2

Расх, кг/с

0,00

20,00

3,24

0,00

0,00

0,10

20,00

3,24

0,00

0,43

0,20

21,00

3,24

0,01

0,68

0,30

22,00

3,24

0,02

0,90

0,40

23,00

3,24

0,03

1,09

0,50

24,00

3,24

0,05

1,27

0,70

25,00

3,24

0,09

1,61

0,90

27,00

3,12

0,15

1,82

1,10

29,00

2,94

0,22

1,92

1,30

31,00

2,75

0,31

1,97

1,50

34,00

2,55

0,41

1,98

1,70

37,00

2,35

0,52

1,96

1,90

40,00

2,14

0,65

1,91

2,10

44,00

1,94

0,80

1,85

2,30

48,00

1,73

0,96

1,78

2,50

52,00

1,53

1,13

1,71

2,70

57,00

1,33

1,32

1,65

2,90

62,00

1,12

1,52

1,54

3,10

67,00

0,93

1,74

1,38

3,30

73,00

0,74

1,97

1,19

3,50

80,00

0,57

2,22

0,98

3,70

87,00

0,41

2,48

0,75

3,90

95,00

0,26

2,75

0,51

4,10

104,00

0,13

3,04

0,28

4,30

114,00

0,02

3,34

0,04

Рис 7.1 График зависимости температуры от времени

Рис 7.2 График зависимости высоты горения от времени

Рис 7.3 График зависимости площади горения от времени

Рис 7.4 График зависимости расхода топлива от времени

Вывод:

1. Максимальная зависимость температуры от времени=118 С0

2. Максимальная зависимость высоты горения от времени=3,25 м.

3. Максимальная зависимость площади горения от времени=3,4 м2.

4. Максимальная зависимость расхода топлива от времени=2 кг/с.

Заключение

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

Список литературы

1. Пузач С.В. методы расчета тепломассообмена при пожаре в помещении и их применение при решении практических задач пожаровзрывобезопасности. - М.: академия ГПС МЧС России, 2011.

2. Рыжов А.М., Хасанов И.Р., Карпов А.В. и др. Применение полевого метода математического моделирования пожаров в помещениях. Методические рекомендации. - М.: ВНИИПО, 2010.

3. Определение времени эвакуации людей и огнестойкости строительных конструкций с учетом параметров реального пожара: Учебное пособие / Пузач С.В., Казеннов В.М., Горностаев Р.П. - М.: Академия ГПС МЧС России, 2009. 147 л.

4. Астапенко В.М., Кошмаров Ю.А., Молчадский И.С., Шевляков А.Н. Термогазодинамика пожаров в помещениях.- М.: Стройиздат, 1986.

5. Мосалков И.Л., Плюсина Г.Ф., Фролов А.Ю. Огнестойкость строительных конструкций. - М.: Спецтехника, 2012.

6. Кошмаров Ю.А. Прогнозирование опасных факторов пожара в помещении: Учебное пособие. - М.: Академия ГПС МВД России, 2009.

7. Драйздейл Д. Введение в динамику пожаров. - М., Стройиздат, 2008.

8. Яковлев А.И. Расчет огнестойкости строительных конструкций. - М.: Стройиздат, 2008.

9. Кошмаров Ю.А. Теплотехника: учебник для вузов. - М.: ИКЦ «Академкнига», 2009. - 501 с.: ил.

10. Задачник по термодинамике и теплопередаче./ Под ред. Кошмарова Ю.А. Часть 3. - М.: Академия ГПС МВД РФ, 2011.

Размещено на Allbest.ru


Подобные документы

  • Изучение вопросов применения теории множеств, их отношений и свойств и теории графов, а также математических методов конечно-разностных аппроксимаций для описания конструкций РЭА (радиоэлектронной аппаратуры) и моделирования протекающих в них процессов.

    реферат [206,9 K], добавлен 26.09.2010

  • Особенности математических моделей и моделирования технического объекта. Применение численных математических методов в моделировании. Методика их применения в системе MathCAD. Описание решения задачи в Mathcad и Scilab, реализация базовой модели.

    курсовая работа [378,5 K], добавлен 13.01.2016

  • Операторы преобразования переменных, классы, способы построения и особенности структурных моделей систем управления. Линейные и нелинейные модели и характеристики систем управления, модели вход-выход, построение их временных и частотных характеристик.

    учебное пособие [509,3 K], добавлен 23.12.2009

  • Возникновение и развитие теории динамических систем. Развитие методов реконструкции математических моделей динамических систем. Математическое моделирование - один из основных методов научного исследования.

    реферат [35,0 K], добавлен 15.05.2007

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

    презентация [255,0 K], добавлен 28.11.2013

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

    контрольная работа [1,6 M], добавлен 22.12.2014

  • Математическое моделирование динамики биологических видов (популяций) Т. Мальтусом. Параметры и основное уравнение модели "хищник-жертва", ее практическое применение. Качественное исследование элементарной и обобщенной модификаций модели В. Вольтерра.

    курсовая работа [158,1 K], добавлен 22.04.2011

  • Определение вероятности брака проверяемых конструкций. Расчет вероятности того, что из ста новорожденных города N доживет до 50 лет. Расчет математического ожидания и дисперсии. Определение неизвестной постоянной С и построение графика функции р(х).

    курсовая работа [290,7 K], добавлен 27.10.2011

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

    курсовая работа [4,5 M], добавлен 15.02.2010

  • Процесс выбора или построения модели для исследования определенных свойств оригинала в определенных условиях. Стадии процесса моделирования. Математические модели и их виды. Адекватность математических моделей. Рассогласование между оригиналом и моделью.

    контрольная работа [69,9 K], добавлен 09.10.2016

Работы в архивах красиво оформлены согласно требованиям ВУЗов и содержат рисунки, диаграммы, формулы и т.д.
PPT, PPTX и PDF-файлы представлены только в архивах.
Рекомендуем скачать работу.