ВЕСТН. САМАР. ГОС. ТЕХН. УН-ТА. СЕР. ТЕХНИЧЕСКИЕ НАУКИ. 2015. № 4 (48)
УДК 621.365.53
ИНДУКЦИОННАЯ СИСТЕМА ДЛЯ СНЯТИЯ ТЕРМОНАПРЯЖЕНИЙ В ТРУБЕ ПОСЛЕ СВАРКИ
А.А. Базаров, Д.Ю. Чигирев
Самарский государственный технический университет Россия, 443100, г. Самара, ул. Молодогвардейская, 244
Рассматривается задача расчета параметров индукционного нагревателя трубы для снятия термонапряжений после сварки. Сварка многих материалов предусматривает предварительный подогрев для исключения появления недопустимых термонапряжений. Несмотря на такие мероприятия, остается проблема снятия остаточных напряжений. При разработке конструкции и режима работы нагревателя требуется учет работы на двух этапах для уменьшения массовых и габаритных показателей и лучшего согласования с источником питания. Использование жесткой конструкции позволяет обеспечить независимость характеристик индуктора от качества изготовления теплоизолирующего слоя между трубой и нагревателем. Предложена разъемная конструкция индуктора, обеспечивающая удобство монтажа и высокую надежность. Выбранные частота и мощность индуктора позволяют решить проблему согласования с источником питания без трансформатора. В конструкции индуктора применен ферритовый магнитопро-вод, что дало возможность использовать двухсекционную петлевую обмотку.
Ключевые слова: индуктор, термонапряжения, магнитопровод, проектирование, температура, конечные элементы.
Эксплуатация оборудования тепловых электростанций происходит в очень тяжелых условиях, предъявляющих повышенные требования к надежности. Одним из таких элементов является паропровод, который используется на промышленной ТЭЦ. Передача большого количества пара от котла к турбогенератору осуществляется под большим давлением (140 кг/см2). Главный паропровод в рассматриваемой ситуации имеет диаметр 325 мм, толщину стенки 38 мм и изготавливается из стали 12Х1МФ.
Для предотвращения возникновения внутренних напряжений при изготовлении и ремонте требуется производить предварительный подогрев трубы до температуры 300 °С [1]. После завершения сварочных работ необходимо прогреть участок паропровода до 720 °С за 1,5 часа, затем выдерживать зону сварного шва при этой температуре в течение 3 часов. После этого зону нагрева следует медленно остужать, понижая мощность индуктора. Снижение температуры до 300 °С в течение 1 часа осуществляется в условиях тепловой изоляции. По окончании этого этапа происходит естественное остывание паропровода до температуры окружающей среды на открытом воздухе.
Следует отметить, что при положительной температуре окружающей среды термообработку следует проводить не позже чем через 72 часа после сварки, а если температура окружающей среды отрицательная, то термообработку следует проводить незамедлительно. Так как паропровод нагревают до температуры,
Александр Александрович Базаров (д.т.н.), доцент кафедры «Электроснабжение промышленных предприятий».
Дмитрий Юрьевич Чигирев, аспирант.
которая находится ниже критической, то в сварном шве не будет происходить структурных изменений. Температура термообработки контролируется с датчиков температуры.
Условия прокладки паропровода предъявляют особые требования к массога-баритным показателям и к системе электроснабжения устройства нагрева. Использование индукционных систем сопровождается повышенной частотой напряжения, низким коэффициентом мощности, что вызывает увеличение потерь напряжения и мощности в сети. С учетом в некоторых случаях большой высоты прокладки паропровода и сложности доступа решение задачи проектирования индукционной системы должно быть ориентировано на достижение таких показателей, как минимальная масса комплекса оборудования, максимум автономности от инженерных систем, высокая точность регулирования температуры.
Конструкция индуктора должна соответствовать кроме быстрой установки требованию охлаждения обмотки. Поэтому простые решения, заключающиеся в намотке катушки гибким изолированным проводом, в данном случае непригодны. Вместо классического соленоидального индуктора [2] более приемлемым вариантом является изготовление его в виде разъемной конструкции из двух половинок, содержащих петлевые обмотки [3]. Для устранения компенсирующего эффекта для встречных токов в расположенных рядом проводниках предусмотрено использование магнитопровода. В качестве материала возможно применение феррита МН4000, исключающего появление токов проводимости. Сравнительно невысокое значение индукции насыщения, примерно равное 0,35 Тл, тем не менее позволяет решить поставленные задачи. Конструкция половинки индуктора без корпуса и теплоизоляции показана на рис. 1. Вторая половина имеет зеркальное отображение и обеспечивает замыкание наведенных токов в виде колец в поверхностном слое трубы. Для простоты не показано наличие каналов водяного охлаждения в проводниках обмотки.
Рис. 1. Фрагмент разъемного петлевого индуктора: 1 - магнитопровод; 2 - обмотка
Для проектирования индукционного нагревателя применены программные пакеты Сош8о1 и Б1си1;, позволяющие проводить расчеты электромагнитных и тепловых процессов как раздельные, так и связанные.
Общая формулировка процесса индукционного нагрева описывается нелинейной системой уравнений Максвелла и Фурье для электромагнитного и теплового полей соответственно:
гоИ = -5 +
rotE = -5 +
8 О
"ВР
8В
йыИ ; йУВ = 0;
(1)
8т "8Г
йТ) - йщ
Е • И
(2)
Здесь
векторы напряженностей магнитного и электрического
полей, электрической и магнитной индукции;
5 = — Е - плотность тока;
Р
1 - коэффициент теплопроводности;
Т - температура.
Объемная плотность внутренних источников тепла, индуцируемых в металле, определяется дивергенцией вектора Пойтинга П = — й1у\Е • И ].
Система уравнений (1), (2) дополняется краевыми условиями, которые будут рассмотрены ниже.
В данной работе исходной математической моделью является дифференциальное уравнение вида
8Т (г, х, t)
8t
= а(Т)
8 2Т (г, х, t) + — 8Т (г, х, t) + _8МГ1Х0
8г2 г 8г 8х
с начальными и граничными условиями
Т (г, х,0) = То = Тс;
Т (г, х, t) х=—х=Тс; Тх, t} х=х=Тс; ■ = Н (Т (г , х, t) -
+ К (г, х, t), (3)
(4)
^^ = Н (Т(Г1, х, 0-Тс); = Н2 (Т(Г2, х, t)-Тс);
8г 8г
где г— и Г2 - внешний и внутренний радиусы трубы;
X,-X - границы расчетной области трубы;
с - удельная теплоемкость;
у - плотность материала трубы;
Н - приведенный коэффициент теплообмена, Н = а;
1
а - коэффициент температуропроводности, а = —;
су
(г, х, t) - функция распределения внутренних источников тепла. Численное решение задачи основывается на минимизации энергетического функционала вида
с =
Я
1»,
2
т
сг
+
дТ_ дх
2 \
дТ
+ с--жг
л
йтйх + — 2
2 |к(г2 - 2ТСТ}и. (5)
Линейный вид уравнения (5) не соответствует зависимости параметров процесса от температуры. Для преодоления этой сложности используется итерационный процесс решения системы уравнений, содержащих функции коэффициентов теплопроводности и теплоемкости от температуры. Для корректного учета влияния температуры на свойства материала размеры конечных элементов, особенно в области тепловыделения, принимаются достаточно малыми.
300
700
600
500
о" га о.
I
га
о- 400
£
Ш
300
200
100
Г
и 2
1
в
-
1000
2000
3000
4000 5000 время, с
6000
7000 3000
Рис. 2. Временные диаграммы температур участка трубы под индуктором:
1 - при мощности индуктора 37,5 кВт; 2- при мощности индуктора 15 кВт;
3 - при мощности индуктора 10,5 кВт
Целью тепловых расчетов является выбор необходимой мощности индуктора, обеспечивающей требуемую скорость нагрева. Создание модели замкнутой системы управления нагревом позволило определить приемлемое значение мощности. Приведенные на рис. 2 диаграммы температур при мощностях тепловыделения от 10,5 до 37,5 кВт показали возможность достижения температуры за разное время. Минимальная мощность, равная 10,5 кВт, обеспечивает нагрев за заданное время, равное 1,5 часа. В дальнейшем, как показано на рис. 3, вся потребляемая мощность расходуется на компенсацию потерь тепла, рассеиваемых в окружающую среду с боковой поверхности трубы. Также учитывается тепло-отвод через теплоизоляцию к индуктору. Соотношение мощностей индуктора в режимах разогрева и стабилизации температуры не превышает двух, что говорит о высоком коэффициенте использования. Кроме того, эта сравнительно небольшая мощность не вызывает проблем с электроснабжением практически в
2
5
к
любом месте производственного помещения. Для достижения требуемых показателей тепловой изоляции по ограничению потерь тепла и механической прочности целесообразно применение современных материалов.
Анализ электромагнитных полей проводится на основе уравнений (1) Максвелла. Учитывая, что длины дуги витка и перемычки между соседними витками отличаются почти в 15 раз, можно с допустимой погрешностью расчетов перейти от трехмерной к двумерной постановке задачи. Такое упрощение оборачивается гораздо более существенной экономией вычислительных ресурсов.
хЮ
2.5
£
"Р Ш
з
о Е
1.5
0,5
1
2
3
1000
2000
3000
4000 время, с
5000
6000
7000
3000
Рис. 3. Временные диаграммы изменения мощности индуктора: 1 - при мощности индуктора 37,5 кВт; 2 - при мощности индуктора 15 кВт; 3 - при мощности индуктора 10,5 кВт
Для практического использования более приспособлена формулировка задачи в терминах векторного магнитного потенциала:
( 1
rot
1 -
— rot A
m a
dA -+ c- = J o;
(6)
divA = 0;
где
rot A = B;
A - векторный магнитный потенциал;
B - магнитная индукция;
J o - плотность наведенных токов;
s - электрическая проводимость;
m a = m om - абсолютная магнитная проницаемость.
Переходя к конечно-элементной формулировке, нелинейную электромагнитную задачу записывают в виде уравнения Пуассона в двумерной осесимметрич-ной области V (г, х) вида
- > А(г, х)- J 0 = 0, (7)
где ш = 2%/ - круговая частота.
Для определенности решения задачи на границе расчетной области векторный магнитный потенциал принимается равным нулю. Линии потока в условиях геометрической симметрии принимаются перпендикулярными плоскостям симметрии:
. I СА I
А^ = 0; —к-5 = 0; (8)
где 51 - удаленная граница расчетной области Q> 5 2 - граница плоскостей симметрии.
Рассмотренная постановка электромагнитной задачи позволяет решить задачи в нелинейной постановке. Нужно отметить, что наличие двух материалов, имеющих существенные отличия в свойствах, вносит некоторые проблемы. Феррит имеет меньшее значение индукции насыщения, что требует больших размеров магнитопровода по сравнению с обычным электротехническим железом. Второй особенностью является прямоугольная петля гистерезиса, что сказывается на более резком переходе от почти вертикального участка к участку с небольшим наклоном. Такой характер нелинейности плохо сказывается на сходимости вычислительного процесса. Для ускорения вычислительного процесса целесообразно переход делать более плавным.
При расчете параметров индуктора, изображенного на рис. 1, были исследованы варианты с разной частотой тока. Результаты расчетов приведены в табл.1. Исходя из располагаемой длины индуктора, равной 0,15 м, были приняты: ширина витка 0,01 м; ширина зубца магнитопровода 0,01 м; число витков w = 8; общая длина индуктора Ь = 0,15 м.
Таблица 1
Параметры петлевого разъемного индуктора с магнитопроводом
F, кГц Р2, кВт ДРи, кВт Ри, кВт Ц,, В 8, кВА I, А 008 ф
2 400 10,55 11,77 22,32 170 244,3 1437 0,091
4 000 10,5 9,96 20,49 197 205,9 1045 0,1
10 000 10,5 7,66 18,16 275 179,3 652 0,101
Наиболее предпочтительным вариантом оказался третий, использующий частоту 10 кГц. Напряжение на индукторе равно 275 В, что сопоставимо с выходным напряжением стандартного преобразователя и не требует применения согласующего трансформатора. Кроме частоты на увеличение напряжения на индукторе влияет наличие магнитопровода [4]. Возможно еще большее приближение к параметрам источника питания, но здесь есть один минус: по мере нагрева трубы увеличивается ее удельное сопротивление, что приводит к снижению мощности, выделяющейся в загрузке. Для компенсации этого явления можно
дх
д А(г, х)
т а (г, х)
дх
д + —
дг
1 1 дг А(г, х)
т а (г, х) г дг
1
поднять напряжение на выходе преобразователя. В случае равенства расчетного напряжения для холодной загрузки и выходного напряжения преобразователя эта возможность утрачивается.
Для сравнения приведены параметры соленоидального индуктора без магни-топровода с числом витков, равным 15 (табл. 2). Учитывая малый зазор между индуктором и трубой, равный 10 мм, энергетические параметры этого варианта можно считать более предпочтительными. Однако сложность изготовления, особенно обеспечения водяного охлаждения, делает такую конструкцию малопригодной.
Таблица 2
Параметры соленоидального индуктора
F, кГц Р2, кВт ЛРи, кВт Ри, кВт Ц,, В 8, кВА I, А сОБ ф
2 400 10,55 0,53 11,08 191 38,37 201 0,289
4 000 10,52 0,366 10,88 227 32,4 143 0,336
10 000 10,6 0,33 10,93 240 24,7 103 0,442
Особенностью распределения мощности тепловыделения в трубе является эффект полос, вызванный протеканием встречно направленных токов в соседних витках. Это выражается в чередовании максимальной мощности и провалов (рис. 4).
Рис. 4. Распределение удельной мощности тепловыделения в трубе под индуктором
При высоких скоростях нагрева такое явление могло бы привести к значительным градиентам температуры в рабочей зоне. В условиях низкой интенсивности нагрева, малой ширины полос и сравнительно высокого коэффициента теплопроводности наблюдаемый эффект практически исчезает, что показано на распределении температуры в трубе в процессе разогрева (рис. 5).
О О а 0.4 0.6 0.8 1 1.2 1.4 1.0 1.8 2
длина, м
Рис. 5. Распределение температуры под индуктором при выходе на режим стабилизации температуры
Полученное решение для индукционного нагревателя позволяет ориентироваться на комплекс оборудования: транзисторный источник питания СЧВ-35В мощностью 35 кВт, массой 31 кг; конденсатор ЭСПВ-0,4-10 мощностью 650 кВАр, массой 25 кг; чиллер для охлаждения воды ЧА-8,7 производительностью по холоду 8,7 кВт, массой 150 кг. За исключением системы охлаждения воды оборудование является очень компактным и легко транспортируемым. При использовании автомобильных гидроподъемников эта проблема полностью решается.
Выполненное исследование по применению современных технических решений в области проектирования индукционных нагревателей показало высокую эффективность и точность расчетов с помощью существующего аппаратного и программного обеспечения и реализации на базе современных материалов и оборудования.
БИБЛИОГРАФИЧЕСКИЙ СПИСОК
1. Лившиц Л. С., Хакимов А.Н. Металловедение сварки и термическая обработка сварных соединений. 2-е изд. - 1989. - 167 с.
2. Йон Инге Асперхайм, Лейф Маркегорд. Оптимизация процесса нормализации сварного шва с помощью 2Б-моделирования // Индукционный нагрев. - 2009. - № 4. - С. 21-25.
3. Слухоцкий А.Е., Рыскин С.Е. Индукторы для индукционного нагрева. - Л.: Энергия, 1974. -264 с.
4. Филиппов К.П., Булатова Л.Ш. Индукторы для нагрева осесимметричных изделий большого диаметра // Индукционный нагрев. - 2012. - № 2. - С. 37-44.
Статья поступила в редакцию 25 сентября 2015 г.
INDUCTION SYSTEM FOR REMOVAL OF THERMAL STRESSES IN THE PIPE AFTER WELDING
A.A. Bazarov, D. Y. Chigirev
Samara State Technical University
244, Molodogvardeyskaya st., Samara, 443100, Russian Federation
The problem of calculating the parameters of a pipe induction heaterto removethermal stresses after welding is discussed. In welding, many materials are supposed to be preheated to avoid the development of harmful thermal stresses.Despite such arrangements, there is a problem of removing residual stresses. In the heater's design and operation mode, the consideration of the two-stage development is required for decreasingtheweight and dimensions, and for better coordination with the power supply. Using the rigid structure allowsensuring the inductor parameters independence from the production quality of the insulating layer between the pipe and heater. The sectional structure of the inductor ensures easy installation, and reliability. The selected frequency and power of the inductor can solve the problem of the coordination with the power supply without a transformer. Magnetic ferrite isapplied in the design of the inductor, which made it possible using a two-section lap winding.
Keywords: inductor, thermal stresses, magnetic, design, temperature, finite elements.
Alexander A. Bazarov (Dr. Sci. (Techn.)), Professor. Dmitry Y. Chigirev, Postgraduate Student.