МЕТАЛЛУРГИЯ
УДК 621.791.722:621.762.04
А.И. Рудской, Ю.А. Соколов, В.Н. Копаев
МОДЕЛИРОВАНИЕ ПРОЦЕССА ЭЛЕКТРОННО-ЛУЧЕВОГО СИНТЕЗА ИЗДЕЛИЙ ИЗ ПОРОШКА ТИТАНОВОГО СПЛАВА ВТ-6 С УЧЕТОМ ТЕПЛОТЫ ФАЗОВОГО ПЕРЕХОДА (ДВУХФАЗНАЯ ЗАДАЧА СТЕФАНА)
Л.1. Rudskoi, Yu.A. Sokolov, V.N. Kopaev
MODELLING THE PROCESS OF ELECTRON BEAM SYNTHESIS OF PRODUCTS FROM GRANULES OF TITANIC ALLOY V^6 ALLOWING FOR HEAT CAPACITY OF PHASE TRANSITION (DIPHASE STEFAN TASK)
В статье рассматривается моделирование процесса электронно-лучевого синтеза изделий из порошка титанового сплава ВТ-6 с учетом теплоты фазового перехода. Решение задачи Стефана позволило определить скорости нагрева и охлаждения, глубину проплавления, скорость движения перемещения твердо-жидкой границы вглубь изделия.
СИНТЕЗ ИЗДЕЛИЙ; ЭЛЕКТРОННЫЙ ЛУЧ; МОДЕЛИРОВАНИЕ; ЗАДАЧА СТЕФАНА; ТЕМПЕРАТУРНОЕ ПОЛЕ.
The present article considers modelling the process of electron beam synthesis of products from granules of titanic alloy B^6 allowing for heat capacity of phase transition. It is established that during the electron beam synthesis of products the process of transforming the powder in a liquid state with its subsequent crystallization can be described by the Stephan two-phase task. The principles of calculating an effective mode of processing are shown at the layer-by-layer melting of powder of the BT-6 titanic alloy on the basis of the mathematical model allowing for processes of heat and mass transfer. The solution to the Stefan task has allowed to define speeds of heating and cooling, fusion depth, traverse speed of moving the firm-liquid boundary into the product.
SYNTHESIS OF PRODUCTS; ELECTRON BEAM; MODELLING; STEFAN TASK; TEMPERATURE FIELD.
Несмотря на широкое применение порошковой металлургии в современной промышленности, работы по созданию новых технологий получения таких материалов остаются актуальными [1—3]. Так, для получения поверхностных слоев с заданной структурой и свойствами перспективно использование лазерной [4, 5] и электронно-лучевой обработки деталей, изготовляемых из различных металлических материалов.
В основу технологии синтеза изделий электронным лучом (СИЭЛ) положена операция послойного плавления порошка тугоплавких сплавов. К преимуществам данной технологии следует отнести способность электронов, достигших поверхности обработки, внедряться в материал на некоторую глубину, называемую пробегом электронов. Длина пробега электронов, зависящая от ускоряющего напряжения и плотности материала, может достигать не-
скольких десятков микрон, что повышает твердость и износостойкость синтезируемого слоя. Высокие скорости нагрева слоя порошка, достигающие значений 105—106 °С/с, создают необходимые условия для формирования изделия с мелкодисперсной структурой.
Целью работы является расчет эффективного режима обработки при послойном плавлении порошка титанового сплава ВТ-6 на базе математической модели, построенной с учетом процессов тепло- и массопереноса.
Постановка задачи
Использование концентрированного потока электронов при процессе СИЭЛ позволяет управлять глубиной проплавления и формированием структуры изделия. К числу важнейших вопросов, которые необходимо решить при реализации технологии СИЭЛ, относится расчет режимов обработки, использующий математическую модель, построенную с учетом процессов тепло- и массопереноса при формировании изделия. Для анализа глубины проплавления и динамики перемещения твердо-жидкой границы вследствие периодического воздействия электронного луча на поверхность необходимо решить задачу Стефана.
Температурное поле при нагреве и теплота фазового перехода, выделяющаяся на границе плавления (кристаллизации), описываются двухфазной задачей Стефана [6, 7]:
срдТ = ^ ( ЕгаёТ) + д - Ъ,ЬУп, (1)
где с — теплоемкость материала; X — коэффициент теплопроводности материала; р — плотность материала; Т — температура; д — плотность теплового источника; 54 — поверхностная 5-функция; Уп — скорость движения границы фазового перехода по нормали; Ь — энтальпия фазового перехода.
Температурную зависимость поглощательной способности металлов можно аппроксимировать полиномом первой степени, коэффициент теплопроводности и объемную теплоемкость — полиномом третьей степени:
А(Т) = а + ЬТ; (2)
Х(Т) = X! + Х2Т + Х3Т2 + Х4Т3; (3) С(Т) = с1 + С2Т+ С3Т2 + С4Т3. (4)
В работе [8] приводятся полиномы для расчета теплопроводности, теплоемкости и плотности титана в зависимости от температуры:
Х(Т = 9,13 + 0,00323 Т - 3,94 10~5Т2 +
+ 1,17 10-8Т3, Вт/(м-К); (5)
C(T = 404,2 + 0,21 Т, Дж/(кг-К); (6) р(Т) = 4574—0,194Т, кг/м3. (7)
Пористость верхнего слоя порошка следует учитывать при операциях предварительного нагрева формируемого изделия, при плавлении порошка — только в начальный малый момент времени.
Построение модели
В качестве инструмента моделирования выбран программный комплекс ANSYS Workbench, пространственный режим — 3D, режим расчета— динамический, тепловой режим — неадиабатный.
Исходные данные: диаметр порошка — 50-10-6 м; высота подложки (изделия) — 0,1 м; температурный интервал между температурами солидус и ликвидус — 100 °С; температура со-лидус Tsol = 1570 °C; температура ликвидус Tliq = = 1670 °C. Подложка моделируется тепловой массой со значением теплового захвата 3,4069-10-2 Дж/°С.
Скорость сканирования электронного луча диаметром 1 мм по поверхности подложки составляет 4530 мм/с, шаг смещения луча — 1 мм. Размер сечения синтезируемого изделия толщиной 0,1 м составлял 70x70 мм. Для моделирования был выбран элементарный объем размером 100,050x1x1 мм.
Источник тепла перемещается по подложке со скоростью Уск по определенной траектории. Способ сканирования электронного луча по поверхности показан на рис. 1, на котором изображен объем элементарной ячейки.
Значения размерных и скоростных параметров процесса электронно-лучевого синтеза образцов из порошка сплава ВТ-6 приведены в таблице.
При моделировании были исследованы два режима электронно-лучевой обработки при постоянном ускоряющем напряжении 60 кВ:
ток луча 15 мА (мощность 900 Вт);
ток луча 20 мА (мощность 1200 Вт).
Нагрев каждой элементарной ячейки происходит за время, равное времени перемещения электронного луча над этой ячейкой. Для нашего случая это эквивалентно приложения тепловой нагрузки треугольной формы:
t = ^пятна/^ск. (8)
Значения размерных и скоростных параметров процесса синтеза образцов из порошка титанового сплава ВТ-6 лучом мощностью 1200 Вт
Start
Finish
-в-
70 мм
70 мм
100,
05 мм
Параметр сканирования Обозначение Единица измерения Значение
Площадь изделия II «2 мм2 4900
Площадь, перекрываемая лучом в секунду Кк do мм2/с 4530
Период сканирования A/VCKd0 с 1,081678
Частота сканирования f Гц 0,92449
Число проходов N 32,35714
Подведенное тепло Qj Дж 42000
Плотность подведенного тепла к изделию за время процесса Jonep Дж/мм2 8,571429
Время сканирования участка 1х1 мм при прямоугольном характере воздействия ¿1x1 = do/VCK с 0,000221
Из таблицы следует, что время сканирования данного участка при прямоугольном воздействии ¿1х1 равно 221 мкс. При треугольном воздействии время увеличится в два раза и составляет 442 мкс.
Граничные и начальные условия математической модели (ММ) показаны на рис. 2:
области А и В (верхняя и нижняя поверхности изделия) — радиационное тепловое излучение на экран вакуумной камеры;
область С — переменная начальная температура по глубине синтезируемого изделия;
область В, Е — начальная температура предварительного нагрева(1000 °С);
область Г — начальная температура подложки; область О — тепловой поток, приложенный к элементарному объему;
область Н — точечная тепловая масса, моделирующая подложку.
Охлаждение на границах «изделие — окружающая среда» задается через лучистый теплообмен (передача тепла посредством электромагнитных волн) поверхности спекания/плавления с внутренней оболочкой вакуумной рабочей
Рис. 2. Граничные и начальные условия модели
камеры. При расчете теплообмена излучением необходимо учитывать потоки тепла в областях, попадающих на экраны, установленные в вакуумной камере установки. С поверхности подложки тепло теряется излучением по закону Стефана — Больцмана.
В верхней области формируемого изделия задаем граничные условия второго рода (условия Неймана), соответствующие заданию на границе распределенного теплового потока. Для уравнения теплопроводности в изотропной среде оно записывается в виде
X дТ/дп = q (t).
(9)
Начальные условия:
приведенная степень черноты излучатель-ных потоков — епр = 0,333;
температура окружающей среды и тепловых экранов — 20 °С;
распределение температуры в зависимости от координаты z:
800 + (z/0,1)0,3 200 при 0< z < 0,1;
1000 при 0,1 < z < 0,10005 м. (10) Начальная температурная нагрузка на уже сплавленную структуру, как функция от Z, представлена на рис. 3.
Рис. 3. Начальная температурная нагрузка
На рис. 4 представлена сетка конечных элементов модели. Как видно из рисунка, сетка имеет сгущение около поверхности нагрева. Размер элемента составляет 1 мкм около верхней поверхности и 8,7 мм в районе подложки. Такое различие в размерах элементов обусловлено минимизацией количества конечных элементов (КЭ), что важно при решении нестационарных задач, и возможностью поймать движение фазового перехода как по времени, так и в пространстве.
Типы конечных элементов, используемых в модели:
LINK 33 — линейный элемент, моделирующий теплопроводность;
SURF152 — поверхностный элемент, позволяющий отобразить нагрузку в виде излучения.
На элементарный объем в модели проведено несколько воздействий треугольной формы. Динамика изменения температур на поверхности порошка (обозначение T_UP), на границе слоев (обозначение T_GR) в интервале времени первого воздействия электронного луча мощностью соответственно 1200 и 900 Вт на ЭЯ приведена на рис. 5.
a)
Т, 'С 2120 2040 1960 1080 1800 1128 1640 1560 1480 1400
T_UP
1
\
1 \
1 \ Tfiq Tsol
J- 1 /
1
1 m 1 -\T_OR __
11
б)
1.2 1.8 2.4 t, 10-3 с
1.5 2.1 2 Л
Т, 'С
1920 1840 1760 16S0 1600 1520 1440 1360 1280 1200
A T_UP
1 niq r~ Tsol
j \ /
"H Г'"
i 1
I / 4k
/ **—
/ T_GR
/
1.2 1.6 t, 10-3 с
1 1.4 1.8
Рис. 5. Динамика изменения температур на верхней и нижней границах слоя порошка в результате первого воздействия электронного луча мощностью 1200 (а) и 900 (б) Вт
Наблюдаются уменьшения скорости изменения температур при прохождении фазового перехода. На рис. 5 хорошо видно влияние скрытой теплоты плавления на процесс нагрева и охлаждения: на графике представлены участки изменения скорости нагрева и охлаждения. При нагреве поглощается скрытая теплота, при охлаждении — выделяется.
При мощности луча 1200 Вт скорость нагрева верхней поверхности слоя порошка в диапазоне температур от 1400 до 2120 °С равна 3,6-106 °С/с, для нижней поверхности в диапазоне температур от 1400 до 1600 °С скорость нагрева равна 1-106 °С/с.
При мощности луча 900 Вт скорость нагрева верхней поверхности слоя порошка в диапазоне температур от 1200 до 1900 °С равна 3,5 106 °С/с, для нижней поверхности в диапазоне температур от 1200 до 1500 °С она равна 1,5 106 °С/с.
Скорость охлаждения в диапазоне температур от ликвидус до солидус при мощности электронного луча 1200 Вт для верхней границы слоя равна 1,4 105 °С/с, при мощности луча 900 Вт — 3-105 °С/с.
Время нагрева верхнего слоя порошка в результате воздействия электронного луча мощностью 900 Вт и 1200 Вт соответственно составляет 0,4-10-3 и 0,2-10-3 с, время кристаллизации — 0,67-10-3 с и 0,65-10-3 с. На рис. 5 видно, что скорости охлаждения верхнего и нижнего слоя порошка практически совпадают.
Для анализа динамики перемещения твердо-жидкой границы задан малый интегральный шаг по времени, равный 2 мкс. Глубина фазового перехода после первого воздействия показана на рис. 6 (а — для луча мощностью 1200 Вт, б— 900 Вт).
На флажке указана температура поверхности соприкосновения порошка с уже сплав-
Рис. 6. Глубина фазового перехода после первого воздействия электронного луча
ленной структурой. Особенности тепловых процессов связаны с эффектами накопления тепла и движения твердо-жидкой границы под воздействием предыдущих периодических импульсов.
Таким образом, решение задачи Стефана позволяет определить важнейшие параметры процесса СИЭЛ (скорости нагрева и охлаждения слоя порошка, время плавления и кристаллизации порошка, перемещение твердо-жидкой границы в слое при электронно-лучевом воздействии). При оптимизации процесса синтеза изделия модель позволяет определить значения составляющих вектора варьируемых параметров (ток 1л луча, диаметр пятна на мишени, скорость Уск сканирования луча, способ сканирования луча, шаг 5см смещения луча при сканировании), которые обеспечивают надежное сплавление текущего слоя порошка с предыдущим, формируют необходимую структуру изделия.
Сложность физико-химических процессов, фазовые переходы при воздействии на металл электронным лучом предопределяют необходимость создания математической модели синтеза изделия.
В ходе процесса СИЭЛ важно размерное формирование синтезируемого слоя порошка,
обеспечение управляемого сплавления текущего слоя с предыдущим.
Математически процесс превращения порошка в жидкое состояние с последующей их кристаллизацией описывается двухфазной задачей Стефана.
Модель учитывает нелинейную зависимость теплофизических характеристик материала порошка от температуры. На графиках отражено влияние скрытой теплоты плавления на процесс нагрева и охлаждения.
Решение задачи Стефана позволило определить скорости нагрева и охлаждения, глубину проплавления, скорость движения перемещения твердо-жидкой границы вглубь изделия при мощности луча 900 и 1200 Вт.
Показана динамика изменения температур на верхней и нижней поверхностях слоя порошка в результате одиннадцати воздействий электронного луча мощностью 900 и 1200 Вт.
Моделирование в среде ANSYS APDL позволяет более гибко реагировать на вариации входных параметров: мощность луча, скорость сканирования, смещение луча и т. д.
Исследования КЭ-модели с мощностью 1200 и 900 Вт показали, что изменения мощности электронного луча до 900 Вт практически не гарантирует проплавления границы верхнего слоя порошка с изделием.
СПИСОК ЛИТЕРАТУРЫ
1. Рудской А.И., Кондратьев С.Ю., Кокорин В.Н., Сизов Н.А. Исследование процесса уплотнения при ультразвуковом воздействии на увлажненную порошковую среду // Научно-технические ведомости СПбГПУ. 2013. № 3 (178). С. 148-155.
2. Рудской А.И., Кокорин В.Н., Кондратьев С.Ю., Филимонов В.И., Кокорин А.В. Прессование гетеро-фазных увлажненных железных порошков при использовании метода интенсивного уплотнения // Наукоемкие технологии в машиностроении. 2013. № 5 (23). С. 13-20.
3. Рудской А.И., Кондратьев С.Ю, Кокорин В.Н. Прессование гетерофазных увлажненных порошковых металлических смесей для повышения качества высокоплотных заготовок с использованием метода интенсивного уплотнения // Справочник. Инженерный журнал с приложением. 2011. № 6. С. 12-16.
4. Kondrat'ev S.Y., Gorynin V.I., Popov V.O. Optimization of the parameters of the surface-hardened layer in laser quenching of components // Welding International. August 2012. Vol. 26. № 8. P. 629-632.
5. Кондратьев С.Ю., Горынин В.И., Попов В.О. Оптимизация параметров поверхностно-упрочненного слоя при лазерной закалке деталей // Сварочное производство. 2011. № 3. С. 11—15.
6. Самарский А.А., Вабищевич П.Н. Вычислительная теплопередача. М.: Едиториал УРСС, 2003. 784 с.
7. Веденов А.А., Гладуш Г.Г. Физические процессы при ЛО материалов. М.: Энергоатомиздат, 1985. 208 с.
8. Углов А.А., Смуров И.Ю., Лашин А.М., Гуськов
А.Г. Моделирование теплофизических процессов импульсного лазерного воздействия на металлы. М.: Наука, 1991. 288 с.
REFERENCES
1. Rudskoy A.I., Kondratyev S.Iu., Kokorin V.N., Sizov N.A. Issledovaniye protsessa uplotneniya pri ultra-zvukovom vozdeystvii na uvlazhnennuyu poroshkovuyu sredu. Nauchno-tekhnicheskiye vedomosti SPbGPU. 2013. № 3 (178). S. 148-155. (rus.)
2. Rudskoy A.I., Kokorin V.N., Kondratyev S.Iu., Filimonov V.I., Kokorin A.V. Pressovaniye geterofaznykh uvlazhnennykh zheleznykh poroshkov pri ispolzovanii metoda intensivnogo uplotneniya. Naukoyemkiye tekh-nologii vmashinostroyenii. 2013. № 5 (23). S. 13-20. (rus.)
3. Rudskoy A.I., Kondratyev S.Iu, Kokorin V.N. Pressovaniye geterofaznykh uvlazhnennykh poroshkovykh me-tallicheskikh smesey dlya povysheniya kachestva vysoko-plotnykh zagotovok s ispolzovaniyem metoda intensivnogo uplotneniya. Spravochnik. Inzhenernyy zhurnals prilozheni-yem. 2011. № 6. S. 12-16. (rus.)
4. Kondrat'ev S.Y., Gorynin V.I., Popov V.O. Optimization of the parameters of the surface-hardened layer in laser quenching of components. Welding International. August 2012. Vol. 26. № 8. P. 629-632. (rus.)
5. Kondratyev S.Yu., Gorynin V.I., Popov V.O. Op-timizatsiya parametrov poverkhnostno-uprochnennogo sloya pri lazernoy zakalke detaley. Svarochnoye proizvod-stvo. 2011. № 3. S. 11-15. (rus.)
6. Samarskiy A.A., Vabishchevich P.N. Vychislitelnaya teploperedacha. M.: Yeditorial URSS, 2003. 784 s. (rus.)
7. Vedenov A.A., Gladush G.G. Fizicheskiye protsessy pri LO materialov. M.: Energoatomizdat, 1985. 208 s. (rus.)
8. Uglov A.A., Smurov I.Iu., Lashin A.M., Guskov A.G. Modelirovaniye teplofizicheskikh protsessov impul-snogo lazernogo vozdeystviya na metally. M.: Nauka, 1991. 288 s. (rus.)
СВЕДЕНИЯ ОБ АВТОРАХ
РУДСКОЙ Андрей Иванович — доктор технических наук ректор Санкт- Петербургского государственного политехнического университета. 195251, Санкт-Петербург, Политехническая ул., 29. E-mail: rector@spbstu.ru
СОКОЛОВ Юрий Алексеевич — кандидат технических наук главный конструктор ООО «НТК Электромеханика». 172386, ул. Заводское шоссе, 2, г. Ржев, Россия. E-mail: s5577@inbox.ru
КОПАЕВ Владимир Николаевич — ведущий инженер ООО «НТК Электромеханика». 172386, ул. Заводское шоссе, 2, г. Ржев, Россия. E-mail: vkopaev_ela@mail.ru
AUTHORS
RUDSKOI Andrei I. — St. Petersburg State Polytechnical University. 29, Politechnicheskaya St., St. Petersburg, 195251, Russia. E-mail: rector@spbstu.ru
SOKOLOV Yurii A. — «Electromechanica». 2, Zavodskoe shosse, Rzhev city Tver region, Russia, 172386. E-mail: s5577@inbox.ru
KOPAEV Vladimir. N. — «Electromechanica». 2, Zavodskoe shosse, Rzhev city Tver region, Russia, 172386. E-mail: vkopaev_ela@mail.ru
© Санкт-Петербургский государственный политехнический университет, 2014