CC BY
21ВЕСТН. САМАР. ГОС. ТЕХН. УН-ТА. СЕР. ТЕХНИЧЕСКИЕ НАУКИ. 2021. Т. 29. № 3_
doi: 10.14498/tech.2021.3.2
УДК 621.785, 669.14, 519.6
ОПТИМИЗАЦИЯ КОНСТРУКТИВНЫХ И РЕЖИМНЫХ ПАРАМЕТРОВ ИНДУКЦИОННОЙ НАГРЕВАТЕЛЬНОЙ УСТАНОВКИ ДЛЯ ЗАКАЛКИ*
А.В. Павлушин
Самарский государственный технический университет Россия, 443100, г. Самара, ул. Молодогвардейская, 244
E-mail: [email protected]
Аннотация. В работе рассматривается проблема оптимизации конструктивных и режимных параметров индукционной нагревательной системы для поверхностной закалки стального цилиндрического вала сложной геометрии, содержащего переход с малого диаметра на большой. Задача оптимального проектирования индуктора формулируется применительно к нелинейной двумерной численной модели взаимосвязанных электромагнитных и температурных полей на стадии индукционного нагрева, разработанной в ППП ANSYS Mechanical APDL. Описаны этапы решения задачи оптимального проектирования на базе альтернансного метода параметрической оптимизации систем с распределенными параметрами. Алгоритм автоматической оптимизационной процедуры реализуется в ППП MATLAB, предоставляющем возможность интеграции численной ANSYS модели. Представлен сравнительный анализ результатов моделирования индукционной нагревательной системы с типовой квадратной формой витков и системы с витками индуктора усложненной геометрической формы, оптимизация которой позволила получить существенно более равномерное температурное распределение по границе закаливаемого слоя.
Ключевые слова: численное моделирование, альтернансный метод, оптимальное проектирование, индукционный нагрев, поверхностная закалка, электромагнитное поле, температурное распределение, ANSYS, MATLAB.
Введение
Шестерни, подшипники, валы и другие металлические детали различного назначения во время работы испытывают постоянные значительные осевые или радиальные нагрузки. Для повышения надежности и срока службы такие детали подвергаются термической обработке. Поверхностная закалка представляет собой один из наиболее распространенных способов термообработки, при котором металлическая деталь нагревается до определенной температуры, выдерживается заданное время при этой температуре и охлаждается со скоростью, необходимой для образования требуемых микроструктурных свойств. Поверхностное упрочнение позволяет получить сочетание требуемой твердости поверхностного слоя и необходимой прочности, ударной вязкости и пластичности сердцевины детали. Стадия нагрева технологии поверхностной закалки часто реализуется индукционным способом при использовании токов высокой частоты.
* Работа выполнена при поддержке гранта РФФИ 19-08-00232.
Павлушин Алексей Владимирович, инженер НИС кафедры «Автоматика и управление в технических системах», аспирант.
Основной целью исследования является оптимизация конструктивных и режимных параметров индукционной нагревательной установки для обеспечения максимально достижимой равномерности температурного распределения в поверхностном закаливаемом слое стального цилиндрического вала, содержащего переход с малого диаметра на большой [1].
Обеспечение равномерного заданного температурного распределения в поверхностном слое заготовок сложной геометрической формы является нетривиальной задачей, решение которой с требуемой точностью позволяет уменьшить растягивающие эквивалентные нагрузки и снизить риск появления микротрещин в процессе эксплуатации, что существенно улучшает качество изготавливаемых деталей в соответствии с постоянно возрастающими технологическими требованиями, предъявляемыми к ним.
Исследование влияния конструктивных и режимных параметров
на равномерность температурного поля
Для исследования влияния конструктивных и режимных параметров на равномерность температурного поля в поверхностном слое заготовки, имеющей переход с одного диаметра на другой, в процессе ее нагрева под закалку были разработаны двумерные нелинейные численные модели в ППП ANSYS Mechanical APDL [1, 2].
Процесс моделирования включает несколько этапов. На первом этапе задаются необходимые физические свойства в ППП ANSYS Mechanical APDL [3, 4]. Далее по исходным данным (табл. 1) строится осесимметричная двумерная геометрия модели, выполняется присвоение заранее внесенных свойств материалов соответствующим областям, накладывается конечно-элементная сетка по специальному алгоритму, состоящему из набора инструментов и функций, предлагаемых программным пакетом. Далее решаются взаимосвязанные электромагнитная и тепловая задачи, которые в общем случае описываются известной системой взаимосвязанных уравнений Максвелла и Фурье [1, 5-9]. Общий вид исходной геометрии системы «индуктор - заготовка» с типовой квадратной формой витков индуктора представлен на рис. 1.
На рис. 1а показаны следующие конструктивные параметры, влияние которых на температурное поле исследуется в работе: pi (размер квадратного сечения витков), рз (расстояние, характеризующее положение витков относительно заготовки, - воздушный зазор), р4 (расстояние между витками 1 и 2), рз (расстояние между витками 2 и 3). Кроме того, исследовалось влияние режимных параметров процесса нагрева: р2 (силы тока на витках 1 и 3) и рв (силы тока на витке 2). Исходные данные для моделирования представлены в табл. 1.
Таблица 1
Исходные данные для моделирования
№ Наименование параметра Значение
1 Начальная температура заготовки Та, °С 20
2 Длительность стадии индукционного нагрева т0, сек 10
3 Требуемая температура нагрева вдоль границы закаливаемого слоя Т*, °С 900
4 Высота заготовки Ь, мм 51
№ Наименование параметра Значение
5 Длина границы закаливаемого слоя стальной заготовки, мм 60
6 Радиус заготовки Я1, мм 11
7 Радиус заготовки Я2, мм 19
8 Глубина поверхностного закаливаемого слоя С, мм 2
9 Типовая углеродистая сталь (марка) С40
б
Рис. 1. Двумерная осесимметричная модель системы «индукционный нагреватель стальная заготовка» с типовой квадратной формой витков:
а - общий вид с обозначением конструктивных параметров; б - наложенная конечно-элементная сетка в ППП ANSYS
Рис. 2. Геометрия угловой зоны стальной заготовки с обозначениями точек в ППП ANSYS
а
Анализ результатов численного моделирования процесса индукционного нагрева заготовки, имеющей переход с одного диаметра на другой, демонстрирует сильное влияние угловой зоны на распределение электромагнитного поля, связанное с существенной вариацией плотности магнитного потока в указанной области, вследствие чего температурное распределение в этой области становится крайне неравномерным. При этом в области внутреннего угла заготовки, образованного границей закаливаемого слоя вдоль точек 13, 19 и 20 (рис. 2), наблюдается недопустимый недогрев, а внешний угол, образованный точками 23, 31 и 32 (рис. 2), недопустимо перегревается. При начальных значениях параметров р — р6 (табл. 2) максимальное отклонение температуры на границе закаливаемого слоя от требуемого значения Т* составило более 200 °С, что совершенно недопустимо.
Варьирование исследуемых параметров рх — р6 в широком диапазоне позволило снизить температурное отклонение только до ±156°С, что также совершенно не удовлетворяет технологическим требованиям. Соответствующие значения исследуемых параметров также представлены в табл. 2 (строка 2).
Таблица 2
Начальные и конечные значения исследуемых параметров
Значения р1, мм Р2, А р3, мм р4, мм р5, мм Р6, А
Исходные 13,00 1054,00 2,37 2,64 6,91 1079,00
Конечные 8,96 913,47 1,92 2,96 9,85 1221,13
Рис. 3. Двумерная осесимметричная модель системы «индукционный нагреватель -стальная заготовка» с усложненной геометрией витка 2: а - общий вид с обозначением конструктивных параметров; б - наложенная конечно-элементная сетка в ППП А№У8
Результаты моделирования привели к выводам, что квадратная форма геометрии витков, представленная на рис. 1, не дает возможности обеспечить требуемую равномерность нагрева закаливаемого поверхностного слоя стальной заготовки сложной формы. Поэтому предлагается перейти к более сложной геометрической форме витка 2 индуктора с выступающей внешней частью для уменьшения резкого отклонения температуры в угловой зоне заготовки. Общий вид системы «индуктор - заготовка» с витком усложненной формы представлен на рис. 3а. На рис. 3б представлена геометрия данной системы с наложенной конечно-элементной сеткой.
Для оптимизации конструктивных и режимных параметров ИНУ, представленной на рис. 3, далее формулируется задача оптимизации, которая решается на базе альтернансного метода оптимизации систем с распределенными параметрами
[9].
Постановка задачи оптимального проектирования
Для формулировки задачи оптимизации стадии нагрева индукционная нагревательная система (рис. 3) рассматривается в качестве объекта управления с распределенными параметрами, который описывается представленной выше нелинейной двумерной численной моделью взаимосвязанных электромагнитных и температурных полей, разработанной в ППП ANSYS Mechanical APDL.
Оптимизируемый вектор P параметров включает следующие конструктивные характеристики (рис. 3): pi - расстояние между витком 1 и заготовкой; р2 -расстояние между витком 1 и координатой Y=0; рз - расстояние между витком 3 и
заготовкой; р4 - высота витка 1. Кроме того, вектор P включает следующие режимные параметры процесса нагрева: рз - силу тока витка 2; рв - силу тока витка 1. Значение силы тока для витка 3 было найдено ранее и определено в качестве константы /3 = 1020,24^ .
Требование к результирующему температурному распределению в поверхностном слое заготовки в конце стадии предварительного нагрева металла под закалку целесообразно формулировать в виде задания максимально допустимой величины s абсолютного отклонения результирующего температурного распределения от заданного значения T*. В соответствии с этим в качестве критерия оптимизации предлагается рассматривать критерий, минимизирующий температурные отклонения в пределах границы закаливаемого слоя, показанного пунктирной линией на рис. 1а и 3а и находящегося на глубине С=2 мм от поверхности, который может быть записан в виде [10, 11]:
T|X,т0,1,,P \-T
J = max
/ e[0,L ]
^ min, (1)
p
¡], Р )
где Т|Х,т0,/у,Р ) - температура в конце стадии нагрева на границе закаливаемого слоя с координатой X, располагающейся на глубине С от поверхности заготовки;
т0 - время нагрева;
Т* - требуемая температура на границе поверхностного слоя [10, 11].
На выбранные оптимизируемые конструктивные и режимные параметры Р накладываются ограничения:
РЫп ^ Р < Р^ Р2тп ^ Р2 ^ Р™, Рзтп ^ Рз ^ Рз™^ Ртп ^ Л ^ Р(2)
Лш, ^ 1 ^ х, ^2тт ^ ^ 12тах *
Задача оптимального проектирования индукционной нагревательной установки заключается в поиске таких значений параметров
поит / „опт „опт „опт „опт т опт г опт \
Р = (р1 ,,р3 ,р4 ,11 ,12 ), на которые наложены ограничения (2), обусловленные характеристиками нагревательной системы, которые за установленное время т0 обеспечат перевод объекта, описываемого нелинейной численной двумерной моделью, из исходного состояния в требуемое конечное состо-
яние, соответствующее минимальному значению критерия оптимальности (1) [10, 11].
Решение задачи оптимизации на основе альтернансного метода
Положительный опыт использования алгоритма на основе альтернансного метода для оптимизации температурного поля модели индукционного нагрева стальной заготовки цилиндрической формы позволяет сделать вывод о возможности его применения для работы с заготовкой сложной геометрической формы [1, 9, 10, 11, 12, 13]. Используемый метод в работе опирается на установленные свойства пространственных распределений температуры по объему заготовки в конце оптимальных процессов индукционного нагрева металла, подобные известным в математике свойствам наилучших приближений заданных функций к нулю. Теория альтернансного метода позволяет произвести процедуру точной редукции исходной задачи оптимизации к решению трансцендентных систем уравнений, замкнутых относительно всех искомых параметров [10, 11].
Согласно альтернансному методу, основное свойство распределения
Т(х,т0,/у,Р | состоит в том, что число К точек /°пт в сечениях обрабатываемой
заготовки вдоль границы закаливаемого слоя с координатой X, в которых достигаются предельно достижимые абсолютные отклонения конечной температуры от требуемой Т*, равные е, всегда оказываются не меньше числа N искомых оптимизируемых параметров процесса, причем эти отклонения Т(х,т0,/у,Р )-Т* являются знакочередующимися в точках /°пг,] = 1,К;^=±1, т. е. обладают альтер-нансным свойством [9].
Рис. 4. Эталонная форма конечного температурного распределения в классе задач
(6)
с шестью оптимизируемыми параметрами по критерию точности нагрева
Для оптимального проектирования конструкции и нахождения оптимальных режимных значений параметров нагревательной установки скомпонован вектор из
шести оптимизируемых параметров, и, следовательно, прогнозируется формирование семи точек с максимальным отклонением температуры от требуемой. С учетом вышеописанного по данному методу за эталонную кривую принята форма конечного температурного распределения, которая представлена на рис. 4 [10, 11].
На основании альтернансного метода для сформулированной задачи составлена система из 14 уравнений, замкнутая относительно всех неизвестных:
Т (х,т°, /;пг, Ропт)-Т* = (-¡У уе^ = 1,7; у = ±1;
ропт _I „опт л опт л опт ^опт т опт т опт \
Р = (Р1 , Р2 , Р3 , Р4 ,11 ,12 )
0 = /1опт < /2опт < -3™ < -40Ш < -5ош < < /7опт = Ь (3)
дТ (- опт ропт)—т *
—-)-= 0; т = 1,7,
д1
;опт ;опт ;опт ;опт ;опт ;опт ;опт (6)
где -1 ,-2 , -з , -4 , -5 ,-6 , -7 - координаты точек экстремума, е^ - значение предельно достижимой точности нагрева в классе задач с шестью неизвестными оптимизируемыми параметрами.
Система уравнений, составленная на основе альтернансного метода, решается в программном пакете МЛТЬЛБ. Нахождение вектора оптимальных значений искомых параметров обеспечивается стандартной функцией МЛТЬЛБ - LSQnonlin, которая позволяет решить задачу нахождения точки минимума некоторой целевой нелинейной функции Дх) методом наименьших квадратов или подгонки кривой наименьших квадратов. При этом целевая функция формируется следующим образом:
тт||/(х)||2 = т/)2 + /2(х)2 +... + /(х)2) , (4)
где в общем случае А(х) - вектор-функция; х - вектор-столбец искомых переменных [14, 15].
Численное решение задачи оптимизации
Исходные данные для численного решения задачи оптимизации представлены в табл. 3.
Таблица 3
Исходные данные для решения задачи оптимизации
Наименование параметра Значение
Виток № 1 Виток № 2 Виток № 3
1 Сила тока, А 823,56 923,23 1020,24
2 Ширина витка, мм 12,72 -(табл. № 4) 10,00
3 Высота витка, мм 11,42 -(табл. № 4) 15,61
4 Расстояние между витками 1 и 2, мм 1,83
5 Расстояние между витками 2 и 3, мм 6,51
6 Частота тока, кГц 200
7 Расстояние между витком 1 и заготовкой, мм 1,95
8 Расстояние между витком 3 и заготовкой, мм 2,27
Значения конструктивных параметров витка 2 сложной геометрической формы представлены в табл. 4.
Таблица 4
Геометрические характеристики витка сложной формы
Конструктивные размеры витка 2
№ Расстояние от заготовки до точки по оси X, мм Значение № Расстояние от заготовки до точки по оси Y, мм Значение
1 до точки х12 2,19 7 до точки у12 15,67
2 до точки x34 2,19 8 до точки у34 6,64
3 до точки х11 14,76 9 до точки у11 15,67
4 до точки х10 9,17 10 до точки у10 6,64
5 до точки х35 0,37 11 до точки у35 0,44
6 до точки х36 3,16 12 до точки у36 0,44
Рис. 5. Геометрия витка № 2 с обозначениями точек в ППП ANSYS
Следующим шагом производится запуск расчета оптимизационной процедуры, которая работает следующим образом (рис. 6). В ПП МАТЛАБ алгоритм создает файлы и прописывает в них значения исходных параметров. Далее ANSYS считывает эти файлы, строит и рассчитывает модель. Формируется температурное распределение. Затем в МАТЛАБ передается этот массив данных из ППП ANSYS, производится программный анализ и численное дифференцирование температурного распределения. Решается система из 14 уравнений по аль-тернансному методу и производится анализ результатов оптимизируещей функцией LSQnonlil.
Геометрия витка № 2 с обозначениями в ППП ANSYS представлена на рис. 5.
В процессе изменения значений координат витков может произойти наложение точек, определяющих геометрию, друг на друга, вследствие чего ANSYS выдаст сообщение об ошибке построения геометрии. По этой причине, а также для
исключения выхода значений силы тока за пределы рабочего диапазона источников питания и корректного определения прогнозируемых координат точек экстремума на оптимизируемые параметры были наложены ограничения вида:
0.01 < p < 10, 0.01 < р < 20, 0.01 < p < 20, 0.01 < р < 10; 10 < / < 1200, 10 < /2 < 1400; 0.01 < /0 < 10.189, 10.189 < ^ < 22.9, 22.9 < ^ < 29.62, 29.62 < /3 < 36.36, 36.36 < /4 < 44.44, 44.44 < /5 < 54.9,
54.9<4 <59.8; 10<£<500.
При нахождении целевой функцией локального минимума производится анализ сформировавшегося температурного распределения вдоль границы закаливаемого слоя, по которому определяется, удовлетворяет ли полученное решение системы уравнений условию существованию глобального оптимума. При этом для получения решения приходилось неоднократно изменять начальные значения компонент вектора оптимизируемых параметров.
Рис. 6. Структура численной оптимизационной процедуры
В табл. 5 представлены оптимальные значения параметров конструкции индуктора со сложной геометрией витка № 2.
Таблица 5
Конечные оптимальные значения параметров
Значения pi, мм p2, мм p3, мм p4, мм p5, А pe, А
Оптимальные 1,91 2,27 2,15 11,31 923,14 831,53
Координаты точек экстремумов температурного распределения при использовании сложной формы витка представлены в табл. 6.
Таблица 6
Координаты точек экстремумов оптимального температурного распределения по длине закаливаемого слоя
Значения /Оопг, мм /1опт, мм /Г, мм /опт „ „ , /3 , мм /40Ш, мм /опт „ „ , /5 , мм /Г, мм
Конечные 6,93 20,13 27,89 32,36 39,93 47,53 59,66
В результате решения системы уравнений по альтернансному методу максимальное температурное отклонение составило =51,87 °С.
На рис. 7 представлено температурное распределение по границе закаливаемого слоя стальной заготовки, имеющей переход с одного диаметра на другой, полученное в конце стадии нагрева в индукторе с квадратной формой витков (кривая 1), в сравнении с температурным распределением, достигаемым при оптимальных значениях конструктивных и режимных параметров (кривая 2) в индукторе с усложненной геометрической конструкцией витка № 2.
2 -
Рис. 7. Результирующее температурное распределение по границе закаливаемого слоя заготовки:
1 - для индуктора с квадратной формой витков; после оптимизации параметров индуктора со сложной с формой витка № 2
Сравнение двумерных температурных полей, полученных в конце стадии нагрева заготовки в продольном ее сечении, представлено на рис. 8.
Температурные поля (рис. 8 а, б) были получены в индукторе с квадратной формой витков с начальными и варьированными значениями конструктивных и режимных параметров.
б
а
в
Рис. 8. Результирующее распределение температуры в продольном сечении заготовки: а - для параметров индуктора, представленных в табл. 2 (строка 1); б - для параметров индуктора, представленных в табл. 2 (строка 2); в - для оптимальных параметров индуктора
Выводы
Анализ результатов моделирования выявил необходимость перехода к более сложной геометрической форме витка индуктора, находящегося в зоне наибольшей неравномерности электромагнитного поля в процессе поверхностного нагрева стальной заготовки, имеющей переход с одного диаметра на другой. Решение задачи оптимального проектирования позволило обеспечить значительное снижение максимального температурного отклонения от заданного значения на границе закаливаемого слоя ±156 °С до ±51,9 °С при отсутствии перегрева поверхности.
Таким образом, оптимизация конструктивных и режимных параметров индуктора на базе альтернансного метода приводит к значительному увеличению равномерности распределения температурных полей в пределах упрочняемого слоя заготовки сложной геометрической формы в конце стадии нагрева, что значительно улучшает качество процесса закалки в целом.
Дальнейшее улучшение равномерности нагрева поверхностного слоя возможно путем решения сформулированной задачи при расширении вектора оптимизируемых параметров, в том числе за счет дальнейшего усложнения геометрической формы витков индуктора, применения локальных магнитных концентраторов и увеличения числа режимных технологических параметров стадии нагрева.
В этом случае оптимизируемый вектор Р параметров может включать в себя следующие конструктивные характеристики:
- пространственно-варьируемый зазор между заготовкой и индуктором;
- оптимальная частота;
- электрический ток индуктора;
- время нагрева;
- возможность применения локальных концентраторов, способных локально
фокусировать электромагнитное поле.
Однако ожидаемые улучшения нетривиальны, нуждаются в подтверждении
и могут являться темой отдельного исследования.
БИБЛИОГРАФИЧЕСКИЙ СПИСОК
1. Павлушин А.В. Оптимизация процесса индукционной закалки металлических заготовок с использованием численной ANSYS-модели // XIV Международная молодежная научная конференция «Тинчуринские чтения». В 3 т. Т. 1. Электроэнергетика и электроника: матер. конф. (Казань, 23-26 апреля 2019 г.) / под общ. ред. ректора КГЭУ Э.Ю. Абдуллазянова. - Казань: Казан. гос. энерг. ун-т, 2019. - Ч. 2. - C. 21-30.
2. Павлушин А.В. Двумерная численная нелинейная модель температурных и термонапряженных состояний стальной заготовки в процессе индукционного нагрева под закалку // Н34 Наука. Технологии. Инновации // Сборник научных трудов в 9 ч. / под ред. А.В. Гадюкиной. - Новосибирск: Изд-во НГТУ, 2020. ISBN 978-5-7782-4288-3 Часть 5: C. 46-51.
3. Engineering Simulation & 3D Design Software | ANSYS. 2019 [Электронный ресурс]. - URL: https://www.ansys.com/ (дата обращения: 14.05.2021).
4. Engineering Simulation & 3D Design Software | ANSYS. 2019 [Электронный ресурс]. - URL: https://www.ansys.com/blog/how-to-make-matlab-apps-for-ansys/ (дата обращения: 14.05.2021).
5. Pleshivtseva Yu., Korshikov S., Wipprecht S., Makarov E., Baake E., Nacke B. Simulation of primary heating stage in resource efficient forging chain // Heat processing, 2015, № 1, Vol. 13, pp. 85-90.
6. Pleshivtseva Y.E., Yakubovich E.A. The influence of induction heating modes on thermal stresses within the billets // Advances in Automation, 2020 / Proceedings of the International Russian Automation Conference, RusAutoCon 2019, September 8-14, 2019, Sochi, Russia. Book Series: Lecture Notes in Electrical Engineering Editors: Prof. Andrey A. Radionov, Alexander S. Karandaev. Publisher: Springer International Publishing Print. ISBN: 978-3-030-39224-6. Electronic ISBN: 978-3030-39225-3. Part of: Springer Professional Wirtschaft+Technik", Springer Professional "Technik".
7. Rudnev V. et al. Handbook of induction heating // Manufacturing Engineering and Materials Processing. 2003. Vol. 61.
8. Источники питания. Математическое моделирование и оптимизация. Интенсивные курсы. Основы II. Учебное пособие / Бааке Э., Барглик Д., Долега Д., Лупи С. и др. - СПБ.: Изд-во СПбГЭТУ «ЛЭТИ», 2013. - 324 с.
9. Попов А.В. Оптимальное проектирование индукционного нагревателя для поверхностной закалки в условиях интервальной неопределенности параметров // Вестник Самарского государственного технического университета. Сер.: Технические науки. - 2020. - Т. 28. - № 3. -C. 139-154. doi: 10.14498/tech.2020.3.9
10. Рапопорт Э.Я., ПлешивцеваЮ.Э. Оптимальное управление температурными режимами индукционного нагрева. - М.: Наука, 2012.
11. Рапопорт Э.Я. Альтернансный метод в прикладных задачах оптимизации. - М.: Наука, 2000. -С. 336-336.
12. ПлешивцеваЮ.Э., Попов А.В., ПоповаМ.А., ДеревяновМ.Ю. Оптимальное проектирование индуктора для поверхностной закалки цилиндрических заготовок на основе численной двумерной модели // Вестник Астраханского государственного технического университета. Сер.: Управление, вычислительная техника и информатика. - 2019. - № 1. - С. 40-50. DOI: 10.24143/2072-9502-2019-1-40-50.
13. Pleshivtseva Y., 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, 2019, ISSN: 03321649, 39 (1), pp. 90-99.
14. MathWorks - Makers of MATLAB and Simulink - MATLAB & Simulink. 2021 [Электронный ресурс]. - URL: https://www.mathworks.com/ (дата обращения: 25.01.2021).
15. MATLAB. Exponenta. Optimization Toolbox 2.2. Руководство пользователя / авт. А.Г. Трифонов [Электронный ресурс]. - URL: http://matlab.exponenta.ru/optimiz/book_1/ ( дата обращения: 02.04.2021).
Статья поступила в редакцию 10 апреля 2021 г.
OPTIMIZATION DESIGN AND OPERATING PARAMETERS OF INDUCTION HEATING SYSTEM FOR HARDENING
A.V. Pavlushin
Samara State Technical University
244, Molodogvardeyskaya st., Samara, 443100, Russian Federation
Abstract. The paper deals with the problem of optimizing the design and operating parameters of an induction heating system for surface hardening of a steel stepped shaft. The problem of optimal design of an inductor is formulated based on a nonlinear two-dimensional numerical model of coupled electromagnetic and temperature fields, developed in the ANSYS Mechanical APDL software. Alternance method of parametric optimization of systems with distributed parameters is used to optimize induction hardening system. MATLAB software has been used for developing parametric optimization subroutine, which was incorporated into the numerical ANSYS model to simulate a process of induction heating. Commonly used a multi-turn solenoid-style coil fabricated from rectangular copper tubing has been used as a hardening inductor. Besides that, an application of profiled copper turns has been investigated. Optimization of induction hardening system described above allows one to substantially improve heating uniformity and enhance metallurgical characteristics of as-hardened stepped shaft. Localized temperature surplus at an upper diameter shoulder has been minimized. At the same time, sufficient austenitization in the fillet area near stepped region (diameter transition) has been obtained.
Keywords: numerical simulation, alternance method, optimal design, induction heating, surface hardening, stepped shaft, induction hardening, electromagnetic field, temperature distribution, ANSYS, MATLAB.
REFERENCES
1. Pavlushin A. V. Optimization of the process of induction hardening of metal blanks using a numerical ANSYS model // XIV International Youth Scientific Conference "Tinchurin Readings". In 3 volumes. V. Power engineering and electronics: mater. conf. (Kazan, April 23-26, 2019) / under total. ed. Rector of KSPEU E.Yu. Abdullazyanov. - Kazan: Kazan. state energ. un-t, 2019. - Part 2. - Pp. 21-30 (In Russian).
2. Pavlushin A.V. Two-dimensional numerical nonlinear model of temperature and thermally stressed states of a steel billet in the process of induction heating for hardening // N34 Nauka. Technologii. Innovations // Collection of scientific papers in 9 hours / ed. Gadyukinoy A.V. -Novosibirsk: Publishing house of NSTU, 2020. ISBN 978-5-7782-4288-3 Part 5: Pp. 46-51 (In Russian).
3. Engineering Simulation & 3D Design Software | ANSYS. 2019. URL: https://www.ansys.com/ (14.05.2021). '
4. Engineering Simulation & 3D Design Software | ANSYS. 2019. URL: https://www.an-sys.com/blog/how-to-make-matlab-apps-for-ansys/ (14.05.2021).
5. Pleshivtseva Yu., Korshikov S., Wipprecht S., Makarov E., Baake E., Nacke B. Simulation of primary heating stage in resource efficient forging chain // Heat processing. 2015. № 1. Vol. 13. Pp. 85-90.
6. Pleshivtseva Y.E., Yakubovich E.A. The influence of induction heating modes on thermal stresses within the billets // Advances in Automation, 2020 / Proceedings of the International Russian Automation Conference, RusAutoCon 2019, September 8-14, 2019, Sochi, Russia. Book Series: Lecture Notes in Electrical Engineering Editors: Prof. Andrey A. Radionov, Alexander S. Karan-daev Publisher: Springer International Publishing Print ISBN: 978-3-030-39224-6. Electronic ISBN: 978-3-030-39225-3. Part of: Springer Professional Wirtschaft+Technik", Springer Professional "Technik".
7. Valery R. et al. Handbook of induction heating // Manufacturing Engineering and Materials Processing. 2003. T. 61.
Aleksey V. Pavlushin, Postgraduate Student.
8. Power supplies. Mathematical modeling and optimization. Intensive Basics Courses II. Tutorial / Baake E., Barglik D., Dolega D., Lupi S. et al. // SPB.: Publishing house of ETU "LETI", 2013. 324 p. (In Russian).
9. Popov A. V. Optimal inductor design for surface hardening under conditions of interval uncertaity of process parameters // Vestnik of Samara State Technical University. Technical Sciences Series. 2020. Vol. 28. № 3. Pp. 139-154. doi: 10.14498/tech.2020.3.9 (In Russian).
10. RapoportE.Ya., Pleshivtseva Yu.E. Optimal control of temperature modes of induction heating. Moscow: Nauka, 2012. (In Russian).
11. Rapoport E.Ya. Alternative method in applied optimization problems. M.: Nauka, 2000. Pp. 336336 (In Russian).
12. Pleshivtseva Yu.E., PopovA.V., PopovaM.A., DerevyanovM.Yu. Optimal design of an inductor for surface hardening of cylindrical blanks on the basis of a numerical two-dimensional model. Bulletin of the Astrakhan State Technical University. Series: Management, Computer Engineering and Informatics. 2019. No. 1. Pp. 40-50. DOI: 10.24143 / 2072-9502-2019-1-40-50 (In Russian).
13. Pleshivtseva Y., 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, 2019. ISSN: 0332-1649, 39 (1). Pp. 90-99.
14. MathWorks - Makers of MATLAB and Simulink - MATLAB & Simulink. 2021. URL: https://www.mathworks.com/ (accessed: 25.01.2021).
15. MATLAB. Exponenta Vvedenie v sovremennye metodu optimizatsii system upravleniya [Introduction to modern methods of control systems optimization]. http://matlab.exponenta.ru/op-tim-robast/book1/index.php (accessed April 2, 2021).
