CC BY
15
УДК 621.791.722:621.762.04
А.И. Рудской, Ю.А. Соколов, В.Н. Копаев
ОПРЕДЕЛЕНИЕ ТЕМПЕРАТУРНОГО ПОЛЯ ПРИ СКАНИРОВАНИИ ПОВЕРХНОСТИ ЭЛЕКТРОННЫМ ЛУЧОМ В ТЕХНОЛОГИИ СИНТЕЗА ПОРОШКОВЫХ ИЗДЕЛИЙ
A.I. Rudskoi, Yu.A. Sokolov, V.N. Kopaev
DETERMINATION OF TEMPERATURE FIELD AT SURFACE SCANNING BY ELECTRONIC BEAM SYNTHESIS TECHNOLOGIES OF POWDER PRODUCTS
На базе математической модели, описывающей процесс термической обработки, решена задача определения температурного поля при сканировании поверхности подложки электронным лучом. Моделирование выполнено на основе метода конечных элементов.
МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ; СКАНИРОВАНИЕ ПОВЕРХНОСТИ; ЭЛЕКТРОННЫЙ ЛУЧ; НАГРЕВ ПОДЛОЖКИ; ТЕПЛОПЕРЕНОС; ТЕМПЕРАТУРНОЕ ПОЛЕ.
In present article on the basis of the mathematical model describing process of thermal treatment, the problem of definition of a temperature field is solved per scanning of surface's substrate by an electronic beam. Modelling is executed on the basis of finite element method.
MATHEMATICAL MODELING; SURFACE SCANNING; ELECTRONIC BEAM; SUBSTRATE HEAT; CONDUCTION; TEMPERATURE FIELD.
В порошковой металлургии интенсивно разрабатываются новые технологии, позволяющие шире использовать ее возможности по созданию материалов с повышенными эксплуатационными свойствами [1—3]. В то же время в промышленности широко применяются технологии, позволяющие формировать на поверхности различных материалов слои с заданной структурой и свойствами [4, 5]. Одна из наиболее перспективных технологий — это синтез изделий с помощью электронного луча (СИЭЛ), где в основу положена операция последовательного наплав-ления слоев из порошка различных сплавов, успешно зарекомендовавшая себя при изготовлении изделий авиационной и медицинской промышленности [6, 7]. Технологический процесс происходит в вакууме, что позволяет проводить синтез изделий из порошка химически активных металлов, таких, как титан, молибден, вольфрам и др.
В процессе электронно-лучевого синтеза изделия нагрев поверхности осуществляют две
операции: предварительный прогрев подложки; предварительный нагрев поверхности формируемого изделия перед плавлением.
Предварительный прогрев подложки, создающий необходимые начальные тепловые условия для проведения процесса синтеза, осуществляется сканированием электронного луча по поверхности. Для реализации данной операции необходимо выбрать способ сканирования и рассчитать параметры нагрева подложки электронным лучом до заданной температуры Тзад: ток луча, скорость сканирования луча, диаметр луча, шаг смещения луча, размеры сканируемой площади х и у1 (угол отклонения луча). Основная задача на данном этапе — достичь равномерного температурного поля по всей поверхности подложки. Одним из инструментов решения этой задачи служит выбор оптимального способа сканирования поверхности обработки.
Целью работы, о которой идет речь в статье, расчет температурного поля подложки и изделия в процессе СИЭЛ, базирующийся на нелиней-
ной математическом модели, учитывающем зависимость теплофизических параметров от температуры.
Описание задачи. Нагрев подложки толщиной £под, по поверхности которой перемещается электронный луч со скоростью Уск, будем рассматривать в движущейся прямоугольной системе координат (ОХУЯ), помещенной в центр пятна нагрева (рис. 1).
Процесс распространения тепла вследствие действия сканирующего луча от кругового источника с мощностью д, перемещающегося по поверхности со скоростью Уск, описывается следующим дифференциальным уравнением теплопроводности [8]:
с(Т)
ЭТ э*
д_ Эх
X (Т)
ЭТ
Эх
+-
Эу
X (Т)
ЭТ Эу.
+
Э
Х(Т %
-*ТУх Эх-
ЭТ
ЭТ
- с(Т)Уу — -с(Т)УЕ — + д (Ух*, Уу*, УЛ, (1)
где Т — температура; Ух, Уу, УЕ — проекции вектора скорости быстродвижущегося источника тепла по осям X, У и 2; с(Т) и ЦТ) — теплоемкость и коэффициент теплопроводности материала; д(У^, Уу*, У* — удельная мощность движущегося источника тепла.
При решении уравнения (1) методом конечных элементов используется следующее граничное условие: заданный поток тепла д имеет вид ХдТ/дп = А(Т д (Ух*, Уу*, УА (2) где А(Т — доля поглощенной мощности от мощности пучка поступающей на поверхность мишени, или эффективный КПД нагрева; п — нормаль к поверхности.
Слой гранул
Подложка
Рис. 1. Система координат при синтезе изделия
Важно отметить, что эффективный КПД нагрева — нелинейная функция, зависящая от температуры и материала. В ряде работ [9, 10] большое значение уделено расчету данного параметра. С увеличением температуры коэффициент поглощения чистой металлической поверхностью растет, а также существенно изменяется теплоемкость и теплопроводность металлов и сплавов, что необходимо учитывать в расчетах.
В работе [10] приведена следующая линейная зависимость от температуры для расчета погло-щательной способности и определения доли поглощенной мощности:
А(Т) = а + ЬТ, (3)
где а, Ь — коэффициенты.
К сожалению, зачастую аналитическое исследование ограничено рамками «линейности» задачи. Поставленная задача — нелинейная, требующая применения численных методов. Нелинейность обусловлена не только зависимостью теплофизических параметров материала от температуры, что само по себе значительно усложняет задачу, но и зависимостью источника тепла от времени (периодическое воздействие луча на формируемое изделие).
Поэтому для ее решения в работе используется метод конечных элементов (МКЭ), реализованный в среде программного комплекса А№У8 АРБЬ. Расчет температурного поля формируемого изделия проводился в трехмерной постановке с применением системы конечно-элементного анализа А№У8.
Построение модели расчета. Процедура проведения нестационарного теплового анализа состоит из четырех основных этапов: создания модели; задания граничных условий и тепловых нагрузок; решения задачи на основе уравнения теплопроводности; анализа результатов.
Процедура создания модели включает задание типа элементов, свойств материала и геометрии модели. Далее выполняется наложение сетки на модель (создание конечно-элементной модели) и задаются граничные и начальные условия. Поскольку в начальный момент времени температура модели совпадает с температурой окружающей среды, то последнюю и принимаем в качестве начальной температуры для всех узлов модели.
Нелинейные зависимости теплопроводности, плотности и теплоемкости титана приведены на рис. 2-4. Для температур до 1100 °С характерен
Вт/(м°С)
Т, 103 °С
Рис. 2. Зависимость коэффициента теплопроводности титана от температуры
р, 103-кг-м-3
О 0 1 02 0.3 0.4 0.5 0.6 07 0.9 0 9 1 1.1 1.2 1.3 1.4 1.5 16 1.7 1.8
Т, 103 °С
Рис. 3. Зависимость плотности титана от температуры
С, 103-Дж/(кг°С)
Т, 103 °С
Рис. 4. Зависимость теплоемкости титана от температуры
Рис. 5. Конечно-элементная модель
участок с малыми значениями коэффициента теплопроводности, что соответствует большей скорости нагрева. При более высоких температурах наступает резкое нелинейное увеличение коэффициента теплопроводности.
Значительные изменения теплопроводности Ц7), плотности р(7), теплоемкости с(Т) показывают, что система уравнений описывающая процесс электронно-лучевого синтеза, существенно нелинейная.
Конечно-элементная модель с нижней (ВОТ) и верхней (ТОР) поверхностями представлена на рис. 5. Электронный луч взаимодействует с верхней поверхностью ТОР.
Рис. 6. Описание свойств элемента SHELL131
Расчетная модель для определения температурного поля формируемого изделия создавалась с использованием многослойных оболочечных конечных элементов SHELL131 с 15 степенью свободы (TEMPBOT, TEMP2, ... TEMP14, TEMPTOP).
В качестве базового принимаем элемента SHELL131 — четырехугольный оболочечный элемент, который используется в трехмерных моделях. Описание свойств элемента SHELL131 приведено на рис. 6: изменение температуры в слое — по квадратичной зависимости, количество слоев равно 7.
Каждый из слоев включает три точки, причем одна точка одновременно общая для двух соседних слоев. Нумерация точек начинается от нижней поверхности BOT (первая точка) к верхней поверхности TOP (15 точка).
Принимаем, что с верхней и нижней поверхностей подложки (изделия) осуществляется радиационный теплоотвод. Плотность электронного излучения задавалась равномерным распределением по пятну нагрева радиусом r0 (r0 = d0/2). Поскольку с поверхности подложки тепло теряется радиационным излучением, то имеем задачу с нелинейным граничным условием.
Исходные данные модели. Режим обработки поверхности подложки электронным лучом определяется током луча, скоростью сканирования луча, диаметром пятна, геометрическими размерами изделия.
К исходным параметрам модели относятся: сторона подложки а = 70 мм; диаметр луча d0 = 1 мм; скорость сканирования Уск = 9060 мм/с; время операции ?опер = 35 c ток луча 1л = 30 мА; ускоряющее напряжение Uy = 60 кВ; смещение луча = 1 мм; время сканирования участка 1x1 мм при прямоугольном характере воздействия, с;
площадь перекрытия лучом в секунду Ап = VCi,d0 = 9060 мм2/с;
период ?пер = а2/(^^0^см) = 1,081678 с (Л^ — коэффициент перекрытия); частота f = 0,92449 Гц; число проходов N = 32; подведенное тепло QJ = 1л Цу?опер = 63000 Дж; плотность подведенного тепла к подложке или изделию за время процесса J^^ = Qj/а2 = = 12, 85 Дж/мм2.
В
^->
^ >
X >
* >
Л >
-------------------------------------------------------------------------------
*-------------------------------------------------------------------------------»
<-
С
см1
в
Рис. 7. Способ сканирования электронного луча по поверхности
А
Источник тепла перемещается по подложке со скоростью Способ сканирования электронного луча по поверхности показан на рис. 7.
Это эквивалентно приложению периодической тепловой нагрузки с длительностью импульса dпятна/VCк (при скорости сканирования луча равной 9060 мм/с и диаметре пятна dпятна = = 1 мм длительность импульса составляет 1,1-10-4 с).
Плотность теплового потока также является поверхностным граничным условием. Модель использует данную тепловую нагрузку, так как известна величина тепла, проходящего через поверхность (тепловой поток, приходящийся на единичную площадку):
HFLU = Q/Aэя = 1л иу/аэя2 = 1800/0,0022 =
= 0,45 109 Вт/м2, (4)
где ^эя — площадь элементарной ячейки; аэя — сторона квадратной элементарной ячейки.
Положительная величина плотности теплового потока означает, что тепло поглощается элементом. Стартовая позиция электронного луча показана темным цветом на рис. 8.
Электронный луч стартует с элемента № 1 и далее перемещается слева направо по узлам (нумерация элементов пропорционально увеличивается). В крайнем правом узле строки следует переход луча верх на один элемент (по столбцу) и перемещение луча по элементам справа налево по второй строке. Далее переход вверх на один элемент и повторение приведенных операций.
Таким образом, сканируется вся поверхность до последнего элемента.
Результаты расчета. Температурное поле на верхней поверхности во время второго прохода электронного луча над элементом, содержащим узел 668, который находится в центре исследуемой поверхности, показано на рис. 9.
Рассмотрим динамику изменения температуры в восьми верхних степенях свободы в узле 668 (точки узла №8, ... №15), которые наиболее сильно реагируют на воздействие электронного луча (рис. 10).
Из рис. 10 видно, что температура в узле на верхней поверхности (ТТОР_668) имеет зависимость от предистории (прохождения электронного луча рядом «до»), от настоящего (прохождение электронного луча «над») и практически не зависит от будущего (прохождение электронного луча «после»).
Рис. 8. Стартовая
позиция электронного луча
Рис. 9. Температурное поле на верхней поверхности во время второго прохода электронного луча над узлом 668
Ж9_668 ТЕ10_661
7Е1166 8 ТЕ12668
ЖЛ 36 68 7X14668
: . :■ ■ ; .Ми: : Г.- с
.667 1.001 1.335 1.669 2.003
Рис. 10. Динамика изменения температуры в 8-ми верхних степенях свободы (температуры) в узле 668 исследуемой модели
Рис. 11. Температурное поле верхней (а) и боковой (б) поверхности формируемого изделия
после двух циклов сканирования
Приращение температуры за один цикл сканирования электронным лучом поверхности размером 70x70 мм составляет примерно 20 °С, что соответствует экспериментальным данным.
Температурное поле верхней и боковой поверхности после двух циклов сканирования электронным лучом показано на рис. 11.
Подведем итоги.
Показана важность операции нагрева подложки для технологии послойного синтеза изделий из порошка различных сплавов, включая тугоплавкие.
Режим термической обработки поверхности электронным лучом должен обеспечивать условия для удержания гранул на поверхности подложки и формируемого изделия.
Для решения нелинейной модели, описывающей нагрев подложки электронным лучом, используется метод конечных элементов, который позволяет учитывать, во-первых, геометри-
ческую форму изделия, во-вторых, зависимость теплофизических параметров от температуры.
Проведенный расчет температурного поля позволяет качественно описывать процесс электронно-лучевого нагрева в вакуумной среде с учетом передачи тепла посредством излучения.
Модель адекватно отражает связь между действием электронного луча и распределением температурного поля.
Анализ результатов показывает, что для качественной оценки температурного поля можно использовать многослойные элементы 8ИБЬЬ131, что позволяет оперативно рассчитать температурное поле для многослойных конструкций.
Модель учитывает историю предыдущих периодических тепловых воздействий электронного луча. Нестационарный тепловой расчет позволяет определить температурное поле, изменяющееся во времени.
Работа выполнена при финансовой поддержке Министерства образования и науки РФ.
СПИСОК ЛИТЕРАТУРЫ
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. Murr L.E., Gaytan S.M., Ceylan A., Martinez E., Martinez J.L., Hernandez D.H., Machado B.I., Ramirez D.A., Medina F., Collins S., Wicker R.B. Characterization of titanium aluminide alloy components fabricated by additive manufacturing using electron beam meltingro // Acta Materialia. 2010. № 58. Р 1887-1894.
7. Biamino S., Penna A., Ackelid U., Sabbadini S., Tassa O., Fino P., Pavese M., Gennaro P., Badini C. Electron beam melting of Ti48Al2Cr2Nb alloy: Microstructure and mechanical properties investigation // Intermetallics. 2011. № 19. P. 776-781.
8. Самарский А.А., Вабищевич П.Н. Вычислительная теплопередача. М.: Едиториал УРСС, 2003. 784 с.
9. Веденов А.А., Гладуш Г.Г. Физические процессы при ЛО материалов. М.: Энергоатомиздат, 1985. 208 с.
10. Углов А.А., Смуров И.Ю., Лашин А.М., Гуськов
А.Г. Моделирование теплофизических процессов импульсного лазерного воздействия на металлы. М.: Наука, 1991. 288 с.
REFERENCES
1. Rudskoy A.I., Kondratyev S.Yu., 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.Yu., 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.Yu, Kokorin V.N. Pressovaniye geterofaznykh uvlazhnennykh poroshkovykh metallicheskikh smesey dlya povysheniya kachestva vyso-koplotnykh zagotovok s ispolzovaniyem metoda inten-sivnogo uplotneniya. Spravochnik. Inzhenernyy zhurnal s prilozheniyem. 2011. № 6. S. 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. 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. Murr L.E., Gaytan S.M., Ceylan A., Martinez E., Martinez J.L., Hernandez D.H., Machado B.I., Ramirez D.A., Medina F., Collins S., Wicker R.B. Characterization of titanium aluminide alloy components fabricated by additive manufacturing using electron beam meltingyu. Acta Materialia. 2010. № 58. S. 1887-1894. (rus.)
7. Biamino S., Penna A., Ackelid U., Sabbadini S., Tassa O., Fino P., Pavese M., Gennaro P., Badini C. Electron beam melting of Ti48Al2Cr2Nb alloy: Microstructure and mechanical properties investigation. Intermetallics. 2011. № 19. P. 776-781. (rus.)
8. Samarskiy A.A., Vabishchevich P.N.Vychislitelnaya teploperedacha. M.: Yeditorial URSS, 2003. 784 s. (rus.)
9. Vedenov A.A., Gladush G.G. Fizicheskiye protsessy pri LO materialov. M.: Energoatomizdat, 1985. 208 s. (rus.)
10. Uglov A.A., Smurov I.Yu., 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: [email protected]
СОКОЛОВ Юрий Алексеевич — кандидат технических наук главный конструктор ООО «НТК Электромеханика»; 172386, ул. Заводское шоссе, 2, г. Ржев, Россия. E-mail: [email protected] КОПАЕВ Владимир Николаевич — ведущий инженер, ООО «НТК Электромеханика». 172386, ул. Заводское шоссе, 2, г. Ржев, Россия. E-mail: [email protected]
AUTHORS
RUDSKOI Andrey I. — 29, Politechnicheskaya St., St. Petersburg, 195251, Russia.E-mail: [email protected] SOKOLOV Yuriy A. — 2, Zavodskoe shosse, Rzhev city Tver region, Russia, 172386. E-mail: [email protected] KOPAEV Vladimir N. — 2, Zavodskoe shosse, Rzhev city Tver region, Russia, 172386. E-mail: vkopaev_ela@ mail.ru
© Санкт-Петербургский государственный политехнический университет, 2014
