MODELING OF ELECTROMAGNETIC PROCESSES IN A MULTILAYER THREE-PHASE INDUCTION CYLINDRICAL SYSTEM

Full Text

Актуальной проблемой сложных технических систем, к которым относятся технологические объекты трубопроводного транспорта, является обеспечение эффективности и надежности их эксплуатации [1-4]. В работах [5, 6, 7] рассматриваются вопросы расчета, проектирования и практического применения систем электрообогрева трубопроводов, резервуаров и технологического оборудования в нефтегазовой промышленности. Приведены описания и характеристики установок электрообогрева, работающих в различных отраслях промышленности. Однако вопросы моделирования и расчета охватывают ограниченный класс нагревателей с источниками энергии в виде греющих кабелей или индукционных систем с однослойной нагрузкой в виде металлической стенки резервуара или трубопровода. В работе [8] выполнены исследования электромагнитных и тепловых процессов в установке технологического нагрева нефти при транспортировке по трубопроводам. Предложена конструкция многосекционного индукционного нагревателя для транспортировки высоковязких жидкостей, в частности высокопарафинистых нефтей. Показано, что вследствие низкой теплопроводности и высокой вязкости нефти при наличии технологических ограничений на максимальную температуру контактирующих поверхностей трубы и жидкости для трубопроводных систем высокой производительности индукционные нагреватели должны иметь несколько автономных секций с развитой поверхностью теплообмена. В работах [9, 10] рассмотрена конструкция индукционного нагревателя для транспортировки высоковязких нефтей по магистральным трубопроводным системам. Предлагаемая конструкция нагревателя с осесимметричными трубами, в которой внутренняя труба служит вытеснителем и позволяет вдвое увеличить площадь теплообмена при минимальном поперечном сечении потока жидкости. Однако и в этом случае, как показывает практика, общая длина нагревателя составляет не менее 8÷12 м. Повысить интенсивность теплообмена и тем самым уменьшить габариты нагревательной системы при той же производительности можно либо путем создания турбулентности потока за счет высоких скоростей, либо путем перемешивания жидкости в нагреваемом потоке. При практически реализуемых скоростях потока нефти в трубопроводных системах имеет место ламинарное течение жидкости. Решение указанной проблемы путем использования внешнего устройства для перемешивания осложняется наличием высокого давления в трубопроводе, которое предполагает полную герметизацию трубопроводной системы. В настоящей работе проблема уменьшения общей длины нагревателя при той же производительности трубопровода решается за счет повышения интенсивности теплообмена между стенкой тепловыделяющего цилиндра и потоком нагреваемой жидкости с помощью индукционного устройства, не нарушающего герметичность трубопроводной системы. Эскиз устройства представлен на рис. 1. Предлагаемая конструкция представляет собой систему, состоящую из трехфазного цилиндрического индуктора с замкнутым магнитопроводом. Внутри расточки индуктора находится труба, встроенная в трубопроводную систему. Внутри трубы располагается полый ферромагнитный цилиндр - ротор с беличьей клеткой. Ротор установлен в подшипниках качения. В полости ротора находится крыльчатка, закрепленная под определенным углом к плоскости сечения цилиндра. Назначением предлагаемой конструкции является преобразование электрической энергии, подводимой к индуктору, в механическую для перемешивания жидкости и в тепловую для нагрева. Перемешивание жидкости осуществляется за счет вращения ротора под воздействием трехфазного электромагнитного поля статора. При этом часть энергии, передаваемой индуктором, выделяется в виде тепла во внешней трубе и в роторе, обеспечивая дополнительный подогрев жидкости. Такой способ, совмещающий нагрев и перемешивание в одном устройстве, позволяет значительно ускорить процесс нагрева жидкости с низкой теплопроводностью и уменьшить общую длину нагревателя. Соотношение между мощностью нагрева и мощностью, идущей на вращение ротора и перемешивание жидкости, зависит от частоты тока источника питания, реологических свойств транспортируемой жидкости, электрофизических свойств материала и конструктивных параметров системы. Рис. 1. Конструкция устройства: 1 - магнитопровод статора; 2 - пазы обмотки статора; 3 - внешняя труба; 4 - ротор; 5 - крыльчатка; 6 - поток нагреваемой жидкости; 7 - слой изоляции Для разработки предлагаемой конструкции, расчета и выбора оптимальных конструктивных параметров нагревателя, частоты тока источника питания, схемы трехфазной обмотки индуктора решается ряд задач, связанных с исследованием электромагнитных и тепловых полей и созданием на этой основе алгоритма и методики расчета конструктивных и режимных параметров многослойной индукционной системы. Основным видом нелинейной среды являются ферромагнитные участки магнитной цепи и стальные конструктивные элементы, для которых связь между индукцией и напряженностью магнитного поля определяется магнитными свойствами среды. Известная неопределенность зависимости связана с проявлением гистерезиса и наличием частных циклов намагничивания, в связи с чем вектор индукции зависит не только от напряженности магнитного поля, но и от предыдущего ее изменения в данной точке, а также от начальной намагниченности. При решении нелинейных уравнений электромагнитного поля основную кривую намагничивания аппроксимируют аналитическими выражениями, которые, с одной стороны, должны достаточно точно описывать эту кривую, а с другой - допускать интегрирование системы уравнений поля в удобном для расчетов виде. Наибольшее распространение получила параболическая зависимость . Математическая модель процесса индукционного нагрева многослойной системы «индуктор - внешняя труба - вращающийся ротор - жидкость» может быть представлена нелинейной взаимосвязанной системой уравнений Максвелла [11] и Фурье [12] соответственно для электромагнитного и теплового полей с соответствующими краевыми условиями: ; ; (1) ; (2) ; (3) . (4) Здесь , - векторы напряженности магнитного поля, магнитной индукции и напряженности электрического поля; - вектор плотности первичного тока; - вектор плотности индуцированного тока; - время; ,, - соответственно температурные поля в стенке немагнитной внешней трубы, в стенке ротора, в нагреваемом потоке жидкости внутри ротора; , , , , , - удельные значения теплоемкости и плотности материалов трубы, ротора и жидкости соответственно; - коэффициенты теплопроводности материала внешней трубы, ротора и нагреваемой жидкости соответственно; - скорость потока. Объемная плотность внутренних источников тепла, индуцируемых в стенках труб, определяется дивергенцией вектора Пойнтинга . Система уравнений (1) - (4) дополняется граничными условиями для электромагнитной и тепловой задач. Для электромагнитной задачи используются условия равенства функции нулю на бесконечно удаленной границе и условие симметрии на осевой линии , которое заключается в равенстве нулю производной от функции. Численный расчет электромагнитных полей в сложной составной структуре тел, содержащей ферромагнитные участки магнитной цепи и стальные конструктивные элементы, производился с помощью программного комплекса ELCUT 5.9 Professional [13]. Для расчета интегральных параметров индукционной системы использовалась двумерная постановка электромагнитной задачи в форме системы дифференциальных уравнений. Перейдя от системы уравнений Максвелла к формулировке с использованием векторного магнитного потенциала, можно записать дифференциальное уравнение с соответствующими граничными условиями ; . Здесь - векторный магнитный потенциал; - относительная магнитная проницаемость материала по осям x, y; j - мнимая единица; - круговая частота тока; g - удельная электрическая проводимость; L - граница расчетной области (линия удаленной поверхности); - плотность стороннего тока. На базе приведенных систем уравнений строится конечно-элементная формулировка для плоской двумерной области [14]. Граница раздела магнитных сред описывается системой соотношений: Последнее выражение учитывает скачкообразное изменение вектора напряженности на границе раздела сред. Кроме условий сопряжения для получения однозначного решения уравнений Максвелла в форме напряженности электрического поля и напряженности магнитного поля в области с границей S необходимо задать: - уравнения поверхностей, отделяющих друг от друга среды i и j, ; - начальные величины E(x,y,z), H(x,y,z) в момент времени в произвольной точке исследуемого объема с границей S; - касательные составляющие вектора или в произвольной точке поверхности в произвольном временном интервале от до t или распределения полей и вне исследуемого объема Q; - функциональные зависимости магнитной проницаемости μ и удельной проводимости γ от пространственных координат или от напряженности соответствующего поля. Решение задачи электромагнитного поля достигается использованием векторного магнитного потенциала {A} и скалярного электрического потенциала V, которые выражаются следующим образом: ; . Чтобы функция была определена, нужно определить значение ее дивергенции. Для этого добавляется условие, которое называется калибровкой Кулона: В результате получим следующую систему уравнений: ; Используя соотношение при μ=const получим уравнение Уравнение Пуассона дополняется граничными условиями Дирихле и Неймана на различных участках границы: на ; на . Такое упрощение условий задачи объясняется тем, что дальнейший переход к конечно-элементной формулировке намного облегчается для линейной задачи. Для учета нелинейной зависимости в ферромагнитных областях используется итерационный алгоритм многократного решения результирующей системы уравнений. В начальной стадии расчета задается значение μ=const по всей области ферромагнитных макроэлементов, затем вычисляются распределенные параметры поля, что позволяет на следующей стадии расчета корректировать μ внутри каждого конечного элемента в зависимости от значения напряженности магнитного поля в данной области. Определение магнитной проницаемости производится с помощью введения в программу расчета полинома, аппроксимирующего кривую намагничивания. Рассмотренная модель электромагнитной задачи реализована в виде пакета программ, предназначенного для решения двухмерных задач. Численный расчет электромагнитных и тепловых полей в сложной составной структуре тел, содержащей ферромагнитные участки магнитной цепи, стальные конструктивные элементы и ферромагнитную загрузку, производился с помощью программного комплекса ELCUT 5.9 Professional - интегрированной диалоговой системы программ, позволяющей решать двумерные краевые задачи математической физики, описываемые эллиптическими дифференциальными уравнениями в частных производных относительно скалярной или однокомпонентной векторной функции. Рис. 2. Распределение плотности тока в элементах системы Исходные данные Параметр Значение Диаметр внешней трубы, мм 169 Материал внешней трубы Нержавеющая сталь Толщина стенки внешней трубы, мм 4,5 Диаметр цилиндрического ротора внешний, мм 158 Толщина стенки цилиндра, мм 12 Диаметр цилиндрического ротора внутренний, мм 134 Материал цилиндрического ротора Сталь ферромагнитная Размер паза (высота×ширина), мм 5×7 Размер зуба (высота×ширина), мм 5×7 Обмотка к.з., алюминий, мм Шинка (4×6) Толщина электрической изоляции, мм 0,5 Количество пазов, шт. 36 Зазор между внешней трубой и цилиндрическим ротором, мм 2 Зазор между внешней трубой и статором, мм 2 Диаметр вала, мм 30 Диаметр статора внутренний, мм 171 Диаметр статора внешний, мм 231 Материал статора Электротехническая сталь Толщина набора магнитопровода, мм 20 Длина статора, мм 400 Размер паза (высота×ширина), мм 17х10 Размер зуба (высота×ширина), мм 17х8 Частота тока индуктора, Гц 50 Кол-во пазов, шт. 36 В качестве исходных данных для решения электромагнитной задачи вводятся свойства сред, источники поля, распределенные и сосредоточенные токи, кривые намагничивания ферромагнитных материалов, граничные условия и др. Основными расчетными параметрами являются изменяющиеся во времени магнитный потенциал, магнитная индукция, напряженность поля, плотности токов, удельная тепловая мощность, вращающий момент. Исходными данными для решения электромагнитной задачи являются: - конструктивные параметры нагревателя; - электрофизические характеристики материала труб, табличные значения зависимостей относительной магнитной проницаемости материала труб от напряженности магнитного поля, параметры индуктирующей катушки, размеры и характеристики магнитопровода; - энергетические параметры; - напряжение питания, частота тока. Численные эксперименты выполнены для устройства с параметрами, приведенными в таблице. На рис. 2 показано распределение плотности тока в элементах электромагнитной системы. На рис. 3, 4, 5, 6 представлены графики момента ротора, плотности тока и тепловыделения в стенках ротора и внешней трубы. Рис. 3. Момент ротора По результатам расчета построены графики значений плотности тока и тепловыделения в элементах индукционной системы. Полученные результаты решения электромагнитной задачи позволяют рассчитать температурные распределения в расчетной области. Результаты расчета для приведенного в работе примера показали следующее распределение мощности: стенка внешней трубы - 18 %, стенка ротора - 55 %, мощность вращения ротора - 27 %, момент ротора - 152,28 Нм. Рис. 4. Мощность тепловыделения в роторе Рис. 5. Мощность тепловыделения в трубе из нержавеющей стали Плотность тока в зубе Плотность тока в пазу Тепловыделение в пазу Тепловыделение в нержавеющей трубе Рис. 6. Плотность тока и тепловыделения в элементах системы Приведенные результаты существенно зависят от частоты источника питания, геометрических параметров системы и электрофизических характеристик материалов. Полученные результаты позволяют выполнить анализ зависимости электромагнитных параметров индуктора - индукции магнитопровода, плотности тока в проводниках, объемной мощности тепловыделения в металлических элементах конструкции, вращающего момента - и других характеристик от геометрических и электрофизических характеристик системы, частоты и напряжения источника питания. На основании всестороннего анализа полученных зависимостей устанавливается закономерность и степень влияния электрических и геометрических параметров на энергетические характеристики системы. Полученные закономерности далее могут быть использованы при решении задачи выбора оптимальных параметров системы, обеспечивающей требуемые эксплуатационные характеристики в стационарных режимах работы установки.

About the authors

Alexander A Bazarov

Samara State Technical University

(Dr. Sci. (Techn.)), Professor. 244, Molodogvardeyskaya st., Samara, 443100, Russian Federation

Alexander I Danilushkin

Samara State Technical University

(Dr. Sci. (Techn.)), Professor. 244, Molodogvardeyskaya st., Samara, 443100, Russian Federation

Vasily A Danilushkin

Samara State Technical University

(Ph.D. (Techn.)), Associate Professor. 244, Molodogvardeyskaya st., Samara, 443100, Russian Federation

Ivan V Vasilyev

Samara State Technical University

Postgraduate Student 244, Molodogvardeyskaya st., Samara, 443100, Russian Federation

References

Supplementary files