Аналитическое доказательство законов подобия для технологии аддитивного лазерного выращивания

Обложка

Цитировать

Полный текст

Аннотация

Показано, что задачу описания технологии аддитивного лазерного выращивания можно рассматривать в рамках автомодельного уравнения теплопроводности. Показано, что при некоторых ограничениях глубина проплава подложки хорошо описывается автомодельным решением. На основе полученного автомодельного решения получены двухпараметрическая зависимость глубины проплава от числа Пекле (отношения скорости сканирования к скорости изменения температуры в материале) и безразмерной энтальпии (отношение удельной энергии, поглощенной материалом, и энергии, необходимой для плавления). Показано, что полученная аналитическая зависимость хорошо описывает экспериментальные данные.

Полный текст

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

Для уменьшения количества натурных и численных экспериментов исследуются законы подобия (или, как их еще называют, законы масштабирования) процессов лазерной наплавки [3–6]. В этих работах показано, что параметры зоны проплавления зависят от двух безразмерных параметров – нормализованной энтальпии B и числа Пекле Pe (отношения скорости сканирования к скорости изменения температуры в материале). В работе [4] на основе простой тепловой модели [7] были получены интерполяционные формулы в виде степенных функций от двух параметров (числа Пекле и нормализованной энтальпии), описывающих форму “ванны” расплава (ширину и глубину). В работе [5] показано, что независимо от объемной доли керамики безразмерные геометрические параметры единичного трека (глубина, ширина и высота) зависят от числа Пекле и нормализованной энтальпии. Эти зависимости были аппроксимированы алгебраическими выражениями.

УРАВНЕНИЯ ТЕПЛОПРОВОДНОСТИ ДЛЯ ЛАЗЕРНОГО НАГРЕВА

Рассмотрим тепловой процесс образования “ванны” расплава, не учитывая слой порошка, движение расплава и испарение материала. Некоторые обоснования такого упрощения задачи лазерной наплавки приведены в работе [4]: 1) теплопроводность порошка намного меньше, чем подложки; 2) порошок имеет более высокую абсорбционную способность, чем расплавленная поверхность; 3) параметры лазерной обработки будут в основном определяться на поздней стадии процесса, а влияние порошка в основном сказывается на начальном этапе. Не будет учитываться теплота плавления подложки, поскольку для металлов она мала по сравнению с энергией, необходимой для нагрева материала до температуры плавления.

Рассмотрим лазерный луч, сканирующий подложку со скоростью V вдоль прямой линии по поверхности. Поверхность при этом нагревается. Рассмотрим одномерную задачу теплопроводности вглубь поверхности (координата r). Одномерное уравнение теплопроводности в эйлеровых координатах для случаев плоской ν=0, осевой ν=1 и сферической ν=2 симметрии имеет вид

cTt=1ρrνrrνW,   W=λTr, (1)

где c – удельная теплоемкость, W – плотность теплового потока,  ρ– постоянная плотность среды,  λ– коэффициент теплопроводности,  T– температура.

Перейдем к массовым переменным Лагранжа:

m=ρν+1rν+1. (2)

Тогда одномерное уравнение теплопроводности в массовых переменных Лагранжа координатах запишется в виде

cTt=W~m,    W~=2λmTm. (3)

В процессе лазерной наT~m,t=Tm,tT0плавки луч движется по прямой вдоль поверхности. Число Пекле процесса обычно больше единицы – луч движется быстрее, чем изменяется температура в подложке, поэтому для описания процесса теплообмена в подложке лучше всего подходит одномерное приближение с осевой симметрией ν=1. В этом случае уравнения (2) и (3) станут следующими:

m=ρ2r2,    cTt=W~m,    W~=2λmTm. (4)

ПОСТАНОВКА АВТОМОДЕЛЬНОЙ ЗАДАЧИ

Будем рассматривать процесс теплообмена при лазерной наплавке в автомодельной постановке. Пусть в момент времени t=0 вдоль оси симметрии m=0 r=0 выделяется энергия E0=const. В последующие моменты времени t>0 тепло распространяется от оси симметрии. В начальный момент времени t=0 среда является холодной:

Tm,0=0. (5)

Если начальная температура Tm,0=T0>0, то замена переменной T~m,t=Tm,tT0 возвращает граничное условие к условию T~m,0=0 с холодной средой.

В данной постановке автомодельная задача для уравнений (4) с условием (5) была решена в работе [8]. Автомодельное решение [8] имеет вид

s=cm2λt,    fs=2λtE0Tm,t,    ωs=W~tE0. (6)

Переходя к размерным переменным Эйлера из решения (6) и первого уравнения системы (4) найдем распределение температуры в среде:

Tr,t=E04πλ1texp14αr2t, (7)

где α = λ/сρ – коэффициент температуропроводности.

СРАВНЕНИЕ С ЭКСПЕРИМЕНТАМИ

Рассмотрим эксперименты по прямому лазерному выращиванию металлокерамических треков [5].

Для описания термодинамических процессов введем число Пекле [9]:

Pe=VD/α, (8)

где V– скорость сканирования, D – диаметр лазерного пучка. Число Пекле представляет собой отношение скорости сканирования к скорости изменения температуры в материале и не учитывает процесс поглощения лазерного излучения и плавление материала.

Введем безразмерный энергетический параметр процесса лазерной наплавки – безразмерную энтальпию B. Она представляет собой отношение удельной энергии, поглощенной материалом, к энергии, необходимой для плавления:

B=ΔHhs=AWhsαπVD3=ATmπcρλWVD3, (9)

где  ΔH– изменение энтальпии, hs=ρcTm – энтальпия плавления [3],  Tm– температура плавления,  A– интегральный коэффициент поглощения,  W– мощность лазерного излучения.

Время действия лазерного излучения на подложку зависит от скорости сканирования V и диаметра лазерного пучка D с некоторым коэффициентом n: t*=nD/V. За это время глубина проплава достигнет значения r*. На границе проплава температура станет равной температуре плавления: Tr*,t*=Tm. Энергия, которая выделяется при сканировании подложки лазером, равна E0=AW/V. Подставляя значения времени t* и энергии E0 в уравнение (7) и используя безразмерные параметры (8) и (9), получим формулу для радиуса проплава:

r*D=k2nPelnPeB216πn2, (10)

где  n– параметр,  k=w/2h– коэффициент формы проплава [5],  w– ширина проплава,  h– глубина проплава.

Формула предполагает, что “ванна” проплава имеет полукруглую форму, т.е. коэффициент k=1. Это связано с одномерным решением с осевой симметрией. В реальности “ванна” проплава не всегда имеет полукруглую форму. На рис. 1 приведены схема и фотографии поперечных сечений единичных треков, показывающие переход от режима теплопроводности к режиму “кинжального” проплавления с характерными значениями безразмерной энтальпии B и безразмерного числа Пекле Pe. Поэтому для расчета некруглой формы проплава в формулу введен коэффициент k.

 

Рис. 1. Схема (а) и фотографии поперечных сечений единичных треков при лазерной наплавке порошковой смеси ВТ-6 + 10% масс. SiC: режим теплопроводности (B=1.46 ,Pe=17.51 ) (б); режим “кинжального” проплавления ( B=14.21, Pe=2.93) (в).

 

Экспериментальные данные работы [5] обсчитаны по формуле (10). Результаты сравнения экспериментов и теоретической формулы приведены на рис. 2 и 3.

 

Рис. 2. Безразмерная полуширина проплава w/2D в зависимости от безразмерной энтальпии B. Сравнение экспериментальных результатов (кружки) с расчетами (линия) по формуле (10) при n = 0.35, k = 1.

 

Рис. 3. Безразмерная полуширина проплава w/2D в зависимости от числа Пекле Pe. Сравнение экспериментальных результатов (кружки) с расчетами (линия) по формуле (10) при n = 0.35, k = 1.

 

Из рис. 2, 3 видно, что формула (10) хорошо описывает экспериментальные данные, за исключением больших чисел безразмерной энтальпии. Это может объясняться тем, что происходит переход от режима теплопроводности к режиму “кинжального” проплавления, для которого характерно наличие парогазового канала и глубокое проникновение излучения внутрь подложки с существенным ее проплавлением. Как результат, происходит изменение коэффициента поглощения лазерного излучения материалом A, который учитывается в расчетах безразмерной энтальпии.

На рис. 4 показано сравнение экспериментальных [5] и расчетных по формуле данных по глубине проплавления. Видно, что формула (10) хорошо описывает эксперименты с зоной теплопроводности, но не описывает зону “кинжального” проплавления.

 

Рис. 4. Безразмерная глубина проплава h/D в зависимости от числа Пекле Pe. Сравнение экспериментальных результатов (кружки) с расчетами (линия) по формуле (10) при n = 0.35, k = 1. Выделены области: I – зона “кинжального” проплавления, II – зона теплопроводности.

 

ЗАКЛЮЧЕНИЕ

Задача описания технологии аддитивного лазерного выращивания рассмотрена в рамках автомодельного уравнения теплопроводности. Показано, что решение автомодельной задачи для уравнения теплопроводности с осевой симметрией хорошо описывает параметры “ванны” расплава. Получена двухпараметрическая зависимость ширины и глубины проплава от числа Пекле и безразмерной энтальпии. Проведено сравнение с экспериментальными данными [5]. Показано, что полученная аналитическая зависимость хорошо описывает экспериментальные данные c режимом теплопроводности, но не описывает “кинжальный” режим проплавления.

ИСТОЧНИК ФИНАНСИРОВАНИЯ

Работа выполнена в рамках государственного задания ИТПМ СО РАН.

×

Об авторах

В. М. Фомин

Институт теоретической и прикладной механики им. С.А. Христиановича Сибирского отделения Российской академии наук

Автор, ответственный за переписку.
Email: fomin@itam.nsc.ru

Академик РАН

Россия, Новосибирск

А. А. Голышев

Институт теоретической и прикладной механики им. С.А. Христиановича Сибирского отделения Российской академии наук

Email: alexgol@itam.nsc.ru
Россия, Новосибирск

А. Е. Медведев

Институт теоретической и прикладной механики им. С.А. Христиановича Сибирского отделения Российской академии наук

Email: medvedev@itam.nsc.ru
Россия, Новосибирск

А. Г. Маликов

Институт теоретической и прикладной механики им. С.А. Христиановича Сибирского отделения Российской академии наук

Email: smalik707@yandex.ru
Россия, Новосибирск

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

  1. Mukherjee T., DebRoy T. Control of asymmetric track geometry in printed parts of stainless steels, nickel, titanium and aluminum alloys // Computational Materials Science. 2020. V. 182. 109791. https://doi.org/10.1016/J.COMMATSCI.2020.109791
  2. Громов В.Е., Иванов Ю.Ф., Ефимов М.О., Шлярова Ю.А. Структура и свойства высокоэнтропийного сплава ALCRFECONI после электронно-ионно-плазменной обработки // Доклады РАН. Физика, технические науки. 2023. T. 511. № 1. С. 5–9. https://doi.org/10.31857/S2686740023040041
  3. Weaver J.S., Heigel J.C., Lane B.M. Laser spot size and scaling laws for laser beam additive manufacturing // J. Manufacturing Processes. 2022. V. 73. № August 2021. P. 26–39. https://doi.org/10.1016/j.jmapro.2021.10.053
  4. Rubenchik A.M., King W.E., Wu S.S. Scaling laws for the additive manufacturing // J. Materials Processing Technology. 2018. V. 257. P. 234–243. https://doi.org/10.1016/j.jmatprotec.2018.02.034
  5. Golyshev A.A., Malikov A.G. Scaling laws for the additive manufacturing of the AISI 316 L deposited by laser surface cladding and direct metal deposition methods // Optik. 2023. V. 295. August. № 171506. https://doi.org/10.1016/j.ijleo.2023.171506
  6. Голышев А.А., Сибирякова Н.А. Законы подобия при прямом лазерном вырашивании металлокерамических треков // Прикладная механика и техническая физика. 2023. V. 64. № 5. P. 102–107. https://doi.org/10.15372/PMTF202315287
  7. Eagar T.W., Tsai N.S. Temperature Fields Produced By Traveling Distributed Heat Sources // Welding Journal (Miami, Fla). 1983. V. 62. № 12. P. 346–355.
  8. Волосевич П.П., Леванов Е.И. Автомодельные решения задач газовой динамики и теплопереноса. М.: Изд-во МФТИ, 1997. 240 с.
  9. Mukherjee T., Manvatkar V., De A., DebRoy T. Dimensionless numbers in additive manufacturing // Journal of Applied Physics. 2017. V. 121. № 064904. https://doi.org/10.1063/1.4976006

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

Доп. файлы
Действие
1. JATS XML
2. Рис. 1. Схема (а) и фотографии поперечных сечений единичных треков при лазерной наплавке порошковой смеси ВТ-6 + 10% масс. SiC: режим теплопроводности ( , ) (б); режим “кинжального” проплавления ( , ) (в).

Скачать (164KB)
3. Рис. 2. Безразмерная полуширина проплава w/2D в зависимости от безразмерной энтальпии B. Сравнение экспериментальных результатов (кружки) с расчетами (линия) по формуле (10) при n = 0.35, k = 1.

Скачать (84KB)
4. Рис. 3. Безразмерная полуширина проплава w/2D в зависимости от числа Пекле Pe. Сравнение экспериментальных результатов (кружки) с расчетами (линия) по формуле (10) при n = 0.35, k = 1.

Скачать (94KB)
5. Рис. 4. Безразмерная глубина проплава h/D в зависимости от числа Пекле Pe. Сравнение экспериментальных результатов (кружки) с расчетами (линия) по формуле (10) при n = 0.35, k = 1. Выделены области: I – зона “кинжального” проплавления, II – зона теплопроводности.

Скачать (99KB)

© Российская академия наук, 2024

Согласие на обработку персональных данных с помощью сервиса «Яндекс.Метрика»

1. Я (далее – «Пользователь» или «Субъект персональных данных»), осуществляя использование сайта https://journals.rcsi.science/ (далее – «Сайт»), подтверждая свою полную дееспособность даю согласие на обработку персональных данных с использованием средств автоматизации Оператору - федеральному государственному бюджетному учреждению «Российский центр научной информации» (РЦНИ), далее – «Оператор», расположенному по адресу: 119991, г. Москва, Ленинский просп., д.32А, со следующими условиями.

2. Категории обрабатываемых данных: файлы «cookies» (куки-файлы). Файлы «cookie» – это небольшой текстовый файл, который веб-сервер может хранить в браузере Пользователя. Данные файлы веб-сервер загружает на устройство Пользователя при посещении им Сайта. При каждом следующем посещении Пользователем Сайта «cookie» файлы отправляются на Сайт Оператора. Данные файлы позволяют Сайту распознавать устройство Пользователя. Содержимое такого файла может как относиться, так и не относиться к персональным данным, в зависимости от того, содержит ли такой файл персональные данные или содержит обезличенные технические данные.

3. Цель обработки персональных данных: анализ пользовательской активности с помощью сервиса «Яндекс.Метрика».

4. Категории субъектов персональных данных: все Пользователи Сайта, которые дали согласие на обработку файлов «cookie».

5. Способы обработки: сбор, запись, систематизация, накопление, хранение, уточнение (обновление, изменение), извлечение, использование, передача (доступ, предоставление), блокирование, удаление, уничтожение персональных данных.

6. Срок обработки и хранения: до получения от Субъекта персональных данных требования о прекращении обработки/отзыва согласия.

7. Способ отзыва: заявление об отзыве в письменном виде путём его направления на адрес электронной почты Оператора: info@rcsi.science или путем письменного обращения по юридическому адресу: 119991, г. Москва, Ленинский просп., д.32А

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

9. Порядок уничтожения персональных данных при достижении цели их обработки или при наступлении иных законных оснований определяется Оператором в соответствии с законодательством Российской Федерации.

10. Я согласен/согласна квалифицировать в качестве своей простой электронной подписи под настоящим Согласием и под Политикой обработки персональных данных выполнение мною следующего действия на сайте: https://journals.rcsi.science/ нажатие мною на интерфейсе с текстом: «Сайт использует сервис «Яндекс.Метрика» (который использует файлы «cookie») на элемент с текстом «Принять и продолжить».