УДК 614.84:681.5
МОДЕЛИРОВАНИЕ ПРОЦЕССА НАГРЕВА ЭЛЕКТРОПРОВОДНИКА
Дмитриченко A.C., к.т.н., доцент, Иванович A.A., Чайчиц Н.И., Зинкевич Г.Н.
В статье предложена математическая модель по определению установившейся температуры нагрева токоведущей жилы изолированного проводника в зависимости от конструкции проводника, теплофизических параметров изоляции и окружающей среды и величины протекающего тока. Приведены экспериментальные данные по исследованию стационарного температурного поля проводников при различных величинах тока. В результате сравнения расчетных и экспериментальных данных обоснована возможность использования разработанной математической модели для уточнения длительно допустимых токовых нагрузок проводников [1, 2], при заданных допустимых температурах нагрева токоведущей жилы, с учетом применения новых материалов изоляции, оболочки и наружных защитных покровов.
Кабельное изделие в процессе эксплуатации подвергается термическому воздействию со стороны электрического тока, протекающего по токоведущей жиле. В настоящее время регламентирована допустимая температура нагрева материала токоведущей жилы проводников при различных видах изоляции, оболочки и наружных защитных покровов.
В соответствии с существующей теорией тепломассообмена кабельное изделие при открытых способах прокладки можно рассматривать как часть термодинамической системы. Совокупность этой термодинамической системы в процессе эксплуатации составляют электрический проводник как источник и окружающая среда как приемник теплоты.
Источник теплоты можно рассматривать как бесконечный элемент цилиндрической формы, имеющий характеристики металла.
Приемник теплоты можно представить в виде двух различных элементов: твердой среды в виде материалов изоляции, оболочки, защитного покрова кабельного изделия и газовоздушной среды.
При анализе термодинамической системы «электрический проводник-среда» можно выделить следующие элементы теплового взаимодействия:
- изменение внутренней энергии источника тепла, обусловленное электротермическими преобразованиями в проводнике при протекании по нему электрического тока. Энергетический (или тепловой) потенциал, который может приобрести проводник, без перехода в другие агрегатные состояния, определяется теплофизическими свойствами материала и геометрическими параметрами токоведущей жилы;
- передача теплового потенциала от источника в первый элемент рассматриваемой среды посредством теплопроводности. В результате часть переданной тепловой энергии источника накапливается в объеме твердого элемента среды. Тепловой потенциал определяется величиной градиента температуры и коэффициентом теплопроводности твердой среды;
- передача тепловой энергии от твердого элемента к газовоздушной среде посредством процесса конвективного теплообмена. Тепловой потенциал, который прошел через первый элемент, поступает в окружающую среду, вызывая увеличение ее энтропии. Температура на границе конвективного теплообмена будет определяться величиной градиента температуры, теплофизическими параметрами газовоздушной среды (влажность, давление, коэффициент объемного расширения, коэффициент кинематической вязкости, коэффициент температуропроводности) и режимами движения газовоздушной среды (направление и относительная скорость).
При достижении высоких температур нагрева необходимо дополнительно учитывать лучистый теплообмен.
Состояние рассматриваемой термодинамической системы «кабель-среда» может быть равновесным (стационарным) и неравновесным (нестационарным).
Равновесное состояние термодинамической системы будет реализовано в случае, если температуры на границах рассматриваемых участков достигнут установившихся значений. Неравновесное состояние термодинамической системы наблюдается при нестационарном температурном поле и происходит, как правило, в начальный момент времени или при изменении температуры источника теплоты.
Рассматривая термодинамическую систему «электрический проводник-среда» в рамках процесса эксплуатации кабельных изделий, прокладываемых открыто, можно сказать, что нормальным условиям эксплуатации соответствует равновесное состояние термодинамической системы, при установленных нормативными документами допустимой температуре токоведущей жилы и величине тока нагрузки. Поэтому научный и практический интерес представляет аналитическое решение задачи по определению установившейся температуры токоведущей жилы проводника в зависимости от заданных условий эксплуатации и величины тока нагрузки. Аналогичная задача рассматривалась в работе [3], но для неизолированного проводника.
Для вывода аналитической зависимости Т=€(1,т) рассмотрим задачу о нагреве изолированного проводника при следующих допущениях:
- в силу высокой теплопроводности металлической жилы она нагревается мгновенно и имеет одинаковую температуру по сечению;
- тепло, выделяющееся в металлической жиле, нагревает ее мгновенно до температуры Т;
- распределение температуры внутри изоляции, оболочки и защитного покрова кабельного изделия определяется из уравнения теплопередачи через многослойную цилиндрическую стенку;
- коэффициенты теплопроводности и теплоемкости элементов кабельного изделия практически не изменяются в диапазоне температур эксплуатации.
Исходя из закона сохранения энергии и второго закона термодинамики, запишем уравнение теплового баланса в дифференциальной форме для рассматриваемой системы:
(1)
где (1С)Ж - количество теплоты, идущее на нагрев токоведущей жилы;
с1С)и - количество теплоты, обусловленное тепловым проявлением тока;
¿РсР - количество теплоты, рассеиваемое в результате процесса теплопередачи.
Уравнение (1) можно представить в виде:
Сжхр хУ хс1Т=:(12чхК~КхЬх(Т-Тср))х(1т5 (2)
Ж Ж Ж г
где Сж - удельная массовая теплоемкость токоведущей жилы, Дж/(кг*К); рж- плотность токоведущей жилы, кг/м3;
"3
Уж - объем токоведущей жилы, м ; I - сила установившегося тока, А; Я - сопротивление токоведущей жилы, Ом; Ь - длина токоведущей жилы, м; Т - текущая температура поверхности жилы, К; Тср - температура окружающей среды, К; т - время нагрева провода, с;
К - линейный коэффициент теплопередачи, Вт/(м2-К)
К = + - ), (3)
где а - коэффициент теплоотдачи от наружной поверхности проводника, Вт/(м2*К);
?ч - коэффициент теплопроводности ьго слоя твердой среды, Вт/(м*К); т{ - радиус ьго слоя твердой среды, м.
Уравнение теплового баланса (2) на единицу длины кабеля с учетом зависимости сопротивления жилы от температуры примет вид:
С хр х8 хс1Т-(12хрпх(1 + рх(Т-Т )/8 )-Кх(Т-Т ))хёт (4)
Ж ГЖ Ж 4 ги ч г \ Ср/ ж/ \ Ср//
где 8Ж - площадь поперечного сечения токоведущей жилы, м2; ро - удельное сопротивление жилы при температуре Тср, Ом-м; Р - коэффициент температурного изменения сопротивления жилы, К"1.
Уравнение (4) представляет собой дифференциальное уравнение с разделяющимися переменными.
Разделим переменные и проинтегрируем правую и левую части уравнения.
С хр х8 хсГГ ... |-Ж-ж-ж-----_ = (5)
(р Хрх12/3 -К)хТ + рп х!2 /8 х(1-рхТ ) + К х Т
После интегрирования с учетом начальных условий т = О, Т = Тср получаем выражение для определения температуры нагрева токоведущей жилы:
РЬх(Зх12/^-К
Таким образом, выражение (6) позволяет определить температуру нагрева поверхности токоведущей жилы проводника в зависимости от величины протекающего тока, времени нагрева, теплофизических и геометрических параметров кабеля и свойств окружающей среды.
Для проверки адекватности полученной модели авторами были проведены экспериментальные исследования температурных характеристик проводников.
Проведение испытаний и отбор образцов проводились в соответствие с методикой экспериментального исследования пожарной опасности элементов электропроводок [4].
В качестве образцов использовались кабели марки ВВГ-0.66 сечени-ем 4 мм и длиной 1,5 м. На исследуемом образце в средней его части (рисунок 1) с интервалом в 0,15 м устанавливались методом зачеканки три
термоэлектрических преобразователя типа ТХА1199 на токоведущую жилу и один на поверхность оболочки. Значение параметров регистрировалось многоканальным измерителем температуры типа Сосна-003. Информация с измерителя в автоматическом режиме передавалась и в дальнейшем обрабатывалась на ПЭВМ.
Для устранения влияния протекающего по токоведутцей жиле тока на показания термоэлектрических преобразователей последние имели гальваническую развязку по линии измерения.
Допустимый ток токоведущей жилы для исследуемых образцов определялся по [5] и составил 44 А. Проводники испытывались при кратностях тока 1х1доп и 1, 2х1доп. Испытания проводились на установке для определения то-ковременных характеристик кабелей [6].
Тоб Тж1 Тж2 ТжЗ
Токоведущая жила Оболочка
(Тоб, Тж1, Тж2? ТжЗ - термоэлектрические преобразователи)
Рисунок 1 - Схема расположения термоэлектрических преобразователей на исследуемом образце электропроводки
При проведении расчетов по выражению (6) коэффициент К определялся, исходя из фактического значения коэффициента теплоотдачи а, полученного в результате обработки экспериментальных данных с учетом критериальных уравнений теории подобия [7, 8]:
N11 = а х (7)
К
Сг = ёх(Тоб-Тср)хБУ(Тсрхд2) (8)
Рг = 3 / а (9)
где Боб - внешний диаметр оболочки кабеля, м;
А,ср- коэффициент теплопроводности воздушной среды, Вт/(м-К);
о
g = 9,81 - ускорение силы тяжести, м /с;
О - коэффициент кинематической вязкости воздушной среды, м2/с;
а - коэффициент температуропроводности воздушной среды, м /с.
Здесь N11, Ог и Рг соответственно критерии Нуссельта, Грасгофа и Прандтля.
Значения параметров окружающей среды, входящих в критерии Нуссельта и Прандтля, определялись при расчетной температуре Тт, равной полусумме температур поверхности оболочки Тоб и окружающей среды Тср при установившемся режиме, т. е.:
Тт=0ДТср+Тоб) (10)
При свободном движении воздуха возле горизонтального кабеля уравнение подобия имеет вид:
Ки = сх(ОгхРг)п, (11)
где сип- коэффициенты, определяемые в зависимости от произведения Gr.Fr [9].
Экспериментальные данные по стационарному режиму для исследованных образцов кабеля ВВГ 2x4 и теплофизические характеристики окружающей среды в условиях экспериментов приведены в таблице 1, а результаты расчета чисел Грасгофа, Прандтля, Нуссельта и фактического коэффициента теплоотдачи - в таблице 2.
Исходные данные для расчета температуры нагрева кабеля ВВГ 2x4, по выражению (6) приведены в таблице 3.
Таблица 1 - Экспериментальные данные
I, А Тж, К Тоб, к Тср? к Ооб, м а, м2/с X, Вт/(м-К) V, м2/с
44 341 330 292 0,0067 2,402Е~05 0,0274 1,675Е-05
52 359 342 292 0,0067 2,486Е-05 0,0279 1/734Е-05
Таблица 2 - Результаты расчета коэффициента а
I, А тт, К вг Рг Ог-Рг с п Ми а, Вт/(м2-К)
44 311 1368,56 0,700 955,14 0,54 0,25 3,00 12,28
52 317 1680,28 0,699 1172,01 0,54 0,25 3,16 13,16
Таблица 3 - Исходные данные для расчета температуры жилы
Сж? Рж, Ро, Р, Кб, Г Ж?
Дж/(кгК) кг/м3 Омм К"1 Вт/(м-К) Вт/(мК) м М м
390 8930 1.75Е-08 0,0039 0,17 0,17 1Д5Е-03 1,87Е-03 3,35Е-03
На рисунках 2, 3 графически показаны результаты расчетов, выполненных по формуле (6), и эксперимента,
Как видно из графиков, имеется сходимость расчетных и экспериментальных данных в области установившихся температур, что свидетельствует об адекватности полученной модели.
Однако на первоначальном этапе в режиме нестационарного теплообмена имеет место расхождение между результатами аналитического расчета и экспериментальными данными. Отклонение экспериментальных значений от расчетных можно объяснить инерционностью теплообменных процессов внутри элементов конструкции кабеля.
Таким образом, предложенное выражение представляет практический интерес и рекомендуется для инженерных расчетов по определению температуры нагрева жилы проводника в стационарном режиме при различных нагрузках и температурах окружающей среды.
Кроме того, полученное выражение можно использовать для уточнения длительно допустимых токовых нагрузок проводников, приведенных в [1,2, 5], при заданных допустимых температурах нагрева токоведущей жилы, с учетом применения новых материалов изоляции, оболочки и наружных защитных покровов.
70 4
10 -i
0 I ■ , ........................................г.......................... ,
О 60 120 180 240 300 360 420 480 540 600 660 720 780 840 900 960 1020 1080
Время испытания, с
-Расчетная кривая Экспериментальная кривая
Рисунок 2 - Сравнительные данные кабеля ВВГ 2x4 при токе 1доп = 44А
Время испытания, с
-Расчетная кривая —•—Экспериментальная кривая
Рисунок 3 - Сравнительные данные кабеля ВВГ 2x4 при токе I = 52А
ЛИТЕРАТУРА
1. Правила устройства электроустановок. - 6-е изд., перераб. и доп. - М.: Энергоатом-издат, 1986.-648 с.
2. Справочник по проектированию электрических сетей и электрооборудования / под ред. В.И. Круповича, Ю.Г. Барыбина, М.Л. Самовера. - 3-е изд. - М.: Энергоиздат, 1981.-408 с.
3. Смелков, Г.И. Пожарная опасность проводов при неполных коротких замыканиях и перегрузках / Г.И. Смелков // Пожарная профилактика в электроустановках: сб. науч. тр. - М.: ВНИИПО, 1981. - С. 3-18.
4. Методика экспериментального исследования пожарной опасности элементов электропроводок в зависимости от вида конструкции электропроводок, величины протекания электрического тока, режимов работы электрической сети. - НИИ ПБ и ЧС МЧС РБ, 2006.
5. ГОСТ 16442 - 80. Кабели силовые с пластмассовой изоляцией. Технические условия.
6. НПБ 9 - 2000. Кабели и провода электрические. Показатели пожарной опасности. Методы испытания.
7. Уонг, X. Основные формулы и данные по теплообмену для инженеров: справочник: пер. с англ. / X. Уонг. - М., Атомиздат, 1979. - 216 с.
8. Кошмаров, Ю.А. Термодинамика и теплопередача в пожарном деле / К).А. Кошмаров, М.П. Башкирцев. - М.: ВПТШ МВД СССР, 1987. - 444 с.
9. Башкирцев, М.П. Основы пожарной теплофизики / М.П. Башкирцев [и др.]. - 3-е изд. - М.: Стройиздат, 1984. - 200 с.
Поступила в редакцию 11 января 2007 г.