МИНИСТЕРСТВО РОССИЙСКОЙ ФЕДЕРАЦИИ ПО ДЕЛАМ
ГРАЖДАНСКОЙ ОБОРОНЫ, ЧРЕЗВЫЧАЙНЫМ СИТУАЦИЯМ
И ЛИКВИДАЦИИ ПОСЛЕДСТВИЙ СТИХИЙНЫХ БЕДСТВИЙ

Федеральное государственное учреждение
«Всероссийский ордена "Знак Почета"
научно-исследовательский институт
противопожарной обороны»

МЕТОДИКА ОБЕСПЕЧЕНИЯ ПОЖАРНОЙ БЕЗОПАСНОСТИ ПЕРЕВОЗКИ САМОВОЗГОРАЮЩИХСЯ ГРУЗОВ

МОСКВА 2006

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

Методика предназначена для региональных подразделений ГПС и специализированных организаций.

Разработана сотрудниками ФГУ ВНИИПО МЧС России: д-ром техн. наук, проф. И.А. Болодьяном; д-ром техн. наук, проф. Ю.Н. Шебеко; д-ром техн. наук, проф. В.И. Горшковым; канд. техн. наук И.А. Корольченко; А.В. Казаковым; Д.Н. Соколовым.

Утверждена ФГУ ВНИИПО МЧС России 29.06.2006 г.

Согласована УГПН МЧС России 21.07.2006 г.

Содержание

I. ОБЩИЕ ПОЛОЖЕНИЯ

II. ОПРЕДЕЛЕНИЕ КИНЕТИЧЕСКИХ ПАРАМЕТРОВ ПРОЦЕССА ТЕРМООКИСЛЕНИЯ МАТЕРИАЛОВ ПО ЭКСПЕРИМЕНТАЛЬНЫМ ДАННЫМ

2.1. Аппаратура

2.2. Подготовка и проведение испытаний

2.3. Расчет параметров кинетики термоокисления материалов

III. РАСЧЕТ УСЛОВИЙ САМОВОЗГОРАНИЯ ГРУЗОВ

3.1. Расчет параметра Франк-Каменецкого δ0 для грузового пространства транспортных средств

3.2. Расчет критической температуры

3.3. Расчет времени индукции

IV. ОПРЕДЕЛЕНИЕ КЛАССА ОПАСНОСТИ САМОВОЗГОРАЮЩИХСЯ ГРУЗОВ И РАЗРАБОТКА ПРОФИЛАКТИЧЕСКИХ МЕРОПРИЯТИЙ

4.1. Определение класса опасности грузов

4.2. Определение безопасных размеров компактной укладки материала

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

ПРИЛОЖЕНИЕ 1 Пример расчета кинетических параметров

ПРИЛОЖЕНИЕ 2 Параметры кинетики процесса термоокисления некоторых материалов

ПРИЛОЖЕНИЕ 3 Пример расчета δ0 грузового пространства

ПРИЛОЖЕНИЕ 4 Значения δ0 для основных типоразмеров грузовых пространств

ПРИЛОЖЕНИЕ 5 Пример расчета критической температуры

ПРИЛОЖЕНИЕ 6 Пример расчета времени индукции

ПРИЛОЖЕНИЕ 7 Пример определения класса опасности самовозгорающегося груза

ПРИЛОЖЕНИЕ 8 Пример определения безопасных размеров компактной укладки материала

I. ОБЩИЕ ПОЛОЖЕНИЯ

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

2. Определение возможности теплового самовозгорания грузов выполняется по результатам расчетов критической температуры и периода индукции процесса для температуры окружающей среды, превышающей верхнюю границу климатического перепада среднесуточных значений для регионов транзита (40 °С).

3. Определение класса опасности грузов выполняется в следующей последовательности:

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

б) для грузового пространства рассматриваемых транспортных средств вычисляется параметр δ0;

в) рассчитывается критическая температура Ткр окружающей среды;

г) при Ткр < 40 °С определяется период индукции процесса для температуры среды 40 °С (минимальное значение для всех условий транспортирования);

д) определяется класс опасности груза при перевозке рассматриваемым транспортным средством. Для других условий транспортирования повторяют расчеты по п.п. 3б - 3г;

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

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

II. ОПРЕДЕЛЕНИЕ КИНЕТИЧЕСКИХ ПАРАМЕТРОВ ПРОЦЕССА ТЕРМООКИСЛЕНИЯ МАТЕРИАЛОВ ПО ЭКСПЕРИМЕНТАЛЬНЫМ ДАННЫМ

2.1. Аппаратура

Аппаратура для определения кинетических параметров процесса термоокисления материалов включает в себя следующие элементы:

2.1.1. Термостат вместимостью рабочей камеры не менее 40 дм3 с терморегулятором, позволяющим поддерживать постоянную температуру от 60 до 500 °С с погрешностью не более ±1 °С.

2.1.2. Корзинки кубической или цилиндрической формы высотой 15, 30, 35, 50, 70, 100, 140 и 200 мм. Диаметр цилиндрической корзинки должен быть равен ее высоте. Материалом для корзинок служит сетка из латуни или нержавеющей стали для сыпучих материалов (с размером ячеек не более 1 мм) или листовая нержавеющая сталь толщиной не более 1 мм - для плавящихся веществ.

2.1.3. Термоэлектрические преобразователи (термопары ТХА и ТХК) с максимальным диаметром рабочего спая не более 0,8 мм.

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

2.1.5. Весы лабораторные с наибольшим пределом взвешивания 1000 г и точностью взвешивания 0,01 г.

2.2. Подготовка и проведение испытаний

2.2.1. К корзинкам крепят по три термоэлектрических преобразователя таким образом, чтобы один конец одной термопары находился внутри корзинки в ее центре, а второй - на расстоянии не более 5 мм от внешней ее стороны (на высоте центра корзинки). Эти термопары соединяют по дифференциальной схеме, с тем чтобы они измеряли разность температур между образцом материала и температурой рабочей камеры. Для фиксирования температуры в термостате (температуры испытаний) рабочий конец третьей термопары располагают на расстоянии (30±1) мм от стенки корзинки на высоте ее центра.

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

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

2.2.4. Корзинку помещают в центр термостата, нагретого до заданной температуры (например, 200 °С) и наблюдают за изменением температуры в центре образца.

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

2.2.6. Если при первом испытании самовозгорание не происходит в течение времени, указанного в табл. 1, то испытание с новым образцом материала того же размера проводят при температуре на 20 °С больше заданной. Если самовозгорание произошло, то испытание проводят при температуре меньшей на 10 °С.

Таблица 1

Размер образца, мм

Продолжительность испытаний, ч

35

6

50

12

70

24

100

48

140

96

200

192

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

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

Таблица 2

Размер образца, мм

Температура самовозгорания Ткр

°С

К

 

 

 

2.3. Расчет параметров кинетики термоокисления материалов

Исходными данными для определения параметров кинетики термоокисления являются:

- данные табл. 2 для критической температуры самовозгорания Tкр º Т0 (К) образцов размером D (м);

- коэффициент теплопроводности материала λ, Вт·м-1·К-1;

- теплоемкость исследуемого материала с, Дж·кг-1·К-1;

- теплота реакции Q, Дж·кг-1.

Расчет выполняется в следующем порядке.

2.3.1. Для каждого размера образца рассчитать число Рэлея по уравнению:

                                                                                                        (1)

где g - ускорение силы тяжести, м·с-2;

v - кинематическая вязкость воздуха при температуре Т0, м2·с-1;

а - температуропроводность воздуха при температуре Т0, м2·с-1;

D - высота образца, м;

R - универсальная газовая постоянная, Дж·моль-1·К-1);

Т0 - температура рабочего пространства термостата, К;

Е - энергия активации реакции окисления. Допускается принимать равной 100 кДж·моль-1.

Для облегчения расчетов зависимость комплекса g/va от температуры в диапазоне Т0 = (350÷800) К может быть рассчитана по формуле

                                                                                                       (2)

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

при 5·102 < Ra 2·107                                                                                                                                                       

                                                                                     (3)

при < Ra > 2·107                                                                                                                                                                      

                                                                                   (4)

где σ=5,67·10-8 - постоянная Стефана-Больцмана, Вт·м-2·К-1.

Коэффициент теплопроводности воздуха при температуре Т0 может быть определен по формуле

λв = 6,98·10-3 + 6,41·10-5 Т0.                                                                                       (5)

2.3.3. По величине α, коэффициенту теплопроводности материала λ и половине высоты r = D/2 вычисляют критерии Био для каждого образца:

                                                                                                                  (6)

2.3.4. Функцию φ(Bi), учитывающую интенсивность теплообмена образца с воздухом, определяют по уравнению

                                                   (7)

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

                                                                                                                    (8)

                                                                                                                   (9)

2.3.6. С учетом интенсивности теплообмена и характеристик реакции для каждого размера образца материала рассчитывают критическое значение параметра Франк-Каменецкого:

δкр = δ0φ(Вi) (1+β) (1+2,4γ2/3),                                                                                   (10)

где δ0 - критическая величина параметра δ при интенсивном теплообмене, равная 2,52 для образцов кубической формы и 2,76 - для цилиндра с высотой, равной диаметру.

Результаты вычислений по формулам (1) - (10) сводят в табл. 3.

Таблица 3

Размер r, м

Т0, К

Ra

α, Вт·м-2·К-1

Bi

φ(Bi)

β

γ

δкр

 

 

 

 

 

 

 

 

 

2.3.7. Зависимость критического значения параметра Франк-Каменецкого δкр от кинетических параметров реакции окисления

                                                                                              (11)

записывают в виде

                                                                                                               (12)

где

                                                                                                              (13)

                                                                                                                 (14)

ρ - плотность упаковки материала, кг·м-3;

k0 - константа скорости реакции, 1·с-1.

2.3.8. По уравнению (13) для каждого размера образца рассчитывают величину М. С учетом значений М и N по уравнению (12) методом наименьших квадратов или программ обработки экспериментальных данных для персональных компьютеров Eureka, Curve Expert 1.3, Mathematica 3.0, Mathematica 4.0 и других определяют численные значения N и энергию активации Е.

2.3.9. Вычисляют предэкспоненциальный множитель реакции окисления Qk0/λ путем деления N на Е. Данные расчетов по уравнениям (12) - (14) сводят в табл. 4.

Таблица 4

Размер r, м

Т0, К

М, Дж·м·К·кг-1·моль-1

N, Дж·м·К·кг-1·моль-1

Е, Дж·моль-1

Qk0/λ, м·К·кг-1

 

 

 

 

 

 

2.3.10. Если величина энергии активации, вычисленная в п. 2.3.8, отличается от ранее принятой и равной 100 кДж·моль-1 более чем на 5 %, расчеты по п.п. 2.3.1 - 2.3.10 необходимо повторить с новым значением энергии активации. Процесс итераций необходимо выполнять до тех пор, пока энергии активации в начале и конце расчета не будут отличаться менее чем на 5 %.

III. РАСЧЕТ УСЛОВИЙ САМОВОЗГОРАНИЯ ГРУЗОВ

3.1. Расчет параметра Франк-Каменецкого δ0 для грузового пространства транспортных средств

Исходными данными для расчета параметра δ0 являются форма и размеры грузового пространства.

Расчет выполняется в следующем порядке.

3.1.1. Вычисляют отношение квадрата характерного размера грузового пространства r2 (минимального размера по одной из осей координат) к квадрату эквивалентной сферы Франк-Каменецкого R02 по одному из соотношений.

Прямоугольный цилиндр (цилиндрические канистры, бочки и т.п.) радиусом r, высотой 2d, p = r/d:

                                                                                       (15)

Прямоугольный брус (мешок, контейнер, вагон, штабель упаковок материала, заполняемое насыпью грузовое пространство) со сторонами 2а, 2b, 2с, р = b/а, q = c/a:

                   (16)

где

3.1.2. Находят радиус эквивалентной сферы Семенова по формуле

                                                                                                                     (17)

где V - объем упаковки материала, м3;

S - ее внешняя поверхность, м2.

3.1.3. Определяют отношение квадратов радиуса эквивалентных сфер Франк-Каменецкого и Семенова:

                                                                                                  (18)

3.1.4. Вычисляют фактор формы для заданной геометрии упаковки материала:

j = 3σ-1.                                                                                                                      (19)

3.1.5. Находят функцию F(j) по формуле

                                                                                                            (20)

3.1.6. Рассчитывают величину параметра Франк-Каменецкого с помощью формулы

                                                                                         (21)

Результаты расчета значений δ0 для основных типоразмеров грузового пространства различных транспортных средств приводятся в прил. 4.

3.2. Расчет критической температуры

Исходными данными для расчета критической температуры при транспортировании веществ и материалов являются:

- плотность упаковки материала, кг·м-3;

- ксоффициент теплопроводности материала λ, Вт·м-1·К-1;

- теплоемкость исследуемого материала с, Дж·кг-1·К-1;

- теплота реакции Q, Дж·кг-1;

- энергия активации реакции окисления Е, Дж·моль-1;

- предэкспоненциальный множитель Qk0/λ, м·К·кг-1.

Расчет выполняется в следующем порядке.

3.2.1. Для заданной формы упаковки материала из прил. 4 выбрать или рассчитать в соответствии с разд. 3.1 величину критерия Франк-Каменецкого δ0.

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

3.2.3. По формуле (1) п. 2.3.1 вычислить значение критерия Рэлея для заданного размера упаковки материала.

3.2.4. Рассчитать коэффициент теплоотдачи по уравнению (4) п. 2.3.2 и величину критерия Био по формуле (6).

3.2.5. Определить численное значение функции φ(Bi) по уравнению (7).

3.2.6. По формулам (8) и (9) найти величину параметров β и γ.

3.2.7. Рассчитать критическое значение параметра Франк-Каменецкого по уравнению (10).

3.2.8. Подставить величину δкр в уравнение (11) и найти новое значение температуры Т0.

3.2.9. Используя это значение Т0, повторить расчет параметров по п.п. 3.2.2 - 3.2.8.

3.2.10. Указанную процедуру расчета продолжать до тех пор, пока предыдущее и последующее значения температуры будут отличаться друг от друга менее чем на 1 °С. За критическую температуру принимается результат последнего расчета.

3.3. Расчет времени индукции

Расчет времени индукции производится в том случае, когда критическая температура для самовозгорания грузового пространства, заполненного материалом, ниже 40 °С (согласно п. разд. I).

Исходными данными для расчета являются:

- температура окружающей среды Т0, К (согласно п. разд. I равна 313 К);

- критическая температура самовозгорания для заданного грузового пространства транспортного средства Ткр, К;

- фактор формы упаковки материала j;

- размер упаковки r, м;

- плотность упаковки материала, кг·м-3;

- коэффициент теплопроводности материала λ, Вт·м-1·К-1;

- теплоемкость исследуемого материала с, Дж·кг-1·К-1;

- теплота реакции Q, Дж·кг-1;

- энергия активации реакции окисления Е, Дж·моль-1;

- предэкспоненциальный множитель Qk0/λ, м·К·кг-1.

Расчет выполняется в следующем порядке.

3.3.1. По температуре Т0 вычислить параметры β и γ с помощью формул (8) и (9).

3.3.2. По уравнению (4) определить коэффициент теплоотдачи α и найти величину параметра Био по формуле (6).

3.3.3. По формуле (11) рассчитать параметр δ, соответствующий температуре окружающей среды Т0, и параметр δкр для критической температуры Ткр.

3.3.4. Вычислить относительное удаление от предела воспламенения

                                                                                                                       (22)

и функции

                                                                             (23)

                                                                      (24)

3.3.5. Рассчитать безразмерное время индукции с помощью уравнений (22) - (24) по выражению

                                                                                   (25)

3.3.6. Определить размерное время индукции tи (с) по формуле

                                                                                                          (26)

IV. ОПРЕДЕЛЕНИЕ КЛАССА ОПАСНОСТИ САМОВОЗГОРАЮЩИХСЯ ГРУЗОВ И РАЗРАБОТКА ПРОФИЛАКТИЧЕСКИХ МЕРОПРИЯТИЙ

4.1. Определение класса опасности грузов

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

Исходными данными для классификации опасности грузов являются:

- факт самовозгорания испытываемых согласно разд. 2.2 образцов материала;

- критическая температура самовозгорания материала в грузовом пространстве транспортного средства Ткр, °С;

- значения периода индукции процесса для температуры среды Т0 = 40 °С.

Определение класса опасности грузов осуществляется в следующем порядке.

4.1.1. Не самовозгорающиеся в условиях испытаний согласно разд. 2.2 и негорючие в соответствии с ГОСТ 12.1.044 [2] материалы следует относить к непожароопасным подклассам и категориям (в соответствии с ГОСТ 19433).

4.1.2. Если значения Ткр ≥ 40 °С, перевозимый груз следует отнести к категории опасности 921.

4.1.3. Если значения критической температуры Ткр < 40 °С, сравнивают период индукции теплового самовозгорания груза τ40 (согласно разд. 3.3) для Т0 = 40 °С со временем осуществления грузоперевозок τпep:

а) при τ40 > τпep груз следует отнести к низкой степени опасности подкласса 4.2;

б) при τ40 ≤ τпep груз следует отнести к средней степени опасности подкласса 4.2.

В случае (б) необходимо обеспечение пожарной безопасности грузоперевозок за счет использования безопасных размеров компактной укладки материала, определенных согласно разд. 4.2 для Т0 = 40 °С. Если реализация этих мер невозможна (малая величина безопасного размера и т.п.), отсутствие самовозгорания материала можно обеспечить при транспортировании груза в течение времени не более 0,8 τ40.

4.2. Определение безопасных размеров компактной укладки материала

Расчеты выполняются для случая, оговоренного в п. 4.1.3б. Исходными данными для расчета критического размера укладки материала являются:

- энергия активации реакции окисления Е, Дж·моль-1;

- предэкспоненциальный множитель Qk0/λ, м·К·кг-1;

- температура среды Т0, К;

- δ0 грузового пространства;

- плотность упаковки материала, кг·м-3;

- теплоемкость исследуемого материала с, Дж·кг-1·К-1;

- теплота реакции Q, Дж·кг-1.

Расчет выполняется в следующем порядке.

4.2.1. Для заданной температуры Т0 (см. пп. 4.1.3б и 4.1.5б) по формулам (8) и (9) вычислить параметры β и γ.

4.2.2. Рассчитать параметр δкр по формуле

δкр = δ0 (l + β)(l + 2,4 γ2/3).                                                                                         (27)

4.2.3. В первом приближении критический размер найти из выражения

                                                                                                     (28)

4.2.4. По уравнению (4) вычислить коэффициент теплоотдачи α.

4.2.5. Найти величину параметра Био по формуле (6) и рассчитать значение функции φ(Bi) по уравнению (7).

4.2.6. Определить величину параметра δкр по уравнению (10).

4.2.7. По формуле (28) вычислить новое значение размера упаковки материала.

4.2.8. Используя это значение r, повторить расчет параметров по п.п. 4.2.4 - 4.2.7.

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

4.2.10. В качестве безопасного значения критического размера в соответствии с требованиями ГОСТ 12.1.004 [3] принимаем 0,8r. Так как определенный критический размер является половиной минимального размера скопления, за безопасный размер укладки материала (требуемый минимальный размер компактной укладки груза) принимается величина 1,6r.

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

1. ГОСТ 19433-88. Грузы опасные. Классификация и маркировка.

2. ГОСТ 12.1.044-89. Пожаровзрывоопасность веществ и материалов. Номенклатура показателей и методы их определения.

3. ГОСТ 12.1.004-91*. ССБТ. Пожарная безопасность. Общие требования.

ПРИЛОЖЕНИЕ 1

Пример расчета кинетических параметров

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

Коэффициент теплопроводности ржаной муки λ = 0,152 Вт·м-1·К-1; теплоемкость с = 1050 Дж·кг-1·К-1; тепловой эффект реакции Q = 1,67·107 Дж·кг-1; плотность упаковки материала ρ = 655 кг·м-3.

Расчет проведем для образца размером D = 35 мм. Данные для других размеров получим, повторяя приведенную ниже последовательность расчета.

1. По формуле (1) вычислим числа Рэлея. Принимая в первом приближении энергию активации Е = 100 000 Дж·моль-1, для образца размером D = 35 мм получим:

2. Коэффициент теплоотдачи α найдем по уравнению (3),

где теплопроводность воздуха определена по формуле (5)

λв = 6,98·10-3 + 6,41·10-5Т0 = 6,98·10-3 + 6,41·10-5·492 = 0,038 Вт·м-1·К-1.

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

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

5. Рассчитаем параметры β и γ:

6. Критическое значение параметра Франк-Каменецкого будет равно:

δкр = δ0 φ(Bi) (1 + β)(1 + 2,4 γ2/3) = 2,52·0,625(1 + 0,04) (l + 2,4·0,001262/3) = 1,68,

где δ0 - критическая величина параметра δ, соответствующая интенсивному теплообмену, для образцов кубической формы равная 2,52. Результаты вычислений для всех образцов представлены в табл. П 1.1.

Таблица П 1.1

Размер r, м

Т0, К

Ra

α, Вт·м-2·К-1

Bi

φ(Bi)

β·10-2

γ·10-3

δкр

0,0175

492

7683

32,5

3,7

0,625

4,00

1,26

1,68

0,025

474

24740

29,2

4,8

0,687

3,94

1,17

1,84

0,035

459

74272

26,5

6,1

0,736

3,82

1,05

1,97

0,05

444

238607

24,4

8,0

0,790

3,69

0,98

2,10

0,07

431

716774

21,9

10,1

0,828

3,58

0,93

2,21

0,10

418

2302797

19,9

13,1

0,863

3,47

0,87

2,47

7. По формуле (13) рассчитаем величину М:

8. С помощью этих значений и уравнения (12) методом наименьших квадратов определим численные значения N и энергию активации Е.

9. Вычислим предэкспоненциальный множитель реакции окисления Qk0/λ путем деления N на Е. Данные расчетов по п.п. 7-9 сведем в табл. П 1.2.

Таблица П 1.2

Размер r, м

Т0, К

м,

Дж·м·К

кг·моль

N,

Дж·м·К

кг·моль

Е, Дж·моль-1

Qk0/λ, м·К·кг-1

0,0175

492

1,68·107

5,76·1016

88 054

6,55·10-11

0,025

474

8,39·106

0,035

459

4,06·106

0,05

444

1,98·106

0,07

431

1,01·106

0,10

418

5,17·105

10. Повторяя расчет по п.п. 1-9 с величиной энергии активации Е = 88 054, найдем новые значения энергии активации Е = 88 068 Дж·моль-1 и предэкспоненциального множителя Qk0/λ = 6,53 1011 м·К·кг-1. Так как последние величины практически не отличаются от предыдущих, процесс итераций следует прекратить и за кинетические параметры реакции окисления ржаной муки принять Е = 88 054 Дж·моль-1 и Qk0/λ = 6,53 1011 м·К·кг-1.

ПРИЛОЖЕНИЕ 2

Параметры кинетики процесса термоокисления некоторых материалов

Материал

Е, Дж·моль-1

Qk0/λ, м·К·кг-1

Дрожжи кормовые

256 212

5,11·1030

Дрожжи товарные (влажн. 8 %)

106 144

5,94·1014

Дрожжи товарные (влажн. 15 %)

92 285

1,67·1013

Мука ржаная

88 054

6,55·1011

Мука костная крупнозернистая

50 740

2,46·108

Пыль костной муки

71 623

8,67·1010

Сено

179 050

8,67·1022

Хлопок

121 825

1,37·1016

Льно-джут

63 134

4,826·109

Метионин кормовой

88 278

2,35·1010

Уголь-сырец марки А

71 280

4,31·1011

Уголь ОУ-А

101 458

2,82·1015

Уголь ОУ-Б

97 650

2,01·1013

Уголь-сырец после сортировки

101 450

2,82·1015

Технический углерод К 354

56 943

4,343·1012

Технический углерод Н 990

90 732

5,665·1014

ПРИЛОЖЕНИЕ 3

Пример расчета δ0 грузового пространства

Рассчитать δ0 грузового пространства в железнодорожном вагоне. Вагон представляет собой параллелепипед шириной 2,75 м, длиной 15,7 м и высотой 4,692 м.

1. Отношение квадратов полувысоты вагона к эквивалентной сфере Франк-Каменецкого рассчитаем как для прямоугольного бруса по выражению (16)

где

а, b, c - половины сторон бруса, а - наименьшая сторона; р = b/а, q = c/a.

Подставляя р, q в эти равенства, получим

2. Средний радиус эквивалентной сферы Семенова равен

где V, S - объем и поверхность грузового пространства вагона.

3. На основании формул (16) - (18) отношение квадратов радиусов эквивалентных сфер Франк-Каменецкого и Семенова будет равно

4. Фактор формы прямоугольного бруса согласно уравнению (19) составит:

j = 3σ - 1 = 3·0,945 - 1 = 1,83.

5. В соответствии с формулой (20) функция F(j) будет равна

6. Из уравнения (21) значение параметра δ0 равно

ПРИЛОЖЕНИЕ 4

Значения δ0 для основных типоразмеров грузовых пространств

Тип грузового пространства (размеры)

δ0

Мягкий контейнер (0,95×0,95×1,3 м)

2,19

Крытый вагон (2,70×2,75×13,80 м)

1,779

Крытый вагон (2,70×2,76×13,80 м)

1,773

Крытый вагон (2,70×2,765×15,724 м)

1,759

Крытый вагон (2,75×4,692×15,70 м)

1,23

Вагон-хоппер (2,64×3,112×13,37 м)

2,178

Вагон-хоппер (3,112×4,722×17,50 м)

1,319

Универсальный контейнер (2,27×2,438×12,12 м)

1,69

Универсальный контейнер (1,95×2,10×2,543 м)

2,11

Универсальный контейнер (2,352×2,70×12,035 м)

1,60

Автофургон - 25 мягких контейнеров в 1 ярус

1,051

Трюм сухогруза (6,0×17,2×21,4 м)

1,112

Трюм сухогруза (6,0×17,2×26,6 м)

1,088

Трюм сухогруза (6,0×17,2×28,5 м)

1,081

ПРИЛОЖЕНИЕ 5

Пример расчета критической температуры

Рассчитать критическую температуру окружающей среды при транспортировании льно-джута в железнодорожных вагонах. Вагон представляет собой параллелепипед шириной 2,75 м, длиной 15,7 м и высотой 4,692 м.

Исходными данными для расчета являются:

- плотность материала ρ = 134 кг·м-3;

- коэффициент теплопроводности материала λ = 0,049 Вт·м-1·К-1;

- теплоемкость исследуемого материала с = 1505 Дж·кг-1·К-1;

- теплота реакции Q = 17 501 000 Дж·кг-1;

- энергия активации реакции окисления Е = 63 134 Дж·моль-1;

- предэкспоненциальный множитель Qk0/λ = 4,826×109 м·К·кг-1.

1. Подставим полученную в прил. 3 величину δ0 для грузового пространства вагона в формулу

вместо δкр и, решив его относительно Т0, получим нулевое приближение для этой температуры равное 302 К.

2. С помощью полученной величины рассчитаем:

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

 - параметр, характеризующий реакцию окисления.

3. Так как для размеров упаковок, превышающих 1 м φ(Bi) » 1, безразмерное значение критического параметра Франк-Каменецкого, учитывающего выгорание вещества и свойства реакции горения, определим по формуле

δкр = δ0 (1+2,4γ2/3) (1+β) = 1,23 (1 + 2,4·0,0012/3) (1 + 0,039) = 1,332.

4. Решая уравнение (11) относительно температуры получим Т0 =304 К или 31 °С. Следовательно, при перевозке льно-джута в железнодорожных вагонах при температурах воздуха больших 31 °С возможно возникновение самовозгорания транспортируемого продукта.

ПРИЛОЖЕНИЕ 6

Пример расчета времени индукции

Рассчитать время индукции при перевозке льно-джута в вагоне при температуре 40 °С (313 К). Фактор формы вагона j = 1,83.

Исходными данными для расчета являются:

- температура перевозки материала Т0 = 313 К;

- критическая температура самовозгорания для заданного размера и формы упаковки материала Ткр = 304 К;

- фактор формы упаковки материала j = 1,83;

- размер упаковки r = 1,375 м;

- плотность упаковки материала ρ = 134 кг·м-3;

- коэффициент теплопроводности материала λ = 0,042 Вт·м-1·К-1;

- теплоемкость исследуемого материала с = 1505 Дж·кг-1·К-1;

- теплота реакции Q = 17 501 000 Дж·кг-1;

- энергия активации реакции окисления Е = 63 134 Дж·моль-1;

- предэкспоненциальный множитель Qk0/λ = 4,826×109 м·К·кг-1.

1. По температуре Т0 вычислим параметры β и γ с помощью формул (8) и (9):

2. По уравнениям (1) и (2) найдем число Рэлея:

3. Коэффициент теплоотдачи α определим по уравнению (4):

где теплопроводность воздуха рассчитана по формуле (5):

λв = 6,98·10-3 + 6,41·10-5Т0 =

= 6,98·10-3 + 6,41·10-5·313 = 0,0269 Вт·м-1·К-1.

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

5. По формуле (11) рассчитаем параметр δ, соответствующий температуре Т0, и параметр δкр для критической температуры Ткр:

6. Вычислим относительное удаление от предела воспламенения:

и функции

7. Рассчитаем безразмерное время индукции:

8. Определим размерное время индукции tи (с), по формуле

ПРИЛОЖЕНИЕ 7

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

Определить класс опасности льно-джута при перевозке железной дорогой в крытых вагонах в течение 15 суток.

1. Рассматриваемый материал самовозгорается в условиях лабораторных испытаний, соответствующих требованиям разд. 2.2. Параметры кинетики процесса термоокисления определены согласно разд. 2.3 (см. прил. 2).

2. Наименьшее значение критической температуры процесса самовозгорания льно-джута в крытых вагонах определяем в соответствии с разд. 3.1 и 3.2 для модификации вагонов с наибольшими размерами. Согласно прил. 3 расчет ведем для вагона с размерами 2,75×4,692×15,70 м. Значение Ткр составляет 304 К или 31 °С.

3. В соответствии с п. 4.1.3 определяем время индукции до самовозгорания льно-джута в вагоне при температуре среды Т0 = 40 °С (согласно разд. 3.3). Это значение равно 53,07 суток.

4. Так как время транспортирования материала не превышает 15 суток (< 51,7 суток), то согласно п. 4.1.3а льно-джут при перевозках в крытых железнодорожных вагонах относится к низкой степени опасности подкласса 4.2 (согласно классификации ГОСТ 19433).

ПРИЛОЖЕНИЕ 8

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

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

Расчет проведем для верхней границы диапазона климатического перепада температур воздуха в средней полосе России равной 40 °С или 313 К.

Исходными данными для расчета критического размера являются:

- плотность упаковки материала ρ = 660 кг·м-3;

- коэффициент теплопроводности материала λ = 0,14 Вт·м-1·К-1;

- теплоемкость исследуемого материала с = 780 Дж·кг-1·К-1;

- теплота реакции Q = 350 000 Дж·кг-1;

- энергия активации реакции окисления Е = 50 740 Дж·моль-1;

- предэкспоненциальный множитель Qk0/λ =2,46·108 м·К·кг-1.

1. Выбираем вагон с наибольшим размером грузового пространства, для которого согласно прил. 4 известно, что δ0 = 1,23.

2. Для температуры Т0 = 313 К вычислим параметры β и γ по уравнениям (8) и (9):

3. Считая φ(Bi) » 1, определим параметр δкр по формуле (10):

δкр = δ0 (1+β) (1+2,4γ2/3) = 1,23 (1 + 0,0513)(1 + 2,4·0,03582/3) = 1,63.

4. В первом приближении минимальный размер найдем из выражения (28):

5. По уравнениям (1) и (2) для полученного размера вычислим число Рэлея:

6. Коэффициент теплоотдачи α найдем по уравнению (4):

где теплопроводность воздуха определена по формуле (5):

λв = 6,98·10-3 + 6,41·10-5Т0 =

= 6,98·10-3 + 6,41·10-5·313 = 0,027 Вт·м-1·К-1.

7. Вычислим критерий Био, соответствующий предварительному размеру компактной укладки материала в вагоне:

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

9. Критическое значение параметра Франк-Каменецкого будет равно:

δкр = δ0 φ(Bi) (1 + β)(1 + 2,4 γ2/3) = 1,23·0,895(1 + 0,0513) (l + 2,4·0,03582/3) = 1,46.

10. По формуле (28) вычислим новое значение размера насыпи материала:

11. Так как последнее значение отличается от предшествующего менее чем на 5 %, считаем его критическим размером компактной укладки согласно п. 4.2.9.

12. В соответствии с п. 4.2.10 безопасный характеристический размер компактной укладки составит 0,8r = 0,165 м. Полный безопасный размер составит 0,33 м.

13. Безопасная для теплового самовозгорания укладка костной муки при перевозке в крытых вагонах будет соответствовать штабелированию на поддонах в несколько ярусов с высотой укладки материала в каждом ярусе не более 0,33 м. Безопасные условия транспортирования костной муки железной дорогой в холодное время года могут определяться дополнительно согласно п. 4 разд. I методики.