Проблемно-ориентированная модель технологического процесса поверхностной индукционной закалки

Обложка


Цитировать

Полный текст

Аннотация

Представлена нелинейная двумерная численная модель температурных полей и полей термических напряжений в процессе поверхностной индукционной закалки стальной заготовки сложной геометрической формы. Проблемно-ориентированная модель разработана в ППП ANSYS Mechanical APDL для интеграции в процедуры численного решения задач оптимального проектирования индуктора и управления стадиями нагрева и охлаждения. Описаны этапы моделирования: задание свойств материалов, построение геометрии и наложение КЭ-сетки, задание краевых условий и настройка алгоритма расчета модели. Для моделирования стадии нагрева решаются взаимосвязанные электромагнитная, тепловая и упругопластическая задачи, в результате чего определяются пространственно-временные распределения электромагнитных источников тепла, температурных полей и полей термических напряжений по объему заготовки. Модель стадии нагрева ориентирована на дальнейшее использование в численной процедуре параметрической оптимизации конструктивных и режимных параметров индукционной установки с целью получения максимально возможной равномерности температурного распределения, от которой существенно зависит качество микроструктуры закаливаемого слоя. Для моделирования стадии охлаждения решаются тепловая и упругопластическая задачи, в результате чего определяются пространственно-временные распределения полей температуры и термонапряжений по объему заготовки. В качестве охлаждающей (закалочной) среды используется водяной душ, применяемый после индукционного нагрева при поверхностной закалке для формирования мартенситной микроструктуры. Расчет стадии охлаждения выполняется при использовании эквивалентного коэффициента поверхностной теплоотдачи, зависящего от температуры, что значительно упрощает процесс моделирования. Модель стадии охлаждения может быть интегрирована в численную процедуру решения задачи оптимального управления спрейерным устройством, которое позволяет контролировать закалку путем регулировки интенсивности подачи воды.

Полный текст

Введение

Использование современных технологий поверхностной индукционной закалки позволяет упрочнять детали различной геометрической формы. Благодаря высокой скорости протекания процесса и концентрации внутренних источников тепла в поверхностном слое этот метод обладает высоким энергосберегающим потенциалом по сравнению с конкурентоспособными термохимическими технологиями.

Однако формирование равномерного распределения температурных полей и полей термических напряжений в процессе поверхностной индукционной закалки деталей сложной геометрической формы с угловыми зонами представляет собой сложную и до сих пор не решенную задачу. Исследование технологического процесса поверхностной индукционной закалки, которое условно может быть разделено на стадии нагрева и охлаждения, удобно провести на основе численной конечно-элементной модели, разработанной в наиболее универсальном мультифизичном ППП ANSYS. Подробное описание процессов, происходящих при индукционном нагреве, представлено в статьях [1, 2].

 В данной статье основное внимание уделяется анализу физических процессов, происходящих при поверхностной закалке с применением высокоинтенсивного способа охлаждения детали. Такой способ позволяет охлаждать заготовки со скоростью не ниже критической Vкр, чтобы предотвратить перлитное превращение и сформировать мартенситную микроструктуру металла. Процесс охлаждения представляет собой интенсивное распыление воды на поверхность заготовки при помощи управляемых разбрызгивающих устройств. Возможность управления позволяет установить один из специальных режимов контроля охлаждения путем регулировки интенсивности подачи воды. По сравнению с растворами полимеров охлаждение распылением водяной струи демонстрирует повышенную экологическую безопасность, поскольку можно отказаться от добавок и других химикатов [3].

Основой целью исследования является моделирование и анализ распределений температуры и термических напряжений при поверхностной закалке заготовок L-образной формы в ППП ANSYS Mechanical APDL. Модель стадии нагрева ориентирована на дальнейшее использование в численной процедуре параметрической оптимизации конструктивных и режимных параметров индукционной установки с целью получения максимально возможной равномерности температурного распределения, от которой существенно зависит качество микроструктуры закаливаемого слоя. Модель стадии охлаждения может быть интегрирована в численную процедуру решения задачи оптимального управления спрейерным устройством, которое позволяет контролировать закалку путем регулировки интенсивности подачи воды.

Моделирование стадии индукционного нагрева

Для исследования процессов индукционного нагрева под поверхностную закалку была разработана двумерная численная осесимметричная модель системы «индуктор – заготовка», включающая взаимосвязанный расчет электромагнитной, тепловой и упругопластической задач в ППП ANSYS Mechanical APDL [4, 5]. Общий вид алгоритма моделирования представлен на рис. 1.

 

Рис. 1. Алгоритм моделирования стадии индукционного нагрева

 

Решение вышеописанных взаимосвязанных задач реализовано при использовании специальных типов конечных элементов для двумерного расчета из библиотеки ANSYS Mechanical APDL: PLANE13, PLANE55, и PLANE182(42) соответственно. В качестве исходных данных задаются свойства материалов в виде табличных нелинейных зависимостей: для электромагнитного анализа – магнитная проницаемость и удельное сопротивление; для теплового анализа – удельная теплоемкость, коэффициент линейного температурного расширения, коэффициент теплопроводности; для прочностного анализа – модуль Юнга, коэффициент Пуассона и билинейная модель пластичности [6, 7]. Затем задается геометрия заготовки сложной формы, типовой квадратной формы витков индуктора, окружающего воздушного пространства, присваиваются свойства отдельным областям и накладывается конечно-элементная сетка. На рис. 2, а представлено схематичное осесимметричное изображение геометрии стальной заготовки сложной геометрической формы, на котором штриховыми линиями указана область поверхностного закаливаемого слоя. На рис. 2, б изображена геометрия системы «индуктор – заготовка» в ANSYS Mechanical APDL c наложенной конечно-элементной сеткой.

 

Рис. 2. Двумерная осесимметричная модель системы «индукционный нагреватель – стальная заготовка»: а – общий вид геометрии заготовки; б – параметры геометрии витков индуктора; в – наложенная конечно-элементная сетка в ППП ANSYS Mechanical APDL

 

На следующем этапе производится настройка опций решателя, задаются начальные и граничные условия. Для исследования технологического процесса индукционной закалки выбраны 18 точек, расположенных в закаливаемом слое: 3 точки – вдоль оси симметрии, 8 точек – на границе закаливаемого слоя, 7 точек –на поверхности заготовки. На рис. 3, а и 3, б слева представлены трехмерные изображения заготовки с указанием расположения выбранных точек. Соответствующие результаты расчета температур в указанных точках представлены на рис. 3, а и 3, б справа.

 

Рис. 3. Температуры закаливаемого слоя: а – зависимость изменения температур от времени на поверхности заготовки; б – зависимость изменения температур от времени вдоль границы закаливаемого слоя и центре заготовки; в – распределение вдоль границы в конце стадии индукционного нагрева

 

Анализ скоростного нагрева выявил, что квадратная форма витков не позволяет достичь требуемой температуры  в области точек Т4 и Т12 и сформировать однородную аустенитную микроструктуру под закалку. При этом область расположения точек Т5 и Т13 перегревается на протяжении всего процесса. Решение демонстрирует очень быстрый нагрев точки Т5 до температуры свыше 1100 °С (приблизительно за 2 секунды), при этом температура в точке Т4 за то же время достигает только 500 °С, за счет чего образуется недопустимый температурный градиент. На рис. 3, в представлено конечное температурное распределение вдоль границы закаливаемого слоя. Перепад температур составил приблизительно 440 °С.

Исходные данные для моделирования стадии нагрева представлены в табл. 1.

 

Таблица 1. Исходные данные для моделирования стадии нагрева

Наименование параметра

Значение

1

Начальная температура заготовки Ta, °С

20

2

Длительность стадии индукционного нагрева τнагр, сек

10

3

Требуемая температура нагрева вдоль границы закаливаемого слоя T*, °С

900

4

Высота заготовки H, мм

50

5

Длина границы закаливаемого слоя стальной заготовки L, мм

55.8

6

Радиус заготовки R1, мм

40

7

Радиус заготовки R2, мм

48

8

Глубина поверхностного закаливаемого слоя C, мм

1,4–2

9

Типовая углеродистая сталь (марка)

40

10

Сторона витка P1, мм

12

11

Расстояние от витков до заготовки P4, мм

2

12

Расстояние между витками № 1 и № 2 P2, мм

2

13

Расстояние между горизонтальной поверхностью выступа и витком № 3 P3, мм

8,17

14

Cила тока витка № 1, 3 I1, А

1432

15

Cила тока витка № 2 I2, А

1832

16

Частота тока источника питания f, кГц

50

17

Шаг расчета модели n, сек

0,2

 

Одним из способов устранения описанных недопустимых неравномерностей температурного распределения является решение задачи оптимального проектирования витка индуктора сложной геометрической формы, частично повторяющей контуры заготовки в угловой зоне [8]. В связи с тем, что в обрабатываемом слое наблюдаются локальные перегревы, при постановке задачи необходимо учесть фазовые ограничения на максимально допустимые значения температуры поверхности Тдоп. Решение соответствующей задачи оптимального проектирования индукционной заготовки представлено в работах [1, 2, 9, 15–17].

При исследовании стадии индукционного нагрева необходимо проанализировать эквивалентные напряжения, для расчета которых применяется критерий наибольшей удельной потенциальной энергии формоизменения. Данный критерий Губера – Мизеса хорошо согласуется с экспериментальными данными, полученными для изотропных материалов. Функционал ANSYS позволяет моделировать напряжения на основе теории малых упругопластических деформаций при помощи билинейной модели пластичности (BKIN) с изотропным упрочнением и не зависящими от температуры параметрами [10, 14]. Результаты расчета эквивалентных напряжений по Мизесу в процессе индукционного нагрева представлены на рис. 4.

 

Рис. 4. Термические напряжение в процессе нагрева: а – эквивалентные напряжения по Мизесу; б – компоненты напряжений по сечению заготовки в точках σ3 и σ7

 

В начале индукционного нагрева наибольшие значения растягивающих напряжений наблюдаются на поверхностном закаливаемом слое заготовки в точках σ1σ6, которые со временем начинают уменьшаться. Напротив, напряжения центра в точке σ7 становятся сжимающими и со временем увеличиваются до значений, сопоставимых с теми, что были на поверхностном слое в начале нагрева.

Для более подробного исследования необходимо рассмотреть осевые компоненты напряжений. На рис. 6 выделены две точки (т. σ3 и т. σ7), в которых проявились наибольшие значения термических напряжений, приближающиеся к пределу текучести материала, который составляет не менее 335 МПа по ГОСТ 1050–2013 для стали 40. В точке σ3, находящейся в угловой области на поверхности заготовки, наибольшие термические напряжения наблюдаются в начале процесса нагрева, а в расположенной в центре точке σ7 – к концу процесса. Приближение значений к пределу текучести подтверждает необходимость учета технологических ограничений на термонапряжения при решении задачи оптимального проектирования и управления процессами индукционного нагрева. Максимальное значение термонапряжений в ходе всего процесса не превышает предела прочности для стали 40, который составляет не менее 570 МПа. Итоговый вектор суммы компонент XY деформаций геометрии заготовки составил около 0,4 мм, что составляет 1 % радиуса R1=40 мм и 0,83 % радиуса R2=48 мм.

Моделирование стадии охлаждения

Одна из главных целей поверхностной закалки – формирование мартенситной микроструктуры, которая достигается высокой скоростью охлаждения детали. Для этого после поверхностного индукционного нагрева заготовок переходят к следующей стадии – высокоскоростному охлаждению закалочной средой, например водяным душем. При таком способе подачи охлаждающей воды теплообмен с окружающей средой происходит значительно интенсивнее благодаря разрушению паровой пленки из-за кинетической энергии струй и большой скорости перемещения воды относительно охлаждаемой поверхности. Интенсивность охлаждения может изменяться в широких пределах за счет изменения количества смеси и ее скорости истечения из разбрызгивающего устройства. Однако у такого способа охлаждения есть очевидный недостаток – необходимость периодической очистки отверстий форсунок, поэтому вода должна быть достаточно чистой, т. к. окалина и другие механические примеси могут очень быстро засорить отверстия [11, 12].

Для исследования стадии охлаждения под поверхностную закалку в ППП ANSYS Mechanical APDL разработана двумерная численная осесимметричная модель системы «устройство подачи охлаждающей среды – заготовка», включающая взаимосвязанный расчет тепловой и упругопластической задач (рис. 5).

 

Рис. 5. Алгоритм моделирования стадии охлаждения

 

В работе рассматривается охлаждение путем распыления воды через разбрызгивающее устройство. Для моделирования процесса охлаждения используется эквивалентный поверхностный коэффициент теплоотдачи, который изменяется в зависимости от температуры поверхности заготовки и параметра расхода воды М, м3/(с∙м2), являющегося постоянной величиной, равной отношению количества охлаждающей жидкости (м3/сек) к площади закаливаемой заготовки (м2) [12]. Общий вид системы охлаждения и применяемая при моделировании зависимость эквивалентного коэффициента поверхностной теплоотдачи αэкв от температуры поверхности охлаждения представлены на рис. 6.

 

Рис. 6. Моделирование стадии охлаждения: а – схема спрейерной системы подачи охлаждающей среды на поверхность металла; б – расчетная зависимость коэффициентов теплоотдачи от температуры с поверхности заготовки при  м3/(с∙м2)

 

Исходные данные для моделирования стадии охлаждения представлены в табл. 2.

 

Таблица 2. Исходные данные для моделирования стадии охлаждения

Наименование параметра

Значение

1

Длительность стадии охлаждения τохл, сек

5

2

Шаг расчета модели m, сек

0,05

3

Температура начала мартенситного превращения Tм, °С

≈345

4

Критическая скорость охлаждения для стали 40, °С/сек

≈520

 

В процессе закалки скорость охлаждения заготовки не должна опускаться ниже критической Vкр, т. к. в интервале температур перлитного превращения (от точки А1 (724 °С) до температур минимальной устойчивости переохлажденного аустенита углеродистых и низколегированных сталей (500–550 °С)) закалочная среда должна обеспечить такую интенсивность охлаждения, при которой будет предотвращен перлитный распад аустенита и температура заготовки опустится ниже обозначенного интервала, в котором устойчивость аустенита возрастает и замедление охлаждения уже не опасно. В противном случае аустенит будет претерпевать перлитное превращение, что снизит качество закалки [11, 13]. В начале стадии охлаждения наибольший температурный перепад возникает между поверхностью и центром, особенно в области выступа. В ходе процесса наблюдается снижение температуры поверхностного слоя заготовки и прогрев центральной части заготовки путем теплопроводности (рис. 7).

 

Рис. 7. Температуры в точках вдоль границы закаливаемого слоя и в центре заготовки в процессе охлаждения

 

Выбор интервала температур мартенситного превращения производится в зависимости от материала. В основном используются диапазоны от 350 до
100 °С или от 300 до 200 °С. В этих интервалах происходит окончательное превращение аустенита в мартенсит, которое сопровождается возникновением высоких внутренних напряжений и опасностью появления закалочных микротрещин [11].

Вследствие высокой пластичности стали деформации при высоких температурах оказываются упругопластическими и поверхностные слои растягиваются относительно центральных. Охлаждение поверхностных слоев приводит к неодинаковому тепловому сжатию, а следовательно, к возникновению термических деформаций и напряжений. Поскольку охлаждаемое тело остается сплошным, то на поверхности в конце охлаждения наблюдаются сжимающие напряжения, а в центре – растягивающие. По окончании процесса охлаждения изначально возникшие временные напряжения затем формируют остаточные напряжения, имеющие противоположный знак [11, 13, 14]. На рис. 8 представлены результаты расчета термических напряжений в процессе быстрого охлаждения.

 

Рис. 8. Термические напряжение на стадии охлаждения: а – эквивалентные напряжения по Мизесу; б – компоненты напряжений в точке σ6 на поверхности; в – компоненты напряжений в точке σ7 в центре

 

Как показывает анализ результатов моделирования, в начале интенсивного охлаждения наибольшие значения термических напряжений образуются на поверхности заготовки, а затем их величина постепенно уменьшается. Со временем в центре заготовки на оси Y наблюдается повышение температуры от прогретой за счет теплопроводности поверхности и одновременное снижение напряжений, после чего все поверхности контакта охлаждаются до температуры воды и происходит перераспределение внутренних напряжений, что влияет на свойства конечного продукта.

Во избежание формирования чрезмерной хрупкости металла и для получения достаточной его твердости необходимо сформулировать и решить задачу поиска такого управления подачей воды спрейерным устройством, которое создает оптимальные физические условия теплообмена с поверхностью заготовки, позволяющие сформировать мартенситную микроструктуру металла с требуемыми свойствами.

Заключение

Рассмотрена проблемно-ориентированная численная нелинейная модель технологического процесса поверхностной индукционной закалки стальной заготовки сложной геометрической формы с угловой зоной. Разработанные и реализованные в ППП ANSYS алгоритмы моделирования позволили получить пространcтвенно-временные распределения внутренних электромагнитных источников тепла, температур и термонапряжений на стадиях нагрева заготовки в индукционном нагревателе и ее охлаждения с помощью водяного душа.

Анализ результатов моделирования позволил выявить области значений температур и термических напряжений, превышающих или приближающихся к предельно допустимым согласно технологическим требованиям. Показано, что представленная конфигурация индуктора с квадратной формой витков не позволяет одновременно прогреть угловую зону и не перегреть выступ, что подтверждает необходимость перехода к более сложной геометрической форме витка индуктора, находящегося в области выступа стальной заготовки L-образной геометрической формы.

Представленная модель стадии индукционного нагрева, позволяющая проанализировать влияние вариаций основных технологических параметров процесса индукционного нагрева на динамику изменения тепловых полей и полей термонапряжений, может быть интегрирована в численную процедуру оптимального проектирования сложной формы витков индуктора и оптимизации его режимных параметров, обеспечивающую требуемую равномерность температурного поля при отсутствии недопустимых перегревов.

Представленная модель стадии охлаждения может быть интегрирована в численную процедуру для поиска алгоритма управления спрейерным устройством, который обеспечит оптимальную интенсивность подачи воды с целью получения максимального качества мартенситной микроструктуры металла в конце процесса закалки.

×

Об авторах

Алексей Владимирович Павлушин

Самарский государственный технический университет

Автор, ответственный за переписку.
Email: Alex1995i79@gmail.com
ORCID iD: 0000-0003-1645-5803

инженер НИС кафедры «Автоматика и управление в технических системах», аспирант

Россия, 443100, г. Самара, ул. Молодогвардейская, 244

Список литературы

  1. Pleshivtseva Y., Pavlushin A., Popov A. Optimal Design of Inductor Coils for Surface Hardening // 2022 International Conference on Industrial Engineering, Applications and Manufacturing (ICIEAM). IEEE, 2022. Pp. 680–685.
  2. Павлушин А.В. Оптимизация конструктивных и режимных параметров индукционной нагревательной установки для закалки // Вестник СамГТУ. Серия: Технические науки. 2021. Т. 29. № 3 (71). С. 38–51.
  3. Арзамасов Б.Н. и др. Материаловедение: Учебник для высших технических учебных заведений. М.: Машиностроение, 1986. 386 с.
  4. Engineering Simulation & 3D Design Software | ANSYS [Электронный ресурс]. 2022. URL: https://www.ansys.com/ (дата обращения: 30.05.2022).
  5. Rudnev V., Loveless D., Cook R.L. Handbook of induction heating. CRC press, 2017.
  6. Сорокин В.Г. и др. Марочник сталей и сплавов. М.: Машиностроение, 1989. Т. 640. С. 36.
  7. Драгунов Ю.Г. Марочник сталей и сплавов. Под ред. А.С. Зубченко. 2016.
  8. Рапопорт Э.Я., Плешивцева Ю.Э. Оптимальное управление температурными режимами индукционного нагрева. М.: Наука, 2012.
  9. Pleshivtseva Y., Popov A., Pavlushin A. Software Solution for Optimal Design of Inductor for Surface Hardening // 2021 3rd International Conference on Control Systems, Mathematical Modeling, Automation and Energy Efficiency (SUMMA). IEEE, 2021. Pp. 646–649.
  10. Коршиков С.Е. Моделирование полей упругих деформаций в процессе индукционного нагрева вращающихся заготовок // Информатика, телекоммуникации и управление. 2012. № 1 (140). С. 94–98.
  11. Башкин Ю.А., Ушаков Б.К., Секей А.Г. Технология термической обработки стали: Учебник для студентов вузов. Металлургия, 1986.
  12. Dossett J., Totten G.E. ASM Handbook, Volume 4A: Steel Heat Treating Fundamentals and Processes, ASM Intern. 2013.
  13. Анализ диаграммы изотермического распада переохлажденного аустенита: метод. указания к практ. и лаб. работам по дисциплинам «Основы термической обработки», «Материаловедение» для студентов всех форм обучения по направлению подготовки 150400 «Металлургия» / Т.В. Нуждина, М.Н. Чеэрова, Т.В. Комарова. Н. Новгород: НГТУ, 2014. 34 с.
  14. 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-3-030-39225-3. Part of: Springer Professional Wirtschaft+Technik", Springer Professional "Technik".
  15. Плешивцева Ю.Э., Попов А.В., Попова М.А., Деревянов М.Ю. Оптимальное проектирование индуктора для поверхностной закалки цилиндрических заготовок на основе численной двумерной модели // Вестник Астраханского государственного технического университета. Серия: Управление, вычислительная техника и информатика. 2019. № 1. С. 40–50.
  16. 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.
  17. Коршиков С.Е. Оптимальное управление температурными режимами индукционного нагрева цилиндрических слитков с учетом технологических ограничений: дис. … канд. техн. наук. Самара: Сам. гос. техн. ун-т, 2015. 157 с.

Дополнительные файлы

Доп. файлы
Действие
1. JATS XML
2. Рис. 1. Алгоритм моделирования стадии индукционного нагрева

Скачать (124KB)
3. Рис. 2. Двумерная осесимметричная модель системы «индукционный нагреватель – стальная заготовка»: а – общий вид геометрии заготовки; б – параметры геометрии витков индуктора; в – наложенная конечно-элементная сетка в ППП ANSYS Mechanical APDL

Скачать (181KB)
4. Рис. 3. Температуры закаливаемого слоя: а – зависимость изменения температур от времени на поверхности заготовки; б – зависимость изменения температур от времени вдоль границы закаливаемого слоя и центре заготовки; в – распределение вдоль границы в конце стадии индукционного нагрева

Скачать (300KB)
5. Рис. 4. Термические напряжение в процессе нагрева: а – эквивалентные напряжения по Мизесу; б – компоненты напряжений по сечению заготовки в точках σ3 и σ7

Скачать (259KB)
6. Рис. 5. Алгоритм моделирования стадии охлаждения

Скачать (130KB)
7. Рис. 6. Моделирование стадии охлаждения: а – схема спрейерной системы подачи охлаждающей среды на поверхность металла; б – расчетная зависимость коэффициентов теплоотдачи от температуры с поверхности заготовки при  м3/(с∙м2)

Скачать (106KB)
8. Рис. 7. Температуры в точках вдоль границы закаливаемого слоя и в центре заготовки в процессе охлаждения

Скачать (123KB)
9. Рис. 8. Термические напряжение на стадии охлаждения: а – эквивалентные напряжения по Мизесу; б – компоненты напряжений в точке σ6 на поверхности; в – компоненты напряжений в точке σ7 в центре

Скачать (208KB)

© Самарский государственный технический университет, 2022

Creative Commons License
Эта статья доступна по лицензии Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International License.

Данный сайт использует cookie-файлы

Продолжая использовать наш сайт, вы даете согласие на обработку файлов cookie, которые обеспечивают правильную работу сайта.

О куки-файлах