Application of the normalized diagram to simulate additive surfacing using energy beam sources
- Authors: Letyagin I.Y.1, Trushnikov D.N1, Belenkiy V.Y.1,2, Erikov A.P2, Olshanskaya T.V1
- Affiliations:
- Perm National Research Polytechnic University
- Perm Military Institute of National Guard Troops of the Russian Federation
- Issue: Vol 22, No 3 (2020)
- Pages: 68-75
- Section: ARTICLES
- URL: https://ered.pstu.ru/index.php/mm/article/view/2955
- DOI: https://doi.org/10.15593/2224-9877/2020.3.09
- Cite item
Abstract
The actively developing additive manufacturing makes it possible to use such beam energy sources as an electron beam and a laser beam for surfacing metal parts from powder materials. The quality of the resulting structure depends on many factors. These factors include the thermo physical properties of the powder and substrate, as well as the physical processes of the interaction of energy beams and materials. Methods for estimating the parameters of additive technologies described in the literature often make it difficult to verify and predict the results of calculations. The proposed method for the use of dimensionless parameters in a simple analytical model for heating samples coated with a layer of powder makes it possible to reduce the number of variables. The thermal model of the process provides a connection between the dimensionless coordinates and the values of the dimensionless temperature. A normalized diagram was built for the processing of a powder material with additive technologies using concentrated energy beams. The diagram is constructed on the basis of a simple analytical model of heating powder layers by a moving linear heat source. Recommended for selective melting is the interval of the diagram for points subject to processing below the line of maximum dimensionless temperature or below the line of maximum thermal efficiency. The range of processing parameters, as well as the ability to predict the optimal energy parameters of an electron beam or a laser beam, with selective melting can be determined depending on the dimensionless beam velocity. The diagram represents useful background information and methodology to select the appropriate processing parameters in the early stages of development not only the sintering additive powder materials, but also to assess the depth of heating of a particular material. It can be used in the heat treatment using the energy beam sources.
Full Text
Введение Аддитивное производство (АП) представляет собой процесс создания трехмерных объектов путем послойного формирования. АП - это быстро развивающийся производственный процесс, имеющий значительный потенциал для снижения стоимости производства, энергопотребления и воздействия на окружающую среду. В качестве концентрированного энергетического пучка в аддитивном производстве может использоваться электронный пучок, лазерный луч или плазменный поток. Несмотря на преимущества в конструкции и стоимости, АП с использованием металлических материалов находится на начальной стадии развития. Это связано с ограниченными знаниями о совокупности физических процессов при АП и более сложных, чем при традиционных производственных процессах. К влияющим на эти процессы факторам относятся: взаимодействие энергетических лучей и материалов, теплообмен и поток расплавленного металла, фазовые превращения, термические напряжения и деформации. Все эти факторы влияют на качество полученного изделия и свойства материала [1-3] полученной конструкции. Комплексный обзор аддитивного производства из металлических изделий приведен в работе [4]. Многие существующие и легко реализуемые программные пакеты могут помочь понять некоторые функции процесса [5-8]. В работе [5] применяются методы теплопроводности и конечных элементов. В источнике [6] компьютерное моделирование основано на алгоритме теплопереноса и течении жидкости. В работе [7] используется метод набора уровней. В источнике [8] объем жидкости моделируется с использованием метода конечных разностей. Применение метода решеточных уравнений Больцмана и произвольного Лагранж - Эйлерова поведения представлено в работе [9]. Наличие различных допущений, таких как отсутствие учета некоторых физических явлений в моделировании и отсутствие данных о некоторых элементах процессов (градиент температур на поверхности жидкой ванны, роль примесей сплава в поверхностном натяжении, вязкость и различная неустойчивость в потоках расплавленной ванны и т.д.), часто затрудняет проверку и прогнозирование результатов вычислений. В этой статье представлена простая аналитическая модель нагрева образцов, покрытых слоем порошка, а также используются известные безразмерные параметры, которые могут уменьшить общее количество рассматриваемых переменных. Таким образом, важно понимать, что зависимости между отдельными переменными процесса не могут быть найдены. Во время АП расплавленная ванна мала и градиенты температур в твердом образце (слой порошка и основной материал) обусловлены размерами расплавленной ванны. Эти градиенты температур объясняют распространение остальной поглощенной энергии, которая не используется для плавления, за счет влияния теплопроводности. В случае селективного плавления последовательно расположенных порошковых слоев при получении конструкций подложка также нагревается. Температура селективно расплавленной области определяется увеличением температуры подложки до температуры плавления обрабатываемого порошка. Термическая модель Решение уравнения теплового баланса при нагреве порошкового слоя толщиной L линейным движущимся тепловым источником с постоянной энергией Р, скоростью V (рис. 1), при условии отсутствия фазовых изменений в образце при теплопередаче и при известных физических параметрах материала: теплопроводность λ, температуропроводность a (a = λ/(Cρ), где C - удельная теплоемкость, ρ - плотность слоя порошка), примет вид [20] (1) где T - температура исследуемой точки; T0 - температура подложки (начальная температура); x - расстояние от движущегося центра; r - радиус-вектор, движущийся вместе с источником тепла; K0 - модифицированная функция Бесселя второго рода нулевого порядка. Рис. 1. Определение параметров геометрии при селективном плавлении первого порошкового слоя: 2rb - диаметр пятна нагрева; L - толщина слоя порошка Заметим, что P - поглощенная энергия луча (энергия пучка Pb с учетом потери энергии за счет обратного рассеяния и вторичных электронов в случае электронного пучка, а в случае лазерного луча - поверхностного поглощения энергии, или коэффициента связи). Чтобы свести к минимуму влияние температурной зависимости тепловых констант, значения λ, C и a взяты для промежуточной температуры (между начальными температурами T0 и температуры плавления Tm), а процесс теплопроводности считается не зависящим от температуры. Из равенства к нулю первой производной теплового цикла для точки выборки видно, что все точки из движущейся системы координат с максимальными температурами задаются (2) или из условий геометрии: (2а) где - модифицированная функция Бесселя второго рода первого порядка и - число Пекле (безразмерное) расстояний до источника тепла, которое имеет значения, пропорциональные скорости движения. Тогда из формул (1) и (2) можно найти безразмерную максимальную температуру qm: (3) Функция безразмерных параметров qm (rV/2a) является более практичной и может быть найдена с помощью итераций из выражения (3) с учетом выражения (2а), при этом учитывается, что y = = b/2 = rsinφ, где b - ширина расплавленного и затвердевшего слоя. Эта функция приведена на рис. 2 для диапазона Vy/2a, который подходит для нагрева луча на малых расстояниях от движущегося линейного источника нагрева [21]. Другой вид Y(X) этого соотношения приведен также на рис. 2 с использованием наклонных безразмерных координат X = P / Lλ(Tm - T0) и Y = Vb / 2a. Рис. 2. Зависимость максимальной безразмерной температуры qm от безразмерных расстояний Vy/2a с безразмерными координатами Y = Y (X) (наклонены на 90°) Термический КПД ηt процесса нагрева порошкового слоя концентрированным пучком энергии, который в нашем случае является отношением энергии, необходимой только для расплавления материала, к суммарной энергии, подаваемой пучком (в P следует принимать во внимание учет эффективности электронов ηe или коэффициент связи лазерного луча), рассчитывается как (4) где S - энтальпия расплавленного материала в сварном шве при повышении температуры между начальной температурой T0 и температурой плавления Tm: (5) где C - средняя удельная теплоемкость для температурного диапазона между данными значениями температуры, а Hf - теплота плавления. Хорошо видно, что термический КПД ηt можно оценить как отношение Y/X. В работе [21] этот подход применялся последовательно для прогнозирования оптимальных режимов при электронно-лучевой сварке. Экспериментальные данные В табл. 1 приведены теплофизические параметры порошковых материалов некоторых более важных металлических сплавов, используемых для моделирования параметров аддитивного производства. В табл. 2 приведены некоторые экспериментальные данные для режимов избирательного плавления электронным пучком и лазерным лучом, при этом значение X определяется как (6) Поглощающая способность поверхности (эффективность электронов для электронных пучков или коэффициент связи для лазерных лучей) принималась равной 0,6. Максимальная безразмерная температура (7) В табл. 1 значения теплопроводности λ для порошков составляют порядка 10 % от значений для сыпучего материала, а значения плотности ρ порошковых слоев предполагаются как 67-75 % от плотности объема наплавки. На рис. 3 приведены функция qm (rV/2a) и экспериментальные точки, представляющие режимы селективного плавления некоторых порошковых слоев промышленных сплавов. Наклонная линия на рис. 3 обозначает максимальный термический КПД 48,45 % [20]. Это теоретический предел для отношения количества поглощенной энергии, использованной для плавления материала при нагревании образцов путем перемещения линейного источника нагрева к общей поглощенной энергии. Рис. 3. Нормализованная диаграмма обработки из функции qm (rV/2a) и экспериментальные данные для селективного плавления: • - Ti-6Al-4V; * - SS 316L; x - Inconel 625; + - CM247 и o - FeCoCrNi. Ссылки цитирования (для этих данных) отмечены цифрами Таблица 1 Теплофизические параметры для порошковых слоев из сплавов, используемых для плавления энергией пучка Порошок l, Вт/(м·К) C, Дж/(кг·К) r, кг/м3 a = l/Cr, м2/с Tm, К Ti-6Al-4V 0,63 530 3000 0,396·10-6 1604 SS 316L 1,63 500 6000 0,543·10-6 1370 Inconel 625 2,52 620 6000 0,677·10-6 1605 CM 247 1,54 320 6000 0,802·10-6 1612 FeCoCrNi 1,8 450 5800 0,69·10-6 1600 Таблица 2 Экспериментальные данные для аддитивного производства металлических деталей в порошковом слое Сплав Пучок T0, К L, м rb / V, м, м/с P, Вт l Vrb/2a qm Ссылка Ti-6Al-4V Эл. пучок 923 - 150·10-6 /- 600 71,54 11,85 0,0878 [10] Ti-6Al-4V Эл. пучок 923 - 150·10-6 /- 600 71,54 16,3 0,0878 [10] Ti-6Al-4V Эл. пучок 923 50·10-6 150·10-6 /0,2 60 559,87 37,87 0,0112 [11] Ti-6Al-4V Эл. пучок 923 50·10-6 150·10-6 /6,4 1400 1306,38 212,21 0,048 [11] Ti-6Al-4V Лазер 298 30·10-6 52·10-6 /1,6 250 3505,94 104,51 0,0018 [12] Ti-6Al-4V Лазер 298 20·10-6 75·10-6 /0,8 150 1458,67 52,52 0,0043 [13] Ti-6Al-4V Лазер 298 20·10-6 75·10-6 /1,5 200 1944,9 142,04 0,003 23 [13] Ti-6Al-4V Лазер 498 30·10-6 70·10-6 /0,7 175 2152,76 92,8 0,002 92 [14] Окончание табл. 2 Сплав Пучок T0, К L, м rb /V, м, м/с P, Вт l Vrb/2a qm Ссылка Ti-6Al-4V Лазер 498 90·10-6 120·10-6 / 1,029 375 2690,76 92,8 0,002 33 [14] Ti-6Al-4V Эл. пучок 730 86·10-6 50·10-6 / 1,6 160 3479,02 101 0,0018 [25] Ti-6Al-4V Эл. пучок 730 86·10-6 200·10-6 / 0,8 240 2609,26 202,02 0,0024 [25] Ti-6Al-4V Эл. пучок 730 86·10-6 150·10-6 / 1,6 320 3219,34 303,03 0,0027 [25] Ti-6Al-4V Эл. пучок 730 86·10-6 200·10-6 / 1,6 480 2609,26 404,04 0,0024 [25] Ti-6Al-4V Эл. пучок 730 86·10-6 250·10-6 / 1,6 640 2782 505,05 0,002 25 [25] Ti-6Al-4V Эл. пучок 730 86·10-6 200´10-6 / 2,4 720 3913,9 606,06 0,0016 [25] Ti-6Al-4V Эл. пучок 730 86·10-6 150·10-6 / 4,0 800 5798,36 757,57 0,011 [25] Ti-6Al-4V Эл. пучок 730 86·10-6 150·10-6 / 5,6 1120 81 17,71 1060 0,000 774 [25] Ti-6Al-4V Эл. пучок 730 86·10-6 150·10-6 / 8,0 1600 11 596,72 1515,15 0,000 542 [25] SS316L Лазер 298 30·10-6 27·10-6 / 0,5 150 1910,83 12,424 0,003 29 [15] SS316L Лазер 298 30·10-6 27·10-6 / 2,5 400 2123,14 62,12 0,002 96 [15] SS316L Лазер 298 50·10-6 100·10-6 / 0,2 97 333,71 18,4 0,019 [16] Inconel625 Лазер 298 100·10-6 850·10-6 / 0,11 900 209,43 69,26 0,03 [17] CM247 Лазер 298 20·10-6 75·10-6 / 0,4 100 395,26 18,7 0,016 [18] CM247 Лазер 298 20·10-6 75·10-6 / 2 200 790,51 93,52 0,007 94 [18] FeCoCrNi Лазер 298 20·10-6 25·10-6 / 0,33 200 3015,075 5,98 0,002 08 [19] FeCoCrNi Лазер 298 50·10-6 25·10-6 / 0,66 200 1205,545 11,96 0,005 209 [19] Обсуждение методики Нормализованная диаграмма обработки на рис. 3 отображает безразмерную максимальную температуру и безразмерную скорость пучка. Безразмерные числа полезны, поскольку они уменьшают общее число переменных, подлежащих изучению, и обеспечивают важное понимание процесса, которое одна технологическая переменная не может обеспечить. Позиции точек измерения безразмерной обработки должны быть выбраны в окне параметров обработки, которое представляет собой область ниже функции qm (rV/2a), показанной на рис. 3. Область режима обработки, расположенная ниже линии максимальной безразмерной температуры (или максимального термического КПД), дает возможность прогнозировать возможные оптимальные параметры, используемые при селективном плавлении энергией пучка. Оптимальное место точки режима для процесса нагрева при селективном плавлении порошковых слоев с помощью энергии пучка - это положение как можно ближе к предельному значению, задаваемому функцией qm (rV/2a). Точки из рис. 3 с безразмерными параметрами режима, расположенными над этой линией, можно объяснить следующим образом. Точки из статьи [10] были рассчитаны авторами этой статьи, также в статье отсутствуют данные, что не позволяет понять причину этих расхождений. Точка для режима лазерного селективного плавления из статьи [17] представляет собой процесс со слишком большим пятном луча. Тепловая модель с движущимся линейным источником нагрева не дает хороших результатов в этом случае. Объяснение положения одной точки из режимов аддитивного производства с использованием электронного пучка из статьи [25], где rb составляет 250·10-6 м, может быть аналогичным. Процессы в жидком металле и размеры расплавленной ванны в этих случаях более важны для полученных результатов и оценки термического КПД процесса. Отметим, что одна и та же аналитическая модель нагрева последовательно применялась также в источниках [22-24] для оценки глубины сварного шва при термической модификации поверхности электронным пучком и при электронно-лучевой сварке тонких пластин при режимах, приводящих к сферическим сечениям сварных швов. С использованием значений известных qm для заданных значений P и L, при этом L - толщина свариваемых тонких пластин или длина сканирования пучком, кривая на рис. 2 дает возможность определить глубину нагрева определенного материала до некоторой температуры (например, температуры плавления при сварке или температуры фазового превращения при закалке). Выводы Аналитическая модель для нагрева твердых образцов движущимся линейным источником тепла, которая уже была изложена авторами для электронно-лучевой сварки [21], адаптирована для образцов, покрытых порошковым слоем. Построена и обсуждена нормализованная диаграмма обработки порошкового слоя концентрированными пучками энергии при аддитивном производстве. В сочетании с данными экспериментальных параметров, доступными в литературе, была исследована применимость этой диаграммы и продемонстрировано окно параметров обработки. Предполагается, что эти диаграммы обработки предоставят полезную информацию для практикующих ученых и инженеров при выборе соответствующих параметров обработки для экспериментов на ранних этапах программ исследований и разработок технологий.About the authors
I. Yu Letyagin
Perm National Research Polytechnic University
D. N Trushnikov
Perm National Research Polytechnic University
V. Ya Belenkiy
Perm National Research Polytechnic University; Perm Military Institute of National Guard Troops of the Russian Federation
A. P Erikov
Perm Military Institute of National Guard Troops of the Russian Federation
T. V Olshanskaya
Perm National Research Polytechnic University
References
- Korner C., Bauereiss A., Attar E. Fundamental consolidation mechanisms during selective beam melting of powders // Modell. and Simul. in Materials Sci. and Eng. - 2013. - Vol. 21, art. no. 085011. - P. 8.
- Klassen A., Bauereiss A., Korner C. Modelling of electron beam absorption in complex geometries // J. of Phys.: Appl. Phys. - 2014. - Vol. 47, iss. 6. - P. 065307.
- Overview of modelling and simulation of metal powder bed fusion process / W. King, A.T. Anderson, R.M. Ferencz, N.E. Hodge, C. Kamath, S.A. Khairallah // Lawrence Livermore National Laboratory Materials Sci. and Techn. (United Kingdom). - 2015. - Vol. 31. - P. 957-968.
- Additive manufacturing of metallic components-process, structure and properties / T. DebRoy, H.L. Wei, J.S. Zuback, T. Mukherjee, J.W. Elmer, J.O. Milewski, A.M. Beese, A. Wilson-Heid, A. De, W. Zhang // Progress in Materials Sci. - 2018. - Vol. 92. - P. 112-224.
- Contuzzi N., Campanelli S.L., Ludovico A.D. 3D-finite element analysis in the selective laser melting process // Int. J. of Simul. Modell. - 2011. - Vol. 10 (3). - P. 113-121.
- Manvatkar V., De A., DebRoy T. Spatial variation of melt pool geometry, peak temperature and solidification parameters during laser assisted additive manufacturing process // Materials Sci. and Eng. - 2015. - Vol. 31(8). - P. 924-930.
- Wen S.Y., Shin Y.C. Modeling of transport phenomena during the coaxial laser direct deposition process // J. Appl. Phys. - 2010. - Vol. 108. - P. 044908.
- Effect of fluid convection on dendrite arm spacing in laser deposition / Y. Lee, M. Nordin, S.S. Babu, D.F. Farson // Metall. and Materials Trans. - 2014. - Vol. B 45 (4). - P. 1520-1529.
- Korner C., Attar E., Heinl P. Mesoscopic simulation of selective beam melting processes // J. of Materials Proc. Techn. - 2011. - Vol. 211(6). - P. 978-987.
- Thomas M., Baxter G.J., Todd I. Normalised model-based processing diagrams for additive layer manufacture of engineering alloys // Acta Mater. - 2016. - Vol. 108. - P. 26-35.
- Processing window and evaporation phenomena for Ti-6Al-4V produced by selective electron beam melting / V. Juechter, T. Scharowsky, R. Singer, C. Korner // Acta Materialia. - 2014. - Vol. 76. - P. 252-258.
- Humbeeck Microstructure and mechanical properties of a novel β titanium metallic composite by selective laser-melting / B. Vrancken, L. Thijs, J.P. Kruth, J. Van // Acta Materialia. - 2014. - Vol. 68. - P. 150-158.
- Qiu C., Adkins N.J.E., Attallah M.M. Microstructure and tensile properties of selectively laser-melted and of HIPed laser-melted Ti-6Al-4V // Materials Sci. and Eng. - 2013. - Vol. A 578. - P. 230-239.
- Additive manufacturing of strong and ductile Ti-6Al-4V by selective laser melting via in situ martensite decomposition / W. Xu, M. Brandt, S. Sun, J. Elambasseril, Q. Liu, K. Latham, K. Xia, M. Qian // Acta Materialia. - 2015. - Vol. 85. - P. 74-84.
- Density of additively-manufactured, 316L SS parts using laser powder-bed fusion at powers up to 400 W. / C. Ka- math, B. El-dasher, G.F. Gallegos, W.E. King, A. Sisto // Int. J. of Adv. Manuf. Techn. - 2014. - Vol. 74. - P. 65-78.
- Application of X-ray CT method for discontinuity and porosity detection in 316L stainless steel parts produced with SLM technology / G. Ziolkowski, E. Chlebus, P. Szymczyk, J. Kurzac // Archives of Civil and Mech. Eng. - 2014. - Vol. 14, iss. 4. - P. 608-614.
- Additive layer manufacture of inconel 625 metal matrix composites, reinforcement material evaluation / D.E. Cooper, N. Blundell, S. Maggs, G.J. Gibbons // J. of Materials Proc. Techn. - 2013. - Vol. 213, no 12. - P. 2191-2200.
- The influence of the laser scan strategy on grain structure and cracking behaviour in SLM powder-bed fabricated nickel superalloy / L.N. Carter, C. Martin, P.J. Withers, M.M. Attallah // J. of Alloys and Comp. - 2014. - Vol. 615. - P. 338-347.
- Brif Y., Thomas M., Todd I. The use of high-entropy alloys in additive manufacturing // Scripta Materialia. - 2015. - Vol. 99. - P. 93-96.
- Рыкалин Н.Н. Расчеты тепловых процессов при сварке. - М.: Гос. науч.-техн. изд-во машиностроит. лит., 1951. - 296 с.
- Koleva E., Mladenov G., Vutova K. Calculation of weld parameters and thermal efficiency in electron beam welding // Vacuum. - 1999. - Vol. 53. - P. 67-70.
- Дворкин И.И., Ледовской В.П., Младенов Г.М. Электронно-лучевая сварка тонкостенных деталей в непрерывном режиме // Электронная техника. - 1972. - Сер. 4, № 3. - C. 53-59.
- Petrov P., Mladenov G. Theoretical analysis of heat flow and structural changes during electron beam irradiation of steel // Vacuum. - 1990. - Vol. 42, no 1/2. - P. 29-32.
- Koleva E., Mladenov G. Process parameter optimization and quality improvement at electron beam welding // Welding: Proc., Quality and Appl. - 2010. - Chapter 2. - P. 101-166. (Seria ″Mechanical Engineering-Theory and Applications″.)
- Macroscopic simulationand experimental measurement of melt pool characteristics in selective electron beam melting of Ti-6Al-4V / D. Riedlbauer, T. Scharowsky, R.F. Singer, P. Steinmann, C. Korner, J. Mergheim // Int. J. of Adv. Manuf. Techn. - 2017. - Vol. 88 (5-8). - P. 1309-1317.
Statistics
Views
Abstract - 54
PDF (Russian) - 29
Refbacks
- There are currently no refbacks.