CC BY
10ВЕСТН. САМАР. ГОС. ТЕХН. УН-ТА. СЕР. ТЕХНИЧЕСКИЕ НАУКИ. 2020. Том 28. № 3_
doi: 10.14498/tech.2020.3.9
УДК 517.977.56, 519.6, 621.785
ОПТИМАЛЬНОЕ ПРОЕКТИРОВАНИЕ ИНДУКЦИОННОГО НАГРЕВАТЕЛЯ ДЛЯ ПОВЕРХНОСТНОЙ ЗАКАЛКИ В УСЛОВИЯХ ИНТЕРВАЛЬНОЙ НЕОПРЕДЕЛЕННОСТИ ПАРАМЕТРОВ*
А.В. Попов
Самарский государственный технический университет Россия, 443100, г. Самара, ул. Молодогвардейская, 244
Аннотация. Работа посвящена оптимальному проектированию индукционной нагревательной установки для поверхностной закалки стальных цилиндрических деталей. Стадия нагрева поверхностной индукционной закалки рассматривается в качестве объекта с распределенными параметрами, для которого неизвестными оптимизируемыми параметрами являются конструктивные характеристики индуктора. Поскольку в реальных промышленных условиях часто информация об объекте известна не полностью, основные параметры могут задаваться с точностью до диапазона значений их возможных изменений. В работе формулируется задача оптимального проектирования нагревателя в условиях интервальной неопределенности начальной температуры обрабатываемой детали и комплексного коэффициента теплоотдачи с поверхности (учитывающего лучистый и конвективный теплообмен). Сформулированная задача решается с помощью альтернансного метода параметрической оптимизации систем с распределенными параметрами на базе численной модели процесса нагрева, разработанной в конечно-элементном программном пакете Altair FLUX.
Ключевые слова: поверхностная закалка, индукционный нагрев, оптимальное проектирование, альтернансный метод, интервальная неопределенность, численная модель, взаимосвязанные электромагнитные и тепловые поля, Altair FLUX.
Введение
Индукционная закалка в настоящее время широко используется в технологических процессах термической обработки металлов. Конечной целью индукционной закалки является получение определенных изменений в микроструктуре обрабатываемого материала с целью его упрочнения, повышения сопротивления изнашиванию или усталостному разрушению. Часто такие изменения структуры металла требуется получить по всему объему обрабатываемой детали, для чего применяется сквозная индукционная закалка. В случаях, когда необходимо упрочнить только поверхностный слой изделия без изменения свойств сердцевины, применяется поверхностная индукционная закалка [1].
В работе рассматривается процесс поверхностной индукционной закалки стальных заготовок, который в самом общем случае состоит из трех стадий: нагрева поверхностного слоя детали выше критической температуры аустенит-ного превращения, выдерживания для получения однородной аустенитной
* Работа выполнена при финансовой поддержке РФФИ в рамках научных проектов №20-08-00240 и №19-08-00232.
Попов Антон Валерьевич, аспирант, старший преподаватель кафедры «Управление и системный анализ теплоэнергетических и социотехнических комплексов».
структуры и резкого охлаждения нагретого слоя ниже температуры образования мартенсита [2]. Основной целью поверхностной закалки является получение в обрабатываемом слое однородной мартенситной структуры равномерной глубины, которая зависит от равномерности нагрева границы закаливаемого слоя на первой стадии. Поэтому главной задачей стадии нагрева процесса поверхностной индукционной закалки можно считать получение максимальной равномерности температурного поля в закаливаемом слое заготовки.
Во многих случаях для обеспечения требуемого температурного профиля в упрочняемом слое оказывается недостаточным использование типовых индукционных нагревательных установок, поэтому актуальной становится задача оптимизации конструктивных характеристик индуктора. Кроме того, в реальных промышленных условиях точные значения основных характеристик процесса индукционной закалки практически всегда неизвестны, однако возможно определить некоторые диапазоны изменения их допустимых значений. В связи с этим существенный интерес приобретает задача оптимального проектирования индукторов для стадии нагрева процесса поверхностного упрочнения стальных деталей в условиях интервальной неопределенности основных параметров, которая решается в данной работе.
Постановка задачи оптимального проектирования в условиях
интервальной неопределенности
В работе рассматривается задача оптимального проектирования двух-виткового индукционного нагревателя для реализации стадии нагрева поверхностной индукционной закалки стальных заготовок цилиндрической формы, неизвестными оптимизируемыми параметрами которой могут являться конструктивные параметры индукционной нагревательной установки и характеристики источника питания, объединенные в вектор Р [3]. В данном случае компонентами этого вектора являются: размер витка индуктора квадратного сечения (Р\ х Р\), расстояние между заготовкой и индуктором (р2), расстояние между витками (р3) и ток источника питания (Р4 = I) (рис. 1).
Рис. 1. Геометрия индукционной системы
На искомые параметры накладываются следующие ограничения:
^тт ~ ~ plmax, p2min ~ p2 ~ p2max, p3min ~ p3 ~ ^тах^тт ~I ~ Imax' (1)
Рассматриваемая стадия нагрева представляет собой объект с распределенными параметрами, функцией состояния которого является пространственно-временное распределение температуры по объему нагреваемой детали Т (г, I, ¿), описываемое двумерной системой нелинейных уравнений Максвелла и Фурье для электромагнитного и температурного полей [4-6]:
rotH = a(T) E,
rotE = _—B, dt
divB = 0, divE = 0,
'{T W)
dl
Л(Т)
ST {r, l, t) _ 1 d_
dt r dr
д T { r, l, t) >
Я(Т )
,5 T {r,l, t)
dl
a
{T)
v dr у
( dH {r, l, t)^2
dr
dH(0,l,t) _o.
dr
= 0;H (R,l,t)=HL;H (r, 0,t)=HR1; H (r, L,t) = H^; = 0; (T^ = a(T)(T(R,l,tК); ( f-^ = a(T )(T (r, 0,t )-Tfl); (T = a(T )(T (r, L,t )-Te).
(2)
(3)
(4)
(5)
(6)
с начальными и граничными условиями:
T(r,l,t)=T(r,l,0)=T(r,l) = T = const,l e[0;L],r e[0;R]; (7)
(8) (9)
(10)
Здесь Н - вектор напряженности магнитного поля; В - вектор магнитной индукции; Е - вектор напряженности электрического поля; ст(Т) - электрическая проводимость; Т - температура;
с(Т),у(Т),1(Т) - соответственно удельная теплоемкость, плотность и коэффициент теплопроводности нагреваемого металла; - время;
ге[0;^],Iе[0;Х]- соответственно радиальная и продольная пространственная координаты;
T - температура окружающей среды; a (T) - коэффициент теплоотдачи в окружающую среду. В качестве оценки точности достижения заданного равномерного распределения температуры по границе упрочняемого слоя глубиной r* предлагается выбрать максимальную величину абсолютного отклонения температуры
T(R — r*,l,t0) вдоль этой границы в конце стадии нагрева длительностью t° от
заданного значения T* [7]. Обычно при закаливании сталей требуется нагреть поверхностный слой изделия до температур выше критических, необходимых для образования аустенита, которые в типовых случаях составляют 850-900 °С [1-2].
В соответствии с вышеизложенным в работе рассматривается следующий критерий оптимизации конструкции индукционной нагревательной установки:
I (P) = max T (R — r*,l,t0 ,P ) —T*\^ min (11)
v ' le[0,i]l V ' I P
при T* = 900°С.
Поскольку в данной задаче информация об объекте оказывается неполной и некоторые основные параметры задаются с точностью до возможного диапазона их изменения, минимизируемый критерий (11) принимает следующий вид:
I (p) = max
V ) ßeQ
max
¡e[0;L]
T (R - r,¡,t0 P,ß)-T*
►min, (12)
p
где Р - вектор неполностью определенных параметров [4-5].
В частном случае рассматриваемой в работе задачи с наличием интервальной неопределенности параметров стадии нагрева процесса поверхностной закалки к вектору Р относятся: начальная температура закаливаемой детали T0 и обобщенный коэффициент теплообмена с поверхности а (конвективного или лучистого), которые заданы с точностью до известных интервалов их изменения
[T0mm';T0max ] и [^n^max ] соответственно, т. е.
Р = (T0 >а); П : T0min ^ T0 ^ T0mx,amn ^a^amax' (13)
Таким образом, задача оптимального проектирования индуктора для стадии нагрева поверхностной индукционной закалки, описываемой уравнениями (2)-(6) с краевыми условиями (7)-(10), сводится к поиску таких значений вектора оптимизируемых параметров P = (px,p2,p3,I), стесненных ограничениями (1), которые обеспечивают нагрев границы закаливаемого поверхностного слоя заготовки толщиной r* до требуемой температуры T* за заданное время t° при минимальном значении критерия оптимальности (12) в условиях наличия интервальной неопределенности начальной температуры детали и коэффициента теплообмена (13).
Сформулированная задача оптимального проектирования может быть решена с помощью альтернансного метода параметрической оптимизации систем с распределенными параметрами.
Методика решения задачи оптимального проектирования с помощью альтернансного метода
Предлагаемая методика решения сформулированной задачи оптимального проектирования базируется на использовании альтернансного метода параметрической оптимизации систем с распределенными параметрами [4-5, 8-9]. Данный метод опирается на установленные универсальные свойства пространственных распределений температуры по объему заготовки в конце оптимальных процессов индукционного нагрева металла, подобные известным в математике свойствам наилучших приближений заданных функций к нулю. На этой основе производится процедура точной редукции исходной задачи оптимизации к решению трансцендентных систем уравнений, замкнутых относительно всех искомых параметров.
Оптимальным значениям параметров вектора Р° =(p10^^,...,p0), являющихся решением рассматриваемой задачи оптимального проектирования при наличии полной информации об объекте, отвечает пространственное распределение температур Т(я — г ,/, г0,Р0), которое должно соответствовать минимальному значению критерия оптимальности (11) [3]. Основное свойство результирующего температурного распределения Т(я — г ,0,Р0) состоит в том, что
число К точек по объему нагреваемого тела (в данном случае вдоль рассматриваемого продольного сечения при г=г ), в которых достигаются предельно допустимые абсолютные отклонения конечной температуры Т(я — г ,/, г0,Р0) от
требуемого значения Т , равные £0 , всегда оказывается не меньше числа N искомых оптимизируемых параметров [4-5]. При этом максимально допустимые отклонения Т (я — г ,/,г0, Р0) — Т оказываются знакочередующимися в точках
/0,] = 1,К, где 0 < </° <... <<Ь , т. е. обладают альтернансным свойством. Причем для температурного отклонения в точках с неизвестными координатами ,т = 1,К < К выполняются необходимые условия экстремума
дТ (я—г*,/, г0, Р0)
—*-'-=0.
д/
Вышеуказанные утверждения позволяют составить следующую систему уравнений:
Т (я—г* 0, Р°)—Т =(—1) V0, ]=1, К, Р° =( А°, Р20^^) ,^ = ±1;
п п п \ N если < £0 < £т„4
0 < /0 < /0 <... < /0 < Ь; К = •!_\ т'„ 0 т'„ ; (14)
~ — N)
'0 £т'„
.0 00 '
(Я — г\?т/Р) _
-А-¿ = 0; т = 1, К < К.
дТ (я—г , /т, г0, Ри
аТ
Здесь г^ - предельно достижимая точность нагрева в задачах с N оптимизируемыми параметрами.
Система уравнений (16) оказывается замкнутой относительно всех неизвестных, включая неизвестные параметры р1,р2°,•••,, координаты ¡], ] = 1,к
и значение г^, поэтому ее решение и будет являться решением рассматриваемой задачи оптимального проектирования.
Поскольку в работе рассматривается вектор Р, содержащий четыре неизвестных конструктивных параметра индукционной нагревательной установки, оказывается, что конечное температурное распределение вдоль границы закаливаемого слоя толщиной г* имеет пять точек с максимальным отклонением температуры от заданного значения (рис. 2).
Рис. 2. Формы конечных температурных распределений для предельно
« —( 4)
достижимой точности нагрева г^ ¡^
Очевидно, что данному температурному распределению соответствует максимально достижимая точность нагрева г*4).
Система уравнений (14) в данном случае принимает следующий вид [3, 7]:
Г (К - ЛI] / )-Г = (-1)] гЦ ] = 1,5, Р =( р0, р0, р0, р0);
о < ¡0 < ¡0 < ¡0 < ¡0 < ¡0 < ц
дт(к - г", ¡0, Г0, Р0 ) _
—*-= 0; т = 0,4.
(15)
д¡
Рассмотрим теперь задачу оптимизации конструкции индукционного нагревателя в случае, когда информация об объекте неполная, т. е. основные параметры стадии нагрева заданы с точностью до возможных интервалов их допустимых значений.
Пусть согласно (13) имеется некоторый векторный параметр Р, который обозначает совокупность любых фиксированных значений не полностью определенных характеристик индукционной системы.
Сложность данной задачи по сравнению с задачей с полной информацией об объекте заключается в необходимости рассматривать максимум температурных
отклонений на множестве допустимых значений вектора Р вместо пространственной координаты /. В [4] показано, что основные свойства конечного температурного поля, описанные в системе уравнений (16), распространяются и на рассматриваемую задачу с не полностью определенными характеристиками исследуемого объекта. Характерной особенностью данной задачи является то, что
для каждого Р представляет собой не просто конечное температурное
Т (р,р)
распределение, а совокупность таких распределений температуры по пространственной координате / для всех возможных реализаций неопределенных факторов Т0 и а.
Вышесказанное означает, что температурные максимумы и минимумы могут достигаться не только в конкретно заданной точке /0 , а в разных точках заданного сечения нагреваемой заготовки для разных реализаций неопределенных факторов. Поскольку значения Т0° и а для каждой из точек /ф заранее не определены, возникает проблема их поиска, а также установления формы температурных распределений Т(р,Р) на множестве значений Р с целью записи расчетных систем уравнений альтернансного метода. Указанная проблема достаточно просто решается для типичных ситуаций, когда форма кривой в задаче с полным объемом информации Т (Я — г* ,/,Р) известна [4].
Пусть для данной задачи максимально достижимая точность нагрева обозначается через величину ё^, ф = 1,N . Тогда каждая из величин ёФ будет определяться с помощью следующего выражения:
ёт„ = т-'„
тах | Т (р, Р)
- Т
(16)
В данном случае означает максимально достижимую точность нагрева заготовки в классе задач с ] оптимизируемыми параметрами на всем множестве допустимых значений начальной температуры и коэффициента теплообмена.
Таким образом, методика решения рассматриваемой задачи с интервальной неопределенностью характеристик объекта может быть описана следующим образом.
На первом этапе решается задача в предположении, что максимально достижимая точность нагрева оказывается равной ё^ . Согласно [4] данное значение достигается при N= 1, то есть при наличии одного неизвестного параметра вектора Р .
Форма кривой Т(Я — г ,/,Р° ) для всех возможных комбинаций значений Т0
и а для наиболее типового случая, когда тепловые потери с поверхности нагреваемого изделия довольно малы, показана на рис. 3.
Согласно (16) температурное распределение, соответствующее значению
точности нагрева ё^1,, имеет две точки Рх = (Т00,а°)и Р0 =(То>а°),/0 предельно допустимого отклонения температуры от заданной в конце оптимального процесса. Поскольку форма температурного распределения Т (Я — г* ,/,Р° ) из-
вестна и имеет вид, показанный на рис. 3, то это означает, что известны и координаты точек с максимальным отклонением ¡^ = 0 и ¡° = ¡э2 для любой фиксированной пары значений Т0 и а . Таким образом, остается только определить неизвестные значения Т0°,а° и Т0о,а0, а следовательно, и кривые, на которых достигаются температурные экстремумы в точках ¡^ = 0 и ¡° = ¡э2 (рис. 4).
Рис. 3. Форма кривой Т (К — г Хр0 ) результирующего температурного распределения при точности нагрева г^
Рис. 4. Температурное распределение по пространственной координате ¡ в конце оптимального по точности нагрева при наличии неполной информации о характеристиках объекта при г = г"*1
Из показанного графика видно, что значения ,ах определяют верхнюю, а значения T00,а0 - нижнюю температурную кривую соответственно. Именно между этими условными границами и будут лежать остальные кривые температурных распределений, которые соответствуют всем другим, отличным от T0 0, а0
и T0 ,а0, значениям, удовлетворяющим ограничениям (13).
Очевидно, что при одинаковой длительности оптимального процесса и значениях неизвестных параметров вектора P для всей совокупности объектов максимум температуры в точке l2 = /э2 достигается при максимально возможной начальной температуре заготовки и минимально возможном коэффициенте тепловых потерь. Минимуму температуры в точке /0 = 0, напротив, соответствуют уже минимальная начальная температура и максимальный коэффициент тепловых потерь.
Таким образом, для точности нагрева S = S^ получаем следующую систему уравнений:
T (R " Л 0,TOmin,amax^JPil )" T
T (r - г,/0 ,T0max,amin^Jp10)-T* =s^; (17)
dT (R - r , l2,70
max * amin,p0 ) dl " '
Данная система оказывается замкнутой и может быть решена относительно всех неизвестных параметров процесса, к которым относятся p°, l2 и .
В случае, когда вектор конструктивных параметров P имеет более одной компоненты, необходимо перейти к решению задачи для нового значения точности нагрева s < ¿^П. Для этого решают последовательность задач для ряда значений s=s', уменьшающихся с малым шагом от ¿Ф. При этом для каждого значения точности нагрева необходимо задаваться формой температурного распределения, соответствующей данному значению. Подробно алгоритм решения данной задачи описан в [4].
Вновь принимая во внимание тот факт, что максимальная температура в конце оптимального процесса достигается при минимальном коэффициенте потерь и максимальной начальной температуре, а минимумы - при максимальном коэффициенте потерь и минимальной начальной температуре, для сформулированной задачи при N = 4 можно составить расчетную систему уравнений альтернансного метода, подобно выражению (15).
Система (18), состоящая из восьми уравнений, вновь оказывается замкнутой относительно всех неизвестных параметров процесса, к которым относятся: че-
D0
тыре компонента вектора оптимизируемых конструктивных параметров P , три координаты точек с максимальным отклонением температуры от заданной
l2,l30,l° и значение точности нагрева s^, и, следовательно, может быть решена
относительно рассматриваемых неизвестных.
T (R - r\ 10 ,T0min,ama^ P0)
a pu \-t' =-ё(4) i = 1 3 5-"mas.P I T fcmin'1 1'3'5'
T (R -/0 .Tomax."min, P° )
5T í R -/./«Tomin.ama, P
5/
5TíR-Л40.Tomas.am. P
5/
-Ts = f(4) j = 2 4-
= 0.0 = /0 < /0 <... < /50 = L;
= 0. k = 2.4.
Численное решение задачи и анализ полученных результатов
Для численного решения сформулированной задачи оптимального проектирования в программном пакете FLUX была разработана нелинейная двумерная модель индукционной нагревательной системы для закалки стальной цилиндрической заготовки [10-12]. Исходные данные для моделирования представлены в табл. 1.
Таблица 1
Исходные данные для построения модели
Параметр Обозначение Значение Ед. изм.
Частота питающего тока f 100000 Гц
Требуемая температура заготовки Т' 900 °C
Температура окружающей среды Т окр 20 °C
Начальная температура заготовки Т Т0 20 °C
Количество витков индуктора N 2
Радиус заготовки RB 14 мм
Длина заготовки LB 32,9 мм
Время нагрева t 8 c
Толщина закаливаемого слоя r" 4 мм
Материал заготовки Сталь (С40) [13]
Общий вид FLUX модели рассматриваемой индукционной нагревательной системы со сгенерированной конечно-элементной сеткой и результирующее температурное поле по длине заготовки показаны на рис. 5.
б
а
Рис. 5. Численная FLUX модель нагревательной системы а - конечно-элементная сетка; б - температурное поле по длине заготовки
в конце нагрева
Рассмотрим сначала решение задачи оптимального проектирования с полной информацией об объекте, когда начальная температура заготовки и коэффициент лучистого теплообмена принимают заранее известные фиксированные значения. В частном случае эти коэффициенты соответственно равны:
Вт
Г0 = 20°С;аи = 0.85-^^.
м • С
На неизвестные оптимизируемые параметры конструкции индукционной нагревательной установки накладываются следующие ограничения:
5 < ^ < 15,2 < p2 < 10,3 < pз < 15,500 < i < 1500. (19)
В табл. 2 приведены заданные начальные и оптимальные значения всех неизвестных параметров, найденные в результате решения системы уравнений (15).
Таблица 2
Результаты решения задачи оптимального проектирования без учета интервальной неопределенности
Параметр p p2 I е® mm
Начальное значение 12 мм 3 мм 4,45 мм 1220 А -
Оптимальное значение 11,9 мм 3 мм 4,4 мм 1382 А 11,5 °С
Температурное распределение вдоль границы закаливаемого слоя в конце стадии нагрева в индукционной нагревательной установке оптимальной конструкции показано на рис. 6.
х°с
915.00
885,00
0 5 10 15 20 25 30 35/= мм
Рис. 6. Конечное температурное распределение вдоль границы закаливаемого слоя
толщиной г
Отметим, что в [14] показано, что полученные при решении альтернансным методом результаты совпадают с результатами, найденными при использовании генетических алгоритмов оптимизации.
Рассмотрим теперь случай, когда начальная температура заготовки и коэффициент лучистого теплообмена задаются с точностью до следующих интервалов возможного изменения их значений:
0.5<а <1.3,15<Т <25. (20)
Ограничения на оптимизируемые параметры дизайна индуктора, а также их начальные приближения совпадают со значениями, использованными при решении задачи с полным объемом информации об объекте.
Значения оптимизируемых параметров, найденные в результате решения системы уравнений (18), представлены в сравнении с решением задачи оптимального проектирования без учета интервальной неопределенности в табл. 3.
Таблица 3
Сравнение оптимальных значений конструктивных параметров индуктора в задаче с полной информацией об объекте и в условиях интервальной неопределенности
Параметр Р1 Р2 Рз I £(4) тт
Задача с полной информацией 11,9 мм 3 мм 4,4 мм 1382 А 11,5 °С
Задача с интервальной неопределенностью 12,2 мм 3,6 мм 4,1 мм 1417 А 16,9 °С
Как уже было отмечено, температурные максимумы и минимумы достигаются в точках /0 ,ф = 1,5, расположенных на границе закаливаемого слоя в разных сечениях нагреваемой заготовки, координаты которых зависят от конкретных
реализаций неопределенных факторов. При этом максимальные отклонения температуры от заданной достигаются при максимальном значении начальной температуры заготовки и минимальном значении коэффициента теплообмена, а минимальные температурные отклонения наблюдаются при минимальном значении температуры и максимальном значении коэффициента тепловых потерь.
Таким образом, в результате решения рассматриваемой задачи получены температурные распределения для сечений с точками максимума и минимума температуры (рис. 7).
На рис. 7 кривая 1 соответствует максимальной начальной температуре
Вт
Т0тах = 25°С и минимальному значению коэффициента потерь аи = 0.5——^ ,
м ■ С
кривая 2 - минимальной начальной температуре Т0тах = 15°С и максимальному
значению коэффициента потерь аи = 1.3—.
м ■ С
Сравнение полученных результатов с результатами, полученными при решении задачи оптимального проектирования с полным объемом информации (см. табл. 3), показало, что введение интервальной неопределенности привело к увеличению максимального отклонения температуры от требуемого на 5 °С. Подобное снижение точности нагрева объясняется тем, что при решении задачи с не полностью определенными параметрами исследуемого объекта получены результаты, соответствующие наихудшей комбинации неопределенных факторов, к которым относятся начальная температура обрабатываемой заготовки и коэффициент тепловых потерь излучением.
Т,°С 920
SSO
О 5 10 15 20 25 30 35 Дым
Рис. 7. Температурные распределения в задаче оптимального проектирования в условиях наличия интервальной неопределенности:
1 - в сечении с максимальными температурами; 2 - в сечении с минимальными температурами
Заключение
В работе рассмотрена задача оптимального проектирования нагревательной установки для реализации стадии нагрева процесса поверхностной индукционной закалки стальных заготовок в условиях интервальной неопределенности начальной температуры детали и коэффициента лучистого теплообмена. Предложена методика решения сформулированной задачи с помощью альтернансного метода параметрической оптимизации. Приведены результаты численных расчетов на базе двумерной численной модели, построенной в программном пакете Altair FLUX.
Представленные расчеты показали, что введение интервальной неопределенности приводит к увеличению максимального отклонения температуры вдоль границы закаливаемого слоя заготовки от требуемого значения на 5 °С, что составляет менее 1 % от заданной температуры 900 °С. Подобное снижение точности нагрева является допустимым для технологии индукционной закалки.
БИБЛИОГРАФИЧЕСКИЙ СПИСОК
1. Rudnev V.I. et al. Handbook of Induction Heating. New York, Basel: Marcel Dekker Inc., 2003. 797 p.
2. БлантерМ.Е. Mеталловедение и термическая обработка. - M.: Mашгиз, 1963. - 416 с.
3. Плешивцева Ю.Э., Попов А.В., Попова М.А., Деревянов М.Ю. Оптимальное проектирование индуктора для поверхностной закалки цилиндрических заготовок на основе численной двумерной модели // Вестн. Астрахан. гос. техн. ун-та. Серия «Управление, вычисл. техн. ин-форм.». - 20i9. - № i. - С. 40-50.
4. Рапопорт Э.Я., Плешивцева Ю.Э. Оптимальное управление температурными режимами индукционного нагрева. - M.: Наука, 2012. - 309 c.
5. Рапопорт Э.Я. Оптимизация процессов индукционного нагрева металла. - M.: Mеталлургия, i993. - 279 c.
6. Плешивцева Ю.Э., Попов А.В., Дьяконов А.И. Двумерная задача оптимального но типовым критериям качества управления процессом сквозного индукционного нагрева // Вестник Самарского государственного технического университета. Сер. Технические науки. - 20i4. -№ 2 (42). - С. 148-í63.
7. Popov A.V. Optimization of Heating Stage for Induction Hardening of Cylindrical Billets // 20i9 XXI International Conference Complex Systems: Control and Modeling Problems (CSCMP). - Samara, Russia, 20i9. - pp. 237-24i.
S. Рапопорт Э.Я. Оптимальное управление системами с распределенными параметрами. - M.: Высшая школа, 2009. - 67S c.
9. Рапопорт Э.Я. Альтернансный метод в прикладных задачах оптимизации. - M.: Наука, 2000. - 33б c.
10. Flux [Электронный ресурс]: сайт программного продукта FLUX. URL: www.altair.com/flux/ (дата обращения 01.09.2020).
11. Шарапова О.Ю. Численное моделирование процесса периодического индукционного нагрева на базе конечно-элементного программного пакета FLUX // Вестник Самарского государственного технического университета. Сер. Технические науки. - 20ii. - № 7 (2S). - C. iS0-iS5.
12. Pleshivtseva Yu., Rogachev G., Popov A. MATLAB-FLUX Coupling for numerical modeling in education // SHS Web of Conferences 29,02033 (2016).
13. Mарочник сталей и сплавов. 2-е изд., дон. и иснр. / А.С. Зубченко, M.M. Колосков, Ю.В. Каширский и др. Под общей ред. А.С. Зубченко. - M.: Mашиностроение, 2003. - 784 с.
14. Pleshivtseva Yu., Baldan M., Popov A., Nikanorov A., Rapoport E., Nacke B. Effective methods for optimal design of induction coils on example of surface hardening // COMPEL - The international journal for computation and mathematics in electrical and electronic engineering, Vol. 39. No. 1, pp. 90-99. https://doi.org/10.1108/COMPEL-05-2019-0216.
Статья поступила в редакцию 1 июля 2020 г.
OPTIMAL INDUCTOR DESIGN FOR SURFACE HARDENING UNDER CONDITIONS OF INTERVAL UNCERTAITY OF PROCESS PARAMETERS
A.V. Popov
Samara State Technical University
244, Molodogvardeyskaya str., Samara, Russian Federation
Abstract. The paper is devoted to the optimal inductor design for surface hardening of steel cylindrical billets. The heating stage of surface induction hardening is considered as an object with distributed parameters, which unknowns are design characteristics of the induction installation. In real industrial conditions the main technological parameters are often defined by the intervals of their possible values. That is why, in the paper the optimal design problem under the conditions of interval uncertainty of initial billet's temperature and thermal exchange coefficient is formulated. Solution of the formulated problem is carried out by alternance method of parametric optimization based on numerical model developed in Altair FLUX software.
Keywords: surface hardening, induction heating, optimal design, alternance method, interval uncertainty, numerical model, electromagnetic and thermal fields, Altair FLUX.
REFERENCES
1. Rudnev V.I. et al. Handbook of Induction Heating. - New York, Basel: Marcel Dekker Inc., 2003. -797 p.
2. Blanter M.E. Metallovedenie i ternicheskaya obrabotka [Metal science and heat treatment]. Moscow: Mashgiz, 1963. 416 s.
3. Pleshivtseva Yu.E., Popov A.V., Popova M.A., Derevyanov M.Yu. Optimalnoye proektirovaniye in-duktora dlya poverkhnostnoi zakalki tsilindricheskikh zagotovok na osnove chislennoi dvumernoi modeli [Optimal inductor design for surface hardening of cylindrical billets] // Bulletin of Astrakhan State Technical Unversity. Series "Control, computer engineering, informatics". - Astrakhan, 2019. - № 1. - pp. 40-50.
4. Rapoport E.Ya., Pleshivtseva Yu.E. Optimalnoe upravleniye temperaturnimi rezhimami induktsion-nogo nagreva [Optimal control of induction heating processes]. Moscow: Nauka, 2012. 309 p.
5. Rapoport E.Ya. Optimizatsiya protsessov induktsionnogo nagreva metalla [Optimization of induction heating processes]. Moscow: Metallurgiya, 1993. 279 p.
6. Pleshivtseva Yu.E., Popov A.V., Diakonov A.I. Dvumernaya zadacha optimalnogo po tipovim kriteriyam kachestva upravleniya protsessom skvoznogo induktsionnogo nagreva [Two-dimensional problem of optimal with respect to typical quality criteria control of through induction heating processes] // Vestnik of Samara State Technical University. Technical Sciences Series. Samara, 2014. № 2(42). pp. 148-163.
7. Popov A.V. Optimization of Heating Stage for Induction Hardening of Cylindrical Billets // 2019 XXI International Conference Complex Systems: Control and Modeling Problems (CSCMP). - Samara, Russia, 2019. pp. 237-241.
Anton V. Popov, Senior Lecturer.
8. Rapoport E.Ya. Optimalnoe upravleniye sistemami s raspredelennimi parameterami [Optimal control of systems with distributed parameters]. - Moscow: Visshaya shkola, 2009. 678 p.
9. Rapoport E.Ya. Alternansniy metod v prikladnikh zadachakh optimizatsii [Alternance method in applied optimization problems]. Moscow: Nauka, 2000. 336 p.
10. Flux [electronic source]: www.altair.com/flux/ (accessed 01.09.2020).
11. Sharapova O.Yu. Chislennoye modelirovaniye protsessa periodicheskogo induktsionnogo nagreva na baze konechno-elementnogo programmnogo paketa FLUX [Numerical simulation of static induction heating process based on finite-element software FLUX] // Vestnik of Samara State Technical University. Technical Sciences Series. Samara, 2011. № 7 (28). pp. 180-185.
12. Pleshivtseva Yu., Rogachev G., Popov A. MATLAB-FLUX Coupling for numerical modeling in education // SHS Web of Conferences 29,02033 (2016).
13. Marochnik staley i splavov 2-e izdaniye [Grade guide of steels and alloys. Second edition] / A.S. Zubchenko, M.M. Koloskov, Yu. V. Kashirskiy et al. - Moscow: Mashinostroyeniye, 2003. 784 p.
14. Pleshivtseva Yu., Baldan M., Popov A., Nikanorov A., Rapoport E., Nacke B. Effective methods for optimal design of induction coils on example of surface hardening // COMPEL - The international journal for computation and mathematics in electrical and electronic engineering, Vol. 39 No. 1, pp. 90-99.
