Брянск, Брянская область, Россия
Рассматривается создание математической модели процесса синтеза изделий из металлической проволоки методом FDM. Вычисляется энергия, необходимая для фазового перехода от твердой проволоки к жидкому расплаву материала. С некоторыми допущениями используется уравнение Навье-Стокса для несжимаемой жидкости для вычисления средней скорости течения расплавленного материала, что позволяет определить общее время, требуемое на формирование капли материала. Вычисляется скорость перемещения сопла подающей головки и скорость подачи проволоки в рабочую зону для получения дорожки материала близкой к идеальной. Расчет изменения скорости перемещения сопла подающей головки и скорость подачи проволоки в рабочую зону для нивелирования дефектов печати, таких как провалы и наплывы материала.
моделирование, управление, математическая модель, плавление металла, скорость подачи проволоки, скорость перемещения сопла, синтеза деталей методом FDM
Введение
В настоящее время существуют несколько методов позволяющих контролировать точность синтезирования деталей, такие как ситуационное управление качеством изготавливаемых изделий [1], нейросетевое управления процессом электродугового послойного формирования [2].
Использование методов, позволяющих контролировать качество синтезируемых деталей часто требует применения большого количества дополнительных устройств и сложных алгоритмов обработки информации и принятия решений.
Одним из альтернативных методов управления качеством синтезируемых изделий является контроль за высотой наплавляемого слоя материала и соответствующая корректировка управляющих воздействия для компенсации обнаруженных дефектов. Такая система обнаруживает разницу в высоте слоя между CAD-файлом и синтезируемым изделием. При обнаружении расхождений корректируется скорость подачи материала или скорость перемещения подающей головки, что приводит к изменению объема наплавленного в единицу времени материала.
Разработка математической модели процесса синтеза металлических изделий из проволоки
Процесс синтеза деталей методом FDM заключается в подаче металлической проволоки в рабочую зону, где проволока расплавляется лазерным лучом (рис. 1, а), в результате чего образуется капля расплавленного материала (рис. 1, б), подающая головка смещается и цикл повторяется.
а) б)
Рис. 1. Схематичное изображение процесса:
а – синтеза детали из металлической проволоки; б – формирования капли металла.
Fig. 1. Schematic representation of the process:
a – synthesis of a part from metal wire; b – formation of a drop of metal.
В первом приближении, энергию, которую необходимо сообщить телу для изменения его агрегатного состояния из твердого в жидкое можно разделить на следующие составляющие [3, 4]:
– энергия для нагрева твердого тела до температуры плавления, определяется как произведение удельной теплоемкости c, 
на массу нагреваемого тела m, кг на величину изменения температуры тела 
, °C:
|
|
(1) |
где 
– изменение температуры тела, °C.
– энергия, необходимая для фазового перехода от твердого к жидкому состоянию (теплота плавления):
|
|
(2) |
где 
– удельная теплота плавления, 
.
В результате получаем следующую формулу:
|
|
(3) |
Определяем массу тела, которая будет расплавлена в единицу времени:
|
|
(4) |
где 
– плотность тела, 
; d – диаметр, см; 
– высота, см.
– энергия, необходимая для плавления, поступает с некоторой мощностью P. Время плавления (
) может быть выражено как отношение энергии к мощности источника тепла:
|
|
(5) |
где 
– мощность теплового источника, нагревающего тело, Вт.
Для дальнейших вычислений рассмотрим плоскость, наклоненную к горизонту под углом θ, °. По плоскости под действием силы тяжести течет жидкость с толщиной слоя h. Для нахождения скорости течения жидкости воспользуемся уравнением Навье-Стокса для несжимаемой жидкости:
|
|
(6) |
где 
– вектор скорости, м/с; 
– градиент давления, Па/м; 
– вязкость жидкости, Па*с; 
– внешние силы (сила тяжести).
Принимаем следующие допущения:
Слой жидкости тонкий, поэтому градиент давления отсутствует 
), с учетом симметрии и несжимаемости жидкости 
.
Получаем уравнение:
|
|
(7) |
где 
– ускорение свободного падения, 
.
Откуда выводим уравнение средней скорости течения жидкости:
|
|
(8) |
Таким образом, приблизительное время стекания расплавленного материала на расстояние L, см под углом 90° к горизонту можно вычислить по формуле:
|
|
(9) |
Тогда общее время, требуемое на формирование капли:
|
|
(10) |
Стоит отметить, что данная математическая модель не учитывает свойства подложки, на которую будет наноситься расплав, что в свою очередь делает невозможным учет эффекта смачивания. А также для упрощения вычислений было принято, что телу сообщается энергии не больше, чем требуется для поддержания его в состоянии расплава [5, 6].
Зная время, за которое проволока расплавится до состояния формирования дорожки и зная длину проволоки L, которая может быть расплавлена в единицу времени можно найти скорость, с которой необходимо подавать проволоку для ее непрерывного плавления:
|
|
(11) |
Необходимо найти скорость перемещения сопла, при которой синтезируемая дорожка будет приближена к идеальной.
Для этого необходимо разделить процесс синтеза дорожки на равные временные участки [7 – 9], где каждый временной участок равен 
. По прошествии времени будет сформирована капля расплавленного материала длины l, как указано на рис. 2.
Рис. 2. Идеальная дорожка материала:
1 – капля материала; 2 – недостающий объем материала
Fig. 2. Ideal material path:
1 – drop of material; 2 – missing volume of material
Объем, получившийся капли (без учета потери объема при испарении металла):
|
|
(12) |
Радиус капли (принимаем допущение, что капля принимает форму полусферы):
|
|
(13) |
Недостающий объем материала (на рис. 2 заштрихованная область):
|
|
(14) |
Недостающий объем 
будет заполняться по мере уменьшения расстояния между центрами капель расплава. Из интегрального уравнения объема тела вращения необходимо найти пределы ограничивающие кривую, объем которой, при вращении ее вокруг оси абсцисс даст объем равный объему 
.
|
|
(15) |
где 
– радиус тела вращения, см; 
– нижний предел интегрирования, см; 
– верхний предел интегрирования, см.
|
|
(16) |
|
|
(17) |
|
|
(18) |
Результат решения интегрального уравнения:
|
|
(19) |
|
|
(20) |
Верхний предел интегрирования b известен и равен 
– необходимо найти нижний предел интегрирования a:
|
|
(21) |
После нахождения верхнего предела интегрирования скорость перемещения головки:
|
|
(22) |
Как видно из рис. 3 из-за того, что капли имеют сферическую форму, их сближении на рассчитанное расстояние приведет к выдавливанию недостаточного для заполнения пустоты объема материала. Следовательно, необходимо ввести поправочный коэффициент, обеспечивающий приемлемую скорость синтеза изделия и их качество (рис. 4).
Рис. 3. Визуализация результата вычислений:
1 – капля материала; 2 – расчетный объем материала для выдавливания; 3 – выдавливаемый объем материала
Fig. 3. Visualization of calculation results:
1 – is a drop of material; 2 – is the calculated volume of the material to be extruded; 3 – is the volume of the material to be extruded
Рис. 4. Визуализация результата вычислений с поправкой:
1 – капля материала; 2 – расчетный объем материала для выдавливания; 3 – выдавливаемый объем материала с учетом поправки
Fig. 4. Visualization of the calculation result with correction:
1 – is a drop of material; 2 – is the calculated volume of the material to be extruded; 3 – the volume of material being squeezed out, taking into account the correction
Используя полученные формулы, можно произвести расчет поправочного коэффициента скорости перемещения сопла или подачи проволоки, для компенсации дефектов синтеза изделий. Для этого в равные интервалы времени измеряем высоту наплавленного слоя материала и при определении расхождении между заданной высотой слоя и измеренными данными производится корректировка скорости подачи проволоки или скорости перемещения сопла, что позволит исправлять дефекты нивелируя их влияние на последующие слои [10].
Опираясь на разработанную математическую модель, было вычислено изменение скорости перемещения сопла и подачи проволоки для корректировки провала и выпуклого дефекта для процесса синтеза изделия со следующими параметрами: материал алюминий, диаметр проволоки 3 мм.
а) б)
Рис. 5. Графики изменения скорости подачи проволоки и перемещения сопла:
а – для корректировки провала; б – для корректировки выпуклого дефекта
Fig. 5. Graphs of changes in wire feed speed and nozzle movement:
a – to correct a dip; b – to correct a convex defect
Таблица 1
Изменение скорости подачи проволоки и перемещения сопла для корректировки провала
Table 1
Changing wire feed speed and nozzle movement to correct dip
|
1 |
2 |
3 |
4 |
5 |
6 |
7 |
8 |
9 |
10 |
|
|
Vп, мм/с |
3,50 |
3,50 |
3,50 |
3,65 |
4,53 |
3,65 |
3,50 |
3,50 |
3,65 |
4,53 |
|
Vг, мм/с |
1,68 |
1,68 |
1,68 |
1,63 |
1,34 |
1,63 |
1,68 |
1,68 |
1,63 |
1,34 |
|
h, мм |
2,16 |
2,16 |
2,16 |
1,91 |
0 |
1,91 |
2,16 |
2,16 |
1,91 |
0 |
Таблица 2
Изменение скорости подачи проволоки и перемещения сопла для корректировки выпуклого дефекта
Table 2
Changing wire feed speed and nozzle movement to correct a convex defect
|
1 |
2 |
3 |
4 |
5 |
6 |
7 |
8 |
9 |
10 |
|
|
Vп, мм/с |
3,50 |
3,50 |
3,50 |
3,08 |
2,34 |
3,08 |
3,50 |
3,50 |
3,65 |
4,53 |
|
Vг, мм/с |
1,68 |
1,68 |
1,68 |
1,81 |
2,05 |
1,81 |
1,68 |
1,68 |
1,63 |
1,34 |
|
h, мм |
2,16 |
2,16 |
2,16 |
2,76 |
3,09 |
2,76 |
2,16 |
2,16 |
1,91 |
0 |
Заключение
Принимая описанные в статье допущения, разработанная математическая модель позволяет моделировать процесс синтеза изделий из металической проволоки методом FDM. Результаты приведенные на рис. 5, табл. 1 и табл. 2 показывают, что математическая модель может быть использована для вычисления скорости подачи проволоки и скорости перемещения сопла, изменяя которые можно нивилировать различные дефекты печати и добиться повышения качества синтезируемых деталей.
1. Колесников А.А. Ситуационное управление обеспечением качества изготовления изделий по аддитивной технологии// Образовательные ресурсы и технологии. – 2017. – Т. 21, № 4. – С. 84-92.
2. Паршин С.Г., Майстро А.С., Переверзев А.Е. Цифровая технология и оборудование для аддитивного производства объемных изделий из цветных сплавов на основе нейросетевого управления процессом электродугового послойного формирования // IV Международная конференция «Аддитивные технологии: настоящее и будущее. – 2018. – С. 285-298.
3. Гаврилин И.В., Плавление и кристаллизация металлов и сплавов / Гаврилин И.В. – Владимир: Владим. гос. ун-т, 2000 – 260 с.
4. Александров Д.В., Зубарев А.Ю., Искакова Л.Ю. Введение в гидродинамику: [учеб. пособие]. – Екатеринбург: Изд-во университета Урал, 2012 – 112 с.
5. Краснова Е.В., Саушкин Б.П. Аддитивное формообразование изделий из металлов и сплавов пучком электронов. Селективное плавление (часть 1) // Аддитивные технологии. – 2021. – № 1. – С. 30-41.
6. Краснова Е.В., Саушкин Б.П. Аддитивное формообразование изделий из металлов и сплавов пучком электронов. прямая подача энергии и материала в зону плавления (часть 2) // Аддитивные технологии. – 2021. – № 2. – С. 30-43.
7. Краснова Е.В., Саушкин Б.П. Гибридные технологии и оборудование // «Аддитивные технологии». – 2021. – № 4. – С. 48-58.
8. Евтодьева М.Г. Аддитивное производство и дополненная реальность как новые производственные технологии в авиационной отрасли // Вестник МГИМО. – 2020. – №13 (5). – С. 307-330.
9. Система планирования и подготовки аддитивного производства / В.Г. Засканов [и др.] // Известия Самарского научного центра Российской академии наук. – 2018. – Т. 20. – № 6. – С.14-23.
10. Аналитическое исследование. Рынок технологий 3D-печати в России и мире. Группа «Деловой профиль» | MGI Worldwide. 2020. 15 с.
