Ratio of polymer temperature at the end of FFF/FDM 3D-printer nozzle to length of active (hot) part for various linear print speeds definition
- Authors: Matveev E.V1, Bezukladnikov I.I1, Trushnikov D.N1, Bogatyrev V.S1, Oskolkov A.A1
- Affiliations:
- Perm National Research Polytechnic University
- Issue: Vol 22, No 1 (2020)
- Pages: 70-78
- Section: ARTICLES
- URL: https://ered.pstu.ru/index.php/mm/article/view/2975
- DOI: https://doi.org/10.15593/2224-9877/2020.1.08
- Cite item
Abstract
This article describes the basic problems of FFF / FDM 3D-printing relates to the behavior of plastic melt in the hot part of the extruder. A more mobile and easily accessible technology for printing plastic products of complex geometric shapes based on FDM technology of 3D-printing developed by the scientific team of Perm National Research Polytechnic University is represented. A comparison of a standard extruder developed by the authors of the article is given. It is concluded that FDM technology has not fully revealed its potential due to a wide range of technical problems. The range of issues to be solved for the successful manufacture of plastic products using this technology is considered. Hypotheses and ways of solving problems are being put forward. Prospects of this technology are considered, as well as an assessment of its utility in production and for society. The methodology and results of modeling the behavior of the polymer melt inside the FFF/FDM 3D printer special shape nozzle are presented. SolidWorksSimulation was used as a modeling environment for determining the eveness of heating of the hot-end of a standard extruder and a new type extruder. To calculate the maximum polymer passage through the nozzle of the developed extruder under the condition of even heating of the hot-end (nozzle), an multi-physical axisymmetric mathematical model was constructed in the ComsolMultiphysics package, which containing a thermal and electromagnetic tasks. The simulation results allow to determine the ratio of print speed, length and temperature of the heated part for 2 types of polymer. This dependence can be used to control the temperature of the nozzle and the feed rate of the polymer thread directly during printing. Such regulation will ensure a higher quality of manufactured objects.
Keywords
FFF, FDM, FFF, FDM, 3D-printing, uneven heating, extruder, mathematical modeling, 3D-printer, printing quality, nozzle, filament.
Full Text
Введение В настоящее время различные технологии 3D-печати не просто активно внедряются в производственный процесс, а уже стали проверенным инструментом на современных малых и больших предприятиях. Одной из самых простых, надежных и старых аддитивных технологий является FDM/FFF 3D-печать, которая была разработана в конце 80-х гг. компанией Stratasys [1-11]. Несмотря на значительную популярность и распространение в производстве сегодня, сам технологический процесс «выращивания» изделия далек от идеала. Связано это с различными проблемами, возникающими при экструдировании расплавленного полимера через сопло [11]. Внешний вид экструдера представлен на рис. 1. Рис. 1. Стандартный экструдер для FDM/FFF 3D-печати Экструдер 3D-принтера состоит из двух компонентов: горячей и холодной частей. Холодная часть экструдера представляет собой подающий механизм, через который осуществляется подача филамента в горячую часть экструдера. Филамент представляет собой полимерную проволоку (пруток) диаметром 1,75 мм в твердом состоянии. Горячая часть экструдера представляет собой совокупность элементов, в которых происходит переход подаваемого твердого филамента в расплавленное состояние. На рис. 2 представлена структурная схема стандартного экструдера для FDM/FFF 3D-печати. Параметры такого экструдера приведены в табл. 1. Таблица 1 Параметры стандартного экструдера № п/п Параметр Значение 1 Тип нагрева Контактный 2 Тип нагревательного элемента Резистивный пальчиковый 3 Пиковая мощность нагревателя 75 Вт 4 Напряжение 12-24 В 5 V нагрева до 200 °С ~60 с 6 Максимальная T сопла 250 °С (при наличии тефлоновой трубки), 380 °С (при цельнометаллическом термобарьере без тефлоновой трубки) 7 Горячая масса 35 г Дополнительно стоит отметить, что масса полимера, которую необходимо постоянно поддерживать в состоянии расплава во время печати, при средней линейной скорости подачи филамента 30 мм/с составляет 1-1,7 г [12]. В среде моделирования Solid Works Simulation была проанализирована тепловая модель стандартного экструдера, представленная на рис. 3. Несмотря на незначительную тепловую массу (35 г), в данном экструдере присутствует нагревательный элемент, тепловая масса которого относительно всей системы в целом достаточно велика [13]. Следствием этого является то, что данная конструкция обладает значительной тепловой инерционностью и имеет ряд недостатков, связанных Рис. 2. Структурная схема стандартного экструдера для FDM/FFF 3D-печати Рис. 3. Тепловая модель стандартного экструдера с процессом экструдирования материала и самим процессом 3D-печати: 1) большая тепловая масса - неконтролируемо вытекающий полимер во время переходов или паузы, подтеки на углах из-за повышения температуры и понижения скорости; 2) низкоскоростной регулятор - отсутствие возможности поддерживать процесс экструзии с заданными параметрами при изменении внешних факторов (неравномерный диаметр прутка, изменение скорости подачи прутка); 3) сложность печати упругими филаментами, филаментами с высокой теплопроводностью и полимерами с высоким ТКЛР (температурным коэффициентом линейного расширения). Поскольку полимер обладает высокой теплопередачей, зона расплава увеличивается, подаваемая пластиковая нить расширяется и застревает в трубке термобарьера. В современных 3D-принтерах основными параметрами, отвечающими за экструзию полимера, являются температура сопла и скорость подачи филамента в экструдер. При увеличении скорости печати, а точнее увеличении объема проходящего через экструдер полимера за единицу времени, изменяется тепловой баланс расплава пластика. Ключевым фактором, обеспечивающим стабильное качество экструзии материала при повышении скорости подачи полимера в горячую часть экструдера и скорости его перехода из твердого состояния в состояние расплава, является сохранение равномерности прогрева в горячей части экструдера. При значительных скоростях и большом диаметре нагреваемой трубки термобарьера и камеры расплава температурное поле расплава не является равномерным [14]. В качестве альтернативы классическому FFF/FDM-экструдеру авторами был разработан экструдер нового типа, изображенный на рис. 4. Новизна данного устройства заключается в использовании бесконтактного способа нагрева горячей части экструдера и большой площади поверхности нагреваемой зоны, низкой тепловой массы горячей части экструдера, низкой массы расплава. Параметры такого экструдера приведены в табл. 2. Рис. 4. Разработанный экструдер Таблица 2 Параметры разработанного экструдера № п/п Параметр Значение 1 Тип нагрева Бесконтактный 2 Тип нагревательного элемента Токи высокой частоты 3 Пиковая мощность нагревателя 700 Вт 4 Напряжение 24 В 5 V нагрева до 200 °С ~2 с 6 Максимальая T сопла 1000 °С 7 Горячая масса 1,16 г Рис. 5. Тепловая модель разработанного экструдера В среде моделирования Solid Works Simulation также была проанализирована тепловая модель разработанного экструдера, представленная на рис. 5. Отличительными особенностями являются: высокий градиент температуры горячей части, увеличенная площадь контакта полимера с зоной нагрева (стенками сопла), равномерный прогрев сопла. Для расчета максимальной скорости прохождения полимера через сопло разработанного экструдера при условии равномерного прогрева горячей части (сопла) была построена мультифизическая осесимметричная математическая модель в пакете Сomsol Multiphysics [15-17], содержащая тепловую и электромагнитную задачи в формулировках, рассмотренных далее. Электромагнитная задача Для упрощения процесса моделирования будем считать, что нагрузка силового источника (индуктор) является согласованной, частота колебаний подстраивается автоматически и, соответственно, в каждый момент времени в витке индуктора протекает гармонический ток I известной величины. В этом случае электромагнитные процессы в такой системе описываются представленной ниже системой дифференциальных уравнений относительно комплексного векторного магнитного потенциала, имеющего единственную составляющую Для области пространства, занимаемой рабочим телом, (1) Для области пространства, содержащей i-й виток индуктора, (2) Для области пространства, занимаемой концентратором магнитного потока, (3) Для окружающего воздуха (4) Во всех вышеприведенных уравнениях σ - электропроводность материала (материала сопла 40Х13, материала катушки индуктора - меди, воздуха); ω - угловая частота; j - мнимая единица; µ0 и ε0 - магнитная и диэлектрическая проницаемости вакуума соответственно, µ0 = 1,256·10-6, ε0 = 8,85·10-12 Ф/м; Jφ - плотность тока в индукторе; r - радиальная координата; Si - площадь поперечного сечения индуктора; Sk - площадь поперечного сечения концентратора; Ui - первоначально неизвестная разность потенциалов, приложенная к i-му витку индуктора; Uk - наведенная разность потенциалов на концентраторе; I - электрический ток, протекающий в витках индуктора (задаваемая величина). Комплексная величина полного напряжения Uн, приложенного к индуктору, определяется как алгебраическая сумма значений напряжения на всех витках. Величина этого напряжения учитывается при расчете активной и реактивной мощностей, потребляемых рассматриваемой системой. Зависящее от частоты эквивалентное комплексное сопротивление такой системы определяется как . Граничными условиями рассматриваемой задачи, описываемой системой уравнений (1)-(4), можно считать условие симметрии относительно оси z и магнитной изоляции на внешних границах расчетной области Тепловая задача Индуктор в рассматриваемой системе охлаждается за счет теплового контакта с внешней цилиндрической оправкой, играющей роль корпуса, которая, в свою очередь, механически соединена с радиатором экструдера, имеющим принудительный обдув. Верхняя часть сопла (выше проточки) также является охлаждаемой за счет теплового контакта с соответствующей частью цилиндрической оправки. Для упрощения решаемой задачи будем считать, что указанные тепловые контакты являются идеальными, а охлаждение экструдера достаточно для обеспечения охлаждения индуктора. Соответственно, температура в их объеме является постоянной и не превышает величины T = 50 °С. При таком допущении математическая модель нестационарного теплового процесса индукционного нагрева заготовки состоит всего из двух уравнений: Здесь p, Cp, λ - плотность, удельная теплоемкость и теплопроводность материала соответствующей среды, зависящие от температуры Т, p = p(T), Cp = Cp(T), λ = λ(T); Q - удельная мощность источника тепла, обеспечивающая индукционный нагрев рабочего узла вихревыми токами, (5) где J - действующее значение плотности наведенного тока в объеме ферромагнитного сопла; Aφ - комплексно-сопряженное значение магнитного потенциала. Передачей тепла в виде конвекции и излучения в системе можно пренебречь в силу малого вклада этих механизмов по сравнению с теплопередачей в результате теплопроводности материалов. Граничным условием для данной задачи можно принять условие Неймана dT/dr = 0 на оси симметрии. В области индуктора и концентратора, имеющих принудительное охлаждение, считаем T = const = 50 °С. На внешних границах расчетной области задаем условие теплообмена с окружающей средой где k - коэффициент теплоотдачи; T0 - температура окружающей среды; n - вектор внешней нормали к границе. Для представления указанных задач в интерфейсе Comsol Multiphysics можно воспользоваться стандартными физическими интерфейсами Magnetic Fields (mf) и Heat Transferin Solids (ht). Обе задачи могут решаться совместно в рамках одной модели за счет их связи посредством джоулева тепла, определяемого по формуле (5) [18, 19]. Исходные данные для моделирования приведены в таблице ниже. Начальные условия приняты в виде выражений . В ходе моделирования решалась параметрическая задача определения соотношения температуры нагрева экструдируемого расплава на выходе сопла и длины активной (горячей) части для различных линейных скоростей печати. Пересчет линейной скорости печати в скорость подачи материала на входе экструдера производился для диаметра выходного отверстия сопла 0,4 мм и диаметра материала на входе в экструдер 1,75 мм. Параметры: полимер - PC, PA6; диаметр прутка 1,75 мм; выходное отверстие 0,4 мм. Результаты моделирования для PС представлены на рис. 6. Результаты моделирования для PA6 представлены на рис. 7. По результатам моделирования были составлены графики. По оси Х представлена длина активной части (нагреваемая часть сопла) в миллиметрах, по оси Y представлена температура расплава полимера в кельвинах. Кривые на графике отражают зависимость сохранения равномерной заданной температуры расплава полимера внутри сопла от определенной длины и скорости экструзии полимера. Исходные данные модели Индуктор Материал медь: электропроводность - 5,998·107 См/м; отн. магнитная проницаемость - 1; отн. диэлектрическая проницаемость - 1 Число витков 26 Внутренний диаметр, мм 10 Внешний диаметр, мм 14 Высота, мм 16 Напряжение источника, В 24 Частота колебаний, кГц 100-150 Ток, А 7 Цилиндрическая оправка, корпус экструдера Материал алюминий: электропроводность - 3,774·107 См/м; отн. магнитная проницаемость - 1; отн. диэлектрическая проницаемость - 1 Внутренний диаметр, мм 15 Внешний диаметр, мм 17 Высота, мм 17 Рабочее тело (сопло) экструдера Материал сталь 40Х13: электропроводность - 1,12·107 См/м; отн. магнитная проницаемость - 380; отн. диэлектрическая проницаемость - 1 Высота активной части, мм 0-100 Внутренний диаметр, мм 2 Внешний диаметр активной части, мм 4 Диаметр выходного отверстия, мм 0,4 Внешний диаметр термобарьера, мм 3 Высота термобарьера, мм 2 Внешний диаметр хвостовика сопла, мм 6 Высота хвостовика сопла, мм 10 Дополнительные данные для решения тепловой задачи Исходная температура окружающей среды, °С 20 Коэффициент теплоотдачи, Вт/(м2·К) 15 Рис. 6. Соотношение скорости печати и температуры нагрева РС во время экструзии Рис. 7. Соотношение скорости печати и температуры нагрева PA6 во время экструзии Заключение В ходе исследования было проведено сравнение стандартного и разработанного авторским коллективном экструдеров. Результатом стало получение зависимости максимальной допустимой скорости выходного потока расплавленного полимера от длины нагреваемой части (сопла) разработанного экструдера при обеспечении равномерности и стабильности заданного поля температуры. Данная зависимость может быть использована для осуществления регулирования температуры сопла и скорости подачи полимерной нити непосредственно во время печати. Такое регулирование позволит обеспечить более высокое качество изготавливаемых объектов.About the authors
E. V Matveev
Perm National Research Polytechnic University
I. I Bezukladnikov
Perm National Research Polytechnic University
D. N Trushnikov
Perm National Research Polytechnic University
V. S Bogatyrev
Perm National Research Polytechnic University
A. A Oskolkov
Perm National Research Polytechnic University
References
- Chennakesava P., Narayan Y.S. Fused deposition modeling-insights // International Conference on Advances in Design and Manufacturing (ICAD&M'14), Tiruchirappalli, 7 December 2014. - Tiruchirappalli, 2014. - P. 1345-1350.
- Fused deposition modelling: a technology evaluation / T. Grimm [et al.] // Time-compression Technologies. - 2003. - Vol. 11, no. 2. - P. 1-6.
- Stratasys. - URL: https://www.stratasys.com/ (accessed 01 February 2020).
- Structural quality of parts processed by fused deposition / M.K. Agarwala, V.R. Jamalabad, N.A. Langrana, A. Safari, P.J. Whalen, S.C. Danforth // Rapid Prototyping J. - 1996. - Vol. 2. - P. 4-19.
- Anisotropic material properties of fused deposition modeling ABS / S.H. Ahn, M. Montero, D. Odell, S. Roundy, P.K. Wright // Rapid Prototyping J. - 2002. - Vol. 8. - P. 248-257.
- Bellini A., Guceri S. Mechanical characterization of parts fabricated using fused deposition modeling // Rapid Prototyping J. - 2003. - Vol. 9. - P. 252-264.
- Ji L.B., Zhou T.R. Finite element simulation of temperature field in fused deposition modeling // Manufacturing Sci. and Eng. - 2010. - Vol. 97/101. - P. 2585-2588.
- Sun Q., Rizvi G.M., Bellehumeur C.T., Gu P. Effect of processing conditions on the bonding quality of FDM polymer filaments // Rapid Prototyping J. - 2008. - Vol. 14. - P. 72-80.
- Thomas J.P., Rodriguez J.F. Modeling the fracture strength between fused deposition extruded roads // Solid freeform Fabrication Proceedings / University of Texas at Austin. - Austin, TX, 2000. - P. 16-23.
- Thermal analysis of fused deposition / M.A. Yardimci, T. Hattori, S.I. Guceri, S.C. Danforth // Solid Freeform Fabrication Proceedings / eds. D.L. Bourell, J.J. Beaman, R.H. Crawford, H.L. Marcus, J.W. Barlow; University of Texas at Austin. - Austin, TX, 1997. - Р. 689-698.
- Turner B.N., Strong R., Gold S.A. A review of melt extrusion additive manufacturing processes: I. Process design and modeling // Rapid Prototyp. J. - 2014. - Vol. 20 (3). - P. 192-204.
- Effects of extrusion speed and printing speed on the 3D printing stability of extruded PEEK filament / P. Geng [et al.] // J. of Manufacturing Processes. - 2019. - Vol. 37. - Р. 266-273.
- Temperature analysis of nozzle in a FDM type 3D printer through computer simulation and experiment / J.H. Park [et al.] // Elastomers and Composites. - 2016. - Vol. 51, no. 4. - Р. 301-307.
- Heating and flow computations of an amorphous polymer in the liquefier of a material extrusion 3D printer / F. Pigeonneau [et al.] // Additive Manuf. - 2020. - Vol. 32. - Р. 1-16.
- Calibration of a numerical model for heat transfer and fluid flow in an extruder / T. Hofstätter, D.B. Pedersen, J.S. Nielsen, R. Pimentel, M. Mischkot // Int. J. of Rapid Manuf. - 2016. - Vol. 6 (1). - Р. 1-16.
- Иродов И.Е. Основные законы электромагнетизма: учеб. пособие для студ. вузов. - 2-e изд., стереотип. - М.: Высш. шк., 1991. - 288 с.
- Comsol A.B. AC/DC Module-user’s guide // COMSOL. - 2011. - Vol. 3. - P. 151.
- The classification of coupled field problems / K. Hameyer, J. Driesen, H. de Gersem, R. Belmans // IEEE Trans. on Magnetics. - 1999. - Vol. 35, no. 3. - P. 1618-1621.
- Applications of coupled field formulations to electrical machinery / G.B. Kumbhar, S.V. Kulkarni, R. Escarela-Perez, E. Campero-Littlewood // Int. J. for Computation and Mathematics in Electrical and Electronic Engineering (COMPEL). - 2007. - Vol. 26, no. 2. - P. 489-523.
Statistics
Views
Abstract - 130
PDF (Russian) - 29
Refbacks
- There are currently no refbacks.