Next Article in Journal
Estimation of the Effects of CO2 and Temperature on the Swelling of PS-CO2 Mixtures at Supercritical Conditions on Rheological Testing
Next Article in Special Issue
Analysis of Thermomechanical Properties and the Influence of Machining Process on the Surface Structure of Composites Manufactured from Metal Chips with a Polymer Matrix
Previous Article in Journal
Effect of Coir Fiber Surface Treatment on Interfacial Properties of Reinforced Epoxy Resin Composites
Previous Article in Special Issue
Compression and Bending Properties of Short Carbon Fiber Reinforced Polymers Sandwich Structures Produced via Fused Filament Fabrication Process
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Polymers
Volume 14
Issue 17
10.3390/polym14173482
polymers-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Experimental Validation of Numerical Model for Thermomechanical Performance of Material Extrusion Additive Manufacturing Process: Effect of Process Parameters

by
Ans Al Rashid
1,* and
Muammer Koç
1,2
1
Division of Sustainable Development, College of Science and Engineering, Hamad Bin Khalifa University, Qatar Foundation, Doha 34110, Qatar
2
Faculty of Engineering, University of Karabük, Karabük 78050, Turkey
*
Author to whom correspondence should be addressed.
Polymers 2022, 14(17), 3482; https://doi.org/10.3390/polym14173482
Submission received: 26 July 2022 / Revised: 11 August 2022 / Accepted: 16 August 2022 / Published: 26 August 2022
(This article belongs to the Special Issue Multifunctional Polymer Molding Processing of Polymers)
Browse Figures
Versions Notes

Abstract

:
The material extrusion additive manufacturing (MEAM) process for polymers seems straightforward. However, several controlled and uncontrolled factors affect the 3D printed product quality, e.g., MEAM process parameters, thermomechanical properties of the material, and part design. Therefore, it is crucial to understand these interlinked factors of part geometry, material properties, and 3D printing (3DP) process parameters to optimize 3D printed product quality. The numerical models and simulation tools can predict the thermomechanical performance of the MEAM process under given input parameters (material, design, and process variables) and reduce the research and development costs significantly. However, the numerical models and tools need further exploration and validation of simulation predictions for their adaptability and reliability. Therefore, in this study, numerical simulations were performed to observe the impact of process parameters on the part quality of MEAM 3D printed components. The two crucial process parameters (i.e., extrusion temperature and layer resolution) were varied while keeping the other process parameters, part geometry (tensile testing coupon), and material properties (acrylonitrile butadiene styrene (ABS)) constant. These two process parameters were sequentially optimized for optimum part quality, first by varying the extrusion temperature and secondly by changing the printing resolution using the optimum printing temperature. The 3DP process quality was evaluated in terms of dimensional accuracy, distortions, and residual stresses. Finally, the specimens were 3D printed under similar process conditions to validate the numerical model predictions.

1. Introduction

The material extrusion additive manufacturing (MEAM), also referred to as the fused filament fabrication (FFF) process, is one of the most widely explored and adopted additive manufacturing (AM) or 3D printing (3DP) processes owing to its simplicity, cheaper hardware, and open-source technology [1,2]. According to a study, in 2016, 96% of the global AM sales market was shared by the FFF-based 3D printers [3]. In this process, the 3D printer is fed with the material in filament form, heated to a semi-solid state, and extruded on a build plate through a nozzle [4]. The material is extruded layer-by-layer until a solid 3D part is fabricated, and the extrudate cools down on the build plate once the process is completed [5]. Thermoplastics are mainly processed using this process. However, the technology has recently evolved to use polymer composites [6] and metals [7].
Although the MEAM process for polymers seems straightforward, several controlled and uncontrolled factors affect the 3D printed product quality [8]. These factors include MEAM process parameters (e.g., nozzle diameter, deposition speed, layer resolution, raster pattern, raster width, infill type and density, material flow rate, raster orientation, building orientation, extrusion temperature, ambient temperature, and build plate temperature), thermomechanical properties of the material (e.g., thermal conductivity, heat capacity, young modulus, coefficient of thermal expansion, and stiffness), and part design [9,10,11]. Therefore, it is crucial to understand these interlinked factors of part geometry, material properties, and 3DP process parameters to optimize 3D printed product quality [12]. Several studies reported the experimental evaluation of the effect of these factors on the 3DP process quality [13,14,15]. Although several studies followed the experimental approach for this cause, the numerical modeling and simulations in this regard are very few.
The numerical models and simulation tools can predict the thermomechanical performance of the MEAM process under given input parameters (material, design, and process variables) and reduce the research and development costs significantly [16]. Computational fluid dynamics (CFD) models have been used to simulate the 3DP process and to analyze the temperature profile [17], material melting, layer adhesion, nozzle pressure, material flow, and heat transfer using different commercial software [18,19,20,21,22,23,24,25,26,27]. However, the numerical models and tools need further exploration and validation of simulation predictions for their adaptability and reliability.
In the view of existing literature, it is vital to understand and evaluate the effect of MEAM process parameters on the 3D printed product quality. The reported study is part of our ongoing research where all aspects of the 3DP process, i.e., material properties, process parameters, and the design of specimens, are being considered. Other aspects will also be covered in our upcoming studies. However, the effect of extrusion temperature and layer height is investigated in this study. Therefore, numerical simulations were performed to observe the impact of process parameters on the part quality of MEAM 3D printed components. Different extrusion temperatures (210 °C, 220 °C, 230 °C, 240 °C, 250 °C, and 260 °C) within the printability window of acrylonitrile butadiene styrene (ABS) material and three different layer heights (0.1 mm, 0.15 mm, and 0.2 mm) were selected for analysis. These two process parameters were sequentially optimized for optimum part quality, first by varying the extrusion temperature and secondly by changing the printing resolution using the optimum printing temperature. The 3DP process quality was evaluated regarding dimensional accuracy, distortions, and residual stresses. Finally, the specimens were 3D printed under similar process conditions to validate the numerical model predictions.

2. Methodology

The current work is divided into two phases: the MEAM process numerical simulations and validation of numerical model results via 3DP experimental measurements. The main aim of this study is to evaluate the effect of process parameters on 3DP process quality quantified in terms of part deflections, distortions, and residual stresses. The extrusion temperature and printing resolution were varied while keeping the other process parameters, part geometry, and material properties constant. A standardized tensile testing coupon (according to ASTM D638 Type-I [28]) was selected as reference geometry to be fabricated using acrylonitrile butadiene styrene (ABS) material (from Ultimaker, The Netherlands). ABS provides flexibility in synthesizing nanocomposites due to its solubility in acetone [29] and possesses significant physical, mechanical, and chemical resistance properties [30,31]. A 3D CAD model for a tensile testing coupon was generated using Solidworks 2021 (Dassault Systèmes, Paris, France), as shown in Figure 1, and imported to slicing software. The slicing information (gcode) was generated using Cura software version 5.0.0 (from Ultimaker, The Netherlands), containing information on toolpath, printing resolution, extrusion temperature, build plate temperature, printing speed, and infill type and density. The subsequent sections present the methods adopted for numerical simulations and experimentation.

2.1. MEAM Process Simulations

The MEAM process simulations were performed using Digimat software (version 2021.3, from e-Xstream engineering, Käerjeng, Luxembourg). Chamber size of 215 mm (X) × 215 mm (Y) × 200 mm (Z) was defined with a moving build plate corresponding to the Ultimaker 3 Extended 3D printer (from Ultimaker, The Netherlands). The software utilizes a thermomechanical numerical model to calculate the temperature variations during the MEAM process [32]. Heat transfer equations evaluate the time-dependent temperature fields layer-by-layer during the material deposition. The built plate and ambient temperatures are used as boundary conditions for bottom surface elements and ambient environment for the overall model, respectively. The first deposited layer is constrained in all degrees of freedom due to the build plate adhesion, which gets released at the end of the 3DP process, and final residual stresses and deformations are estimated in this state.
The thermal gradients due to convective and radiative heat transfer cause the material shrinkage, which results in part deformations and distortions [33]. However, the residual stresses arise due to restricted molecular movements due to structural constraints [34]. The 3D model of the tensile testing coupon was imported in stl format, and thermomechanical properties (for ABS material) were defined from the software database. The numerical model requires toolpath information for 3D printed specimens imported as a gcode file. Process parameters were specified, including extrusion temperature (210 °C, 220 °C, 230 °C, 240 °C, 250 °C, and 260 °C), build plate temperature (80 °C), bead width (0.4 mm), chamber temperature (25 °C), convection coefficient (0.015 mW/mm2·°C), and ambient temperature (25 °C). Except for extrusion temperature and printing resolution, the other process parameters were constant for all the simulations. Voxel size (0.1, 0.15, or 0.2) was selected per the layer heights. Digimat software provides two options for geometry discretization, i.e., filament or layer-by-layer. Although filament discretization offers more precise results, the layer-by-layer technique was used to lower the computational time and cost.

2.2. MEAM Experiments

Ultimaker 3 Extended 3D printer and ABS material in filament form (from Ultimaker, The Netherlands) were used for 3DP experiments. The slicing information (gcode) was generated using Cura software (from Ultimaker, The Netherlands), containing information on toolpath, printing resolution, extrusion temperature, build plate temperature, printing speed, and infill type and density. Only two process parameters, i.e., extrusion temperature and printing resolution, were varied for their sensitivity to product quality. Different extrusion temperatures (210 °C, 220 °C, 230 °C, 240 °C, 250 °C, and 260 °C) within the printability window of ABS material and three different layer heights (0.1 mm, 0.15 mm, and 0.2 mm) were selected to fabricate 3D printed samples. These two process parameters were sequentially optimized for optimum part quality, first by varying the extrusion temperature and secondly by changing the printing resolution using the optimum printing temperature. The specimens were manufactured at a speed of 60 mm/s, using the build plate temperature of 80 °C, with two wall layers (corresponding to 0.8 mm of thickness as bead width was 0.4 mm). The top and bottom layers were 1.2 mm thick (where the number of layers varied depending upon the layer height), having 100% infill density with a zig-zag infill pattern.
The physical measurements were performed on the 3D printed samples for their dimensions and distortions. The dimensions were recorded using a vernier caliper, and at least three measurements were recorded for each targeted dimension. The measurements were compared with the numerical model predictions and designed part dimensions. Finally, the variations in residual stresses due to changes in 3DP process parameters were observed, and their effect on part quality was analyzed.

3. Results and Discussions

The numerical model performs the MEAM process simulation layer-by-layer to mimic the actual process, where a set of elements in the printing direction are activated in a time step. Depending upon the discretization approach, a group of elements is activated at a temperature selected for the extrusion process (ranging from 210 °C to 260 °C in this case). Due to convective and radiative heat losses, the extrudate temperature is significantly higher than the already deposited material [22]. As the printing process progresses, the deposited material cools down close to the printing platform temperature; though, the upcoming extrudate re-heats the already deposited material [35]. Due to the continuous thermal cycling of the material and heat transfer phenomena in the ambient environment, the thermal gradients are observed in the build direction, resulting in thermal stresses within the structure [33]. The temperature variations and material properties govern the thermal strains within the 3D printed part due to induced thermal stresses. The thermal strains result in dimensional variations of the 3D printed part from the designed dimensions [36]. Hence, the impact of 3DP process parameters, design, and material properties need to be investigated and understood to obtain desired structural integrity. The process parameters are sequentially optimized for ABS material to achieve optimum 3DP product quality by adjusting the extrusion temperature and printing resolution. The numerical model predictions were validated with 3DP experiments performed under the same process conditions.

3.1. Effect of Extrusion Temperature

The numerical model estimates the deviation of 3D printed sample dimensions from the designed values in three principal directions. The overall maximum displacements for the 3D printed samples are reported in Figure 2. The extrusion temperature varied from 210 °C to 260 °C within the processability window of the ABS material. The overall maximum displacement of 1.478 mm, 1.556 mm, 1.634 mm, 1.711 mm, 1.789 mm, and 1.867 mm was observed for specimens 3D printed at 210 °C, 220 °C, 230 °C, 240 °C, 250 °C, and 260 °C, respectively. In numerical predictions, a maximum deflection from designed dimensions was observed for specimen 3D printed at 260 °C. However, a minimum deviation was achieved for extrusion temperature of 210 °C. The coefficient of thermal expansion plays a vital role in material shrinkage, and the numerical simulation results also showed a prominent effect of printing temperature. In addition, the deformation trend for all the printing temperatures was the same, i.e., maximum at the specimen edges and minimum at the central region, owing to non-uniform cooling of the specimen edges.

3.1.1. Dimensional Analysis

The numerical simulation tool does not directly provide the dimensions of the simulated parts. Therefore, warped geometries were exported for dimensional analysis using Solidworks software 2021 (Dassault Systèmes, Paris, France). The targeted dimensions were measured as estimated by the numerical model and compared with the physical measurements taken for the 3D printed samples, as reported in Figure 3.
Table 1 comprehensively compares measurements recorded for the samples fabricated at different extrusion temperatures. Generally, the material shrinkage was observed in all the numerical simulation results in three principal directions. In numerical simulation results, a maximum deviation in overall length from the designed value of the specimen was 3.47 mm for 260 °C extrusion temperature. However, a minimum variation of 2.76 mm was achieved for specimen 3D printed at 210 °C. A maximum deviation of 0.32 mm in grip width from the designed value of the specimen was observed for 260 °C extrusion temperature. However, a minimum deviation of 0.25 mm was achieved for specimen 3D printed at 210 °C. Similarly, 0.15 mm and 0.18 mm of variation from the designed value of gauge width was measured for extrusion temperatures of 210 °C and 260 °C. Finally, the specimen thickness of 3.89 mm was achieved at printing temperatures of 210 °C and 220 °C, although a thickness of 3.86 mm was recorded for samples 3D printed at 260 °C. However, the numerical model predicted higher deflections from the designed values than the experimental measurements. In addition, a different trend in the dimensional measurements was observed, where the overall 3D printed specimen at 240 °C revealed better dimensional control over the rest of the extrusion temperatures considered.
In experimental measurements, a maximum deviation in overall length from the designed value of the specimen was 1.24 mm for 210 °C extrusion temperature. However, a minimum variation of 0.04 mm was achieved for specimen 3D printed at 240 °C. A maximum deviation of 0.12 mm in grip width from the designed value of the specimen was observed for 260 °C extrusion temperature. However, a minimum deviation of 0.05 mm was achieved for specimen 3D printed at 240 °C. Similarly, 0.05 mm and 0.15 mm of variation from the designed value of gauge width was measured for extrusion temperatures of 240 °C and 210 °C. Finally, the specimen thickness of 3.98 mm was achieved at printing temperatures of 240 °C. However, a thickness of 3.90 mm was recorded for samples 3D printed at 210 °C.
Overall, the numerical model predicted the measurements near the experimental dimensions. However, the higher deflections from the designed values with increasing extrusion temperatures are attributed to the material model limitation in the existing model. The material model currently does not consider the temperature dependence of the material properties, resulting in a linear trend with decaying product quality due to increased extrusion temperature. Secondly, higher variations in numerical model predictions are observed due to higher material shrinkages. This can be attributed to the stress relaxation phenomena in the 3D printed samples, which lowers the residual stresses and deflections; nevertheless, the stress relaxation is not considered in the current numerical model; hence higher shrinkages are observed. Secondly, the printing temperature also played an essential role in dimensional control of the 3D printed specimens, as higher variations are observed at higher extrusion temperatures. The experimental analysis concluded that optimum control over dimensions is achieved at an extrusion temperature ranging from 240 °C to 250 °C.

3.1.2. Warpage Analysis

Besides dimensional analysis, the out-of-plane deformations of the 3D printed parts are essential to analyze and are termed as distortions or warpage. The warpages come into action due to the thermal gradient in the build direction. In addition, the specimen edges go through non-uniform cooling, resulting in higher thermal stresses in these regions. The warpage measurements recorded from the numerical model predictions and 3D printed samples are also reported in Table 1. A visual comparison of warpages for simulation and experimental results is reported in Figure 4. The numerical model results recorded maximum and minimum warpage of 0.17 mm and 0.13 mm for specimens fabricated at 260 °C and 210 °C, respectively.
Although the differences in warpage measurements are a fraction of millimeters, they are important to consider where tighter control over dimensions is desired. The measured warpage in the 3D printed samples revealed higher values than the numerically predicted results. The imperfect build plate adhesion resulted in this variation, as the numerical model assumes perfect adhesion to the build plate. However, thermal stresses developed due to non-uniform cooling of the specimen edges, resulting in part detachment from the build plate. Therefore, higher z-direction deflections are observed. The numerical model predictions show that minimum warpage is achieved at an extrusion temperature of 210 °C. However, in experimental measurements, minimum warpage is achieved at 240 °C. The numerical model revealed a similar trend to deflections, where the warpages increase with increased extrusion temperatures. This is once again attributed to the existing material model limitations. Contrary to experimental measurements, higher warpages at higher or lower temperatures than the optimum temperature are due to two main reasons. Once the material is deposited at a temperature lower than the optimum temperature of 240 °C, it tends to cool down to the build plate or ambient temperature, depending on the layer number [37,38]. As the material cools down, the phase-change process from semi-solid to solid starts due to lower extrusion temperature, and inter-bead diffusivity of polymer chains is inhibited, resulting in imperfect adhesion. However, higher thermal cycling increases residual strains and warpages at higher temperatures.

3.1.3. Residual Stresses

The numerical model also predicted the residual stresses generated within the 3D printed part due to thermal cycling and process constraints. The numerically predicted residual stresses for the specimens manufactured at different extrusion temperatures are presented in Figure 5. The thermal gradient along the z-direction and non-uniform cooling of the samples results in the residual stresses.
Generally, higher residual stresses are generated along the specimen edges, and significantly lower residual stresses are observed within the specimen structure. In addition, the specimen 3D printed at an extrusion temperature of 210 °C revealed lower residual stresses than all the other printing temperatures under consideration. This observation also supports numerical predictions of the higher dimensional control at an extrusion temperature of 210 °C, as higher residual stresses result in higher dimensional variations and distortions in the specimen. The residual stresses can be limited by controlling the cooling process. Hence, the 3DP process within the temperature-controlled closed chamber can provide better quality.

3.2. Effect of Layer Height

In the second phase, only numerical simulations were performed to evaluate the effect of layer height on 3DP quality. Three different layer heights (0.10 mm, 0.15 mm, and 0.20 mm) were selected with an optimum extrusion temperature of 240 °C, while all the other process parameters were kept constant. The voxel sizes varied accordingly to mimic the layering effect and height. The process simulation for 0.20 mm layer height was performed in the previous phase. Only two simulations for 0.10 mm and 0.15 mm of layer resolution were performed, and results are reported in the subsequent sections.
The dimensional and distortion analysis was performed on the numerically predicted results, as reported in Table 2. A significant improvement in dimensional control was achieved by reducing the layer height. A comparison of dimensions for specimens virtually fabricated at layer heights of 0.15 mm and 0.10 mm is presented in Figure 6. An overall length of 162.26 mm and 165.00 mm, grip width of 18.77 mm and 18.95 mm, the gauge width of 12.88 mm and 12.92 mm, and thickness of 3.91 mm and 3.96 mm was achieved for 0.15 mm and 0.10 mm of layer heights, respectively. Finally, a warpage of 0.13 mm and 0.08 mm was observed for specimens modeled at 0.15 mm and 0.10 mm of printing resolution, as reported in Figure 7.
The results show that, at higher printing resolutions (lower layer heights), higher dimensional accuracy of the product is achieved. The 0.10 mm layer height produced specimens with final dimensions close to desired ones. As the layer height reduces (increasing resolution), the number of layers producing a 4 mm thickness increases. Therefore, the cooling of a layer with 0.10 mm has enough time to adequately bond to the already deposited layer and solidifies at the deposited position before the deposition of the incoming layer. Consequently, improved dimensional accuracy is achieved.

4. Conclusions

In this study, numerical simulations were performed to observe the impact of process parameters on the part quality of MEAM 3D printed components. Different extrusion temperatures (210 °C, 220 °C, 230 °C, 240 °C, 250 °C, and 260 °C) within the printability window of ABS material and three different layer heights (0.1 mm, 0.15 mm, and 0.2 mm) were selected for analysis. These two process parameters were sequentially optimized for optimum part quality, first by varying the extrusion temperature and secondly by changing the printing resolution using the optimum printing temperature. The 3DP process quality was evaluated regarding dimensional accuracy, distortions, and residual stresses. Finally, the specimens were 3D printed under similar process conditions to validate the numerical model predictions.
The numerical model predicted the actual 3DP process outcomes with a reasonable difference from the experimental results. However, the higher deflections from the designed values with increasing extrusion temperatures are attributed to the material model limitation in the existing model. The higher dimensional variations in numerical model results can also be attributed to the stress relaxation effect, which is not considered in the numerical model. However, higher warpages/distortions recorded in 3D printed samples were due to the numerical model considering perfect build plate adhesion of the specimen. The 3DP experimental results conclude that the optimum print quality can be achieved at an optimum extrusion temperature of 240 °C. Two main reasons are that higher-dimensional variations and warpages at higher or lower temperatures than the optimum temperature. Once the material is deposited at a temperature lower than the optimum temperature of 240 °C, it tends to cool down to the build plate or ambient temperature, depending on the layer number. As the material cools down, the phase-change process from semi-solid to solid starts due to lower extrusion temperature, and inter-bead diffusivity of polymer chains is inhibited, resulting in imperfect adhesion. However, higher thermal cycling increases residual strains and warpages at higher temperatures. The higher thermal gradients within the printed part can also lead to higher distortions in the build direction.
In addition, the specimen 3D printed at an extrusion temperature of 210 °C revealed lower residual stresses than all the other printing temperatures under consideration. This observation also supports numerical predictions of the higher dimensional control at an extrusion temperature of 210 °C, as higher residual stresses result in higher dimensional variations and distortions in the specimen. Finally, the printing resolution also plays a vital role in the dimensional control of 3D printed parts. Based on the numerical model predictions, it is suggested that higher printing resolutions should be used to produce higher quality parts, especially in functional components where higher dimensional control is desired.
It is concluded that extrusion temperature and process resolution are of paramount concern while achieving optimum product quality. The optimum printing temperature varies with material compositions. However, lower layer heights are suggested for higher dimensional accuracy of 3D printed parts. The presented numerical model can be further improved to optimize the printing temperature and layer heights for other materials. In future studies, other polymeric materials will be considered to optimize the process parameters for the MEAM process. The proposed methodology can be implemented at a lower cost and limit the need for experimental optimization of 3DP process parameters.

Author Contributions

Conceptualization, A.A.R. and M.K.; methodology, A.A.R.; software, A.A.R.; validation, A.A.R.; formal analysis, A.A.R.; investigation, A.A.R.; resources, M.K.; data curation, A.A.R.; writing—original draft preparation, A.A.R.; writing—review and editing, A.A.R. and M.K.; visualization, A.A.R.; supervision, M.K.; project administration, M.K.; funding acquisition, M.K. All authors have read and agreed to the published version of the manuscript.

Funding

APC was funded by the College of Science and Engineering, Hamad Bin Khalifa University, Qatar Foundation, Doha 34110, Qatar.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Yilmaz, B.; Al Rashid, A.; Mou, Y.A.; Evis, Z.; Koç, M. Bioprinting: A review of processes, materials and applications. Bioprinting 2021, 23, e00148. [Google Scholar] [CrossRef]
  2. Thiam, B.G.; El Magri, A.; Vanaei, H.R.; Vaudreuil, S. 3D Printed and Conventional Membranes—A Review. Polymers 2022, 14, 1023. [Google Scholar] [CrossRef] [PubMed]
  3. Torrado, A.R.; Shemelya, C.M.; English, J.D.; Lin, Y.; Wicker, R.B.; Roberson, D.A. Characterizing the effect of additives to ABS on the mechanical property anisotropy of specimens fabricated by material extrusion 3D printing. Addit. Manuf. 2015, 6, 16–29. [Google Scholar] [CrossRef]
  4. Al Rashid, A.; Khan, S.A.; Al-Ghamdi, S.G.; Koç, M. Additive manufacturing: Technology, applications, markets, and opportunities for the built environment. Autom. Constr. 2020, 118, 103268. [Google Scholar] [CrossRef]
  5. Ikram, H.; Al Rashid, A.; Koç, M. Additive manufacturing of smart polymeric composites: Literature review and future perspectives. Polym. Compos. 2022. [Google Scholar] [CrossRef]
  6. Al Rashid, A.; Ikram, H.; Koç, M. Additive manufacturing and mechanical performance of carbon fiber reinforced Polyamide-6 composites. Mater. Today Proc. 2022, 62, 6359–6363. [Google Scholar] [CrossRef]
  7. Pellegrini, A.; Palmieri, M.E.; Guerra, M.G. Evaluation of anisotropic mechanical behaviour of 316L parts realized by metal fused filament fabrication using digital image correlation. Int. J. Adv. Manuf. Technol. 2022, 120, 7951–7965. [Google Scholar] [CrossRef]
  8. Wu, H.; Fahy, W.; Kim, S.; Kim, H.; Zhao, N.; Pilato, L.; Kafi, A.; Bateman, S.; Koo, J. Recent developments in polymers/polymer nanocomposites for additive manufacturing. Prog. Mater. Sci. 2020, 111, 100638. [Google Scholar] [CrossRef]
  9. Al Rashid, A.; Khan, S.A.; Al-Ghamdi, S.G.; Koç, M. Additive manufacturing of polymer nanocomposites: Needs and challenges in materials, processes, and applications. J. Mater. Res. Technol. 2021, 14, 910–941. [Google Scholar] [CrossRef]
  10. Popescu, D.; Zapciu, A.; Amza, C.; Baciu, F.; Marinescu, R. FDM process parameters influence over the mechanical properties of polymer specimens: A review. Polym. Test. 2018, 69, 157–166. [Google Scholar] [CrossRef]
  11. Jiang, J.; Lou, J.; Hu, G. Effect of support on printed properties in fused deposition modelling processes. Virtual Phys. Prototyp. 2019, 14, 308–315. [Google Scholar] [CrossRef]
  12. Imran, R.; Al Rashid, A.; Koç, M. Review on computational modeling for the property, process, product and performance (PPPP) characteristics of additively manufactured porous magnesium implants. Bioprinting 2022, 28, e00236. [Google Scholar] [CrossRef]
  13. Al Rashid, A.; Qadir, S.A.; Koç, M. Microscopic analysis on dimensional capability of fused filament fabrication three-dimensional printing process. J. Elastomers Plast. 2021, 54, 385–403. [Google Scholar] [CrossRef]
  14. Dey, A.; Yodo, N. A Systematic Survey of FDM Process Parameter Optimization and Their Influence on Part Characteristics. J. Manuf. Mater. Process. 2019, 3, 64. [Google Scholar] [CrossRef]
  15. Vishwas, M.; Basavaraj, C.; Vinyas, M. Experimental Investigation using Taguchi Method to Optimize Process Parameters of Fused Deposition Modeling for ABS and Nylon Materials. Mater. Today Proc. 2018, 5, 7106–7114. [Google Scholar] [CrossRef]
  16. Al Rashid, A.; Koç, M. Fused Filament Fabrication Process: A Review of Numerical Simulation Techniques. Polymers 2021, 13, 3534. [Google Scholar] [CrossRef] [PubMed]
  17. Vanaei, H.R.; Khelladi, S.; Deligant, M.; Shirinbayan, M.; Tcharkhtchi, A. Numerical Prediction for Temperature Profile of Parts Manufactured using Fused Filament Fabrication. J. Manuf. Process. 2022, 76, 548–558. [Google Scholar] [CrossRef]
  18. Roy, M.; Wodo, O. Data-driven modeling of thermal history in additive manufacturing. Addit. Manuf. 2020, 32, 101017. [Google Scholar] [CrossRef]
  19. Lepoivre, A.; Boyard, N.; Levy, A.; Sobotka, V. Heat Transfer and Adhesion Study for the FFF Additive Manufacturing Process. Procedia Manuf. 2020, 47, 948–955. [Google Scholar] [CrossRef]
  20. D’Amico, T.; Peterson, A.M. Bead parameterization of desktop and room-scale material extrusion additive manufacturing: How print speed and thermal properties affect heat transfer. Addit. Manuf. 2020, 34, 101239. [Google Scholar] [CrossRef]
  21. Wang, L.; Chen, X.; Kang, S.; Deng, X.; Jin, R. Meta-modeling of high-fidelity FEA simulation for efficient product and process design in additive manufacturing. Addit. Manuf. 2020, 35, 101211. [Google Scholar] [CrossRef]
  22. Costa, S.; Duarte, F.; Covas, J. Thermal conditions affecting heat transfer in FDM/FFE: A contribution towards the numerical modelling of the process. Virtual Phys. Prototyp. 2015, 10, 35–46. [Google Scholar] [CrossRef]
  23. Comminal, R.; Serdeczny, M.P.; Pedersen, D.B.; Spangenberg, J. Numerical modeling of the strand deposition flow in extrusion-based additive manufacturing. Addit. Manuf. 2018, 20, 68–76. [Google Scholar] [CrossRef]
  24. Serdeczny, M.P.; Comminal, R.; Mollah, M.T.; Pedersen, D.B.; Spangenberg, J. Numerical modeling of the polymer flow through the hot-end in filament-based material extrusion additive manufacturing. Addit. Manuf. 2020, 36, 101454. [Google Scholar] [CrossRef]
  25. Costa, S.; Duarte, F.; Covas, J. Estimation of filament temperature and adhesion development in fused deposition techniques. J. Mater. Process. Technol. 2017, 245, 167–179. [Google Scholar] [CrossRef]
  26. Phan, D.D.; Horner, J.S.; Swain, Z.R.; Beris, A.N.; Mackay, M.E. Computational fluid dynamics simulation of the melting process in the fused filament fabrication additive manufacturing technique. Addit. Manuf. 2020, 33, 101161. [Google Scholar] [CrossRef]
  27. D’Amico, A.; Peterson, A.M. An adaptable FEA simulation of material extrusion additive manufacturing heat transfer in 3D. Addit. Manuf. 2018, 21, 422–430. [Google Scholar] [CrossRef]
  28. ASTM (D638); American Society for Testing and Materials. Standard Test Method for Tensile Properties of Plastics (D 638-02a)—SCAN VERSION. ASTM: West Conshohocken, PA, USA, 2006; Volume 8, pp. 46–58. [CrossRef]
  29. Palmero, E.M.; Casaleiz, D.; de Vicente, J.; Hernández-Vicen, J.; López-Vidal, S.; Ramiro, E.; Bollero, A. Composites based on metallic particles and tuned filling factor for 3D-printing by Fused Deposition Modeling. Compos. Part A Appl. Sci. Manuf. 2019, 124, 105497. [Google Scholar] [CrossRef]
  30. Kulich, D.M.; Gaggar, S.K.; Lowry, V.; Stepien, R. Acrylonitrile–Butadiene–Styrene (ABS) Polymers. In Kirk-Othmer Encyclopedia of Chemical Technology; American Cancer Society: Atlanta, GA, USA, 2003. [Google Scholar] [CrossRef]
  31. Cress, A.K.; Huynh, J.; Anderson, E.H.; O’Neill, R.; Schneider, Y.; Keleş, Ö. Effect of recycling on the mechanical behavior and structure of additively manufactured acrylonitrile butadiene styrene (ABS). J. Clean. Prod. 2021, 279, 123689. [Google Scholar] [CrossRef]
  32. Shaikh, M.Q.; Singh, P.; Kate, K.H.; Freese, M.; Atre, S.V. Finite Element-Based Simulation of Metal Fused Filament Fabrication Process: Distortion Prediction and Experimental Verification. J. Mater. Eng. Perform. 2021, 30, 5135–5149. [Google Scholar] [CrossRef]
  33. Cattenone, A.; Morganti, S.; Alaimo, G.; Auricchio, F. Finite Element Analysis of Additive Manufacturing Based on Fused Deposition Modeling: Distortions Prediction and Comparison With Experimental Data. J. Manuf. Sci. Eng. 2018, 141, 011010. [Google Scholar] [CrossRef]
  34. Al Rashid, A.; Koҫ, M. Creep and Recovery Behavior of Continuous Fiber-Reinforced 3DP Composites. Polymers 2021, 13, 1644. [Google Scholar] [CrossRef]
  35. Zhou, Y.; Nyberg, T.; Xiong, G.; Liu, D. Temperature Analysis in the Fused Deposition Modeling Process. In Proceedings of the 2016 3rd International Conference on Information Science and Control Engineering (ICISCE), Beijing, China, 8–10 July 2016; pp. 678–682. [Google Scholar] [CrossRef]
  36. De Macedo, R.Q.; Ferreira, R.T.L.; Jayachandran, K. Determination of mechanical properties of FFF 3D printed material by assessing void volume fraction, cooling rate and residual thermal stresses. Rapid Prototyp. J. 2019, 25, 1661–1683. [Google Scholar] [CrossRef]
  37. Sun, Q.; Rizvi, G.; Bellehumeur, C.; Gu, P. Effect of processing conditions on the bonding quality of FDM polymer filaments. Rapid Prototyp. J. 2008, 14, 72–80. [Google Scholar] [CrossRef]
  38. Bellehumeur, C.; Li, L.; Sun, Q.; Gu, P. Modeling of Bond Formation Between Polymer Filaments in the Fused Deposition Modeling Process. J. Manuf. Process. 2004, 6, 170–178. [Google Scholar] [CrossRef]
Polymers 14 03482 g001 550
Figure 1. CAD Model for ASTM D638 type I specimen with dimensions.
Figure 1. CAD Model for ASTM D638 type I specimen with dimensions.
Polymers 14 03482 g001
Polymers 14 03482 g002 550
Figure 2. Maximum displacements in numerical model simulations at different extrusion temperatures: (a) 210 °C, (b) 220 °C, (c) 230 °C, (d) 240 °C, (e) 250 °C, (f) 260 °C.
Figure 2. Maximum displacements in numerical model simulations at different extrusion temperatures: (a) 210 °C, (b) 220 °C, (c) 230 °C, (d) 240 °C, (e) 250 °C, (f) 260 °C.
Polymers 14 03482 g002
Polymers 14 03482 g003 550
Figure 3. Specimen dimensions for numerically predicted and fabricated samples 3D printed at different temperatures: (a) 210 °C, (b) 220 °C, (c) 230 °C, (d) 240 °C, (e) 250 °C, (f) 260 °C.
Figure 3. Specimen dimensions for numerically predicted and fabricated samples 3D printed at different temperatures: (a) 210 °C, (b) 220 °C, (c) 230 °C, (d) 240 °C, (e) 250 °C, (f) 260 °C.
Polymers 14 03482 g003
Polymers 14 03482 g004 550
Figure 4. Numerically predicted (above) and physically measured (below) warpage for 3D printed samples at different temperatures: (a) 210 °C, (b) 220 °C, (c) 230 °C, (d) 240 °C, (e) 250 °C, (f) 260 °C.
Figure 4. Numerically predicted (above) and physically measured (below) warpage for 3D printed samples at different temperatures: (a) 210 °C, (b) 220 °C, (c) 230 °C, (d) 240 °C, (e) 250 °C, (f) 260 °C.
Polymers 14 03482 g004
Polymers 14 03482 g005 550
Figure 5. Residual stresses in numerical model simulations at different extrusion temperatures: (a) 210 °C, (b) 220 °C, (c) 230 °C, (d) 240 °C, (e) 250 °C, (f) 260 °C.
Figure 5. Residual stresses in numerical model simulations at different extrusion temperatures: (a) 210 °C, (b) 220 °C, (c) 230 °C, (d) 240 °C, (e) 250 °C, (f) 260 °C.
Polymers 14 03482 g005
Polymers 14 03482 g006 550
Figure 6. Specimen dimensions for numerically predicted 3D printed at different layer heights: (a) 0.15 mm, (b) 0.10 mm.
Figure 6. Specimen dimensions for numerically predicted 3D printed at different layer heights: (a) 0.15 mm, (b) 0.10 mm.
Polymers 14 03482 g006
Polymers 14 03482 g007 550
Figure 7. Numerically predicted warpage for 3D printed samples at different layer heights: (a) 0.15 mm, (b) 0.10 mm.
Figure 7. Numerically predicted warpage for 3D printed samples at different layer heights: (a) 0.15 mm, (b) 0.10 mm.
Polymers 14 03482 g007
Table
Table 1. Effect of extrusion temperature on dimensional variations in numerical model and 3DP parts.
Table 1. Effect of extrusion temperature on dimensional variations in numerical model and 3DP parts.
RunExtrusion Temperature (°C)ModelLength
(mm)		Grip Width
(mm)		Gauge Width
(mm)		Thickness
(mm)		Warpage
--Design165.0019.0013.004.00-
1210Numerical162.24 ± 0.0118.75 ± 0.0212.85 ± 0.023.89 ± 0.020.13 ± 0.01
3DP163.76 ± 0.0118.89 ± 0.0112.85 ± 0.013.90 ± 0.020.25 ± 0.01
2220Numerical162.10 ± 0.0118.74 ± 0.0212.85 ± 0.013.89 ± 0.020.14 ± 0.01
3DP164.21 ± 0.0118.91 ± 0.0212.88 ± 0.023.93 ± 0.020.23 ± 0.01
3230Numerical161.96 ± 0.0118.73 ± 0.0212.84 ± 0.013.88 ± 0.020.15 ± 0.01
3DP164.89 ± 0.0118.90 ± 0.0212.91 ± 0.023.95 ± 0.020.20 ± 0.01
4240Numerical161.81 ± 0.0218.71 ± 0.0212.83 ± 0.023.87 ± 0.020.15 ± 0.01
3DP164.96 ± 0.0118.95 ± 0.0212.95 ± 0.023.98 ± 0.020.19 ± 0.01
5250Numerical161.67 ± 0.0218.69 ± 0.0212.83 ± 0.013.87 ± 0.020.16 ± 0.01
3DP164.90 ± 0.0118.93 ± 0.0212.92 ± 0.023.97 ± 0.020.21 ± 0.01
6260Numerical161.53 ± 0.0218.68 ± 0.0212.82 ± 0.013.86 ± 0.020.17 ± 0.01
3DP163.85 ± 0.0118.88 ± 0.0212.90 ± 0.013.94 ± 0.020.24 ± 0.01
Table
Table 2. Effect of Layer Height on Dimensional Variations in Numerical Model and 3DP Parts.
Table 2. Effect of Layer Height on Dimensional Variations in Numerical Model and 3DP Parts.
RunLayer Height
(mm)		ModelLength
(mm)		Grip Width
(mm)		Gauge Width
(mm)		Thickness
(mm)		Warpage
--Design165.0019.0013.004.00-
10.20Numerical161.81 ± 0.0218.71 ± 0.0212.83 ± 0.023.87 ± 0.020.15 ± 0.01
20.15Numerical162.26 ± 0.0118.77 ± 0.0112.88 ± 0.013.91 ± 0.020.13 ± 0.01
30.10Numerical165.00 ± 0.0118.95 ± 0.0112.92 ± 0.013.96 ± 0.010.08 ± 0.01
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Al Rashid, A.; Koç, M. Experimental Validation of Numerical Model for Thermomechanical Performance of Material Extrusion Additive Manufacturing Process: Effect of Process Parameters. Polymers 2022, 14, 3482. https://doi.org/10.3390/polym14173482

AMA Style

Al Rashid A, Koç M. Experimental Validation of Numerical Model for Thermomechanical Performance of Material Extrusion Additive Manufacturing Process: Effect of Process Parameters. Polymers. 2022; 14(17):3482. https://doi.org/10.3390/polym14173482

Chicago/Turabian Style

Al Rashid, Ans, and Muammer Koç. 2022. "Experimental Validation of Numerical Model for Thermomechanical Performance of Material Extrusion Additive Manufacturing Process: Effect of Process Parameters" Polymers 14, no. 17: 3482. https://doi.org/10.3390/polym14173482

APA Style

Al Rashid, A., & Koç, M. (2022). Experimental Validation of Numerical Model for Thermomechanical Performance of Material Extrusion Additive Manufacturing Process: Effect of Process Parameters. Polymers, 14(17), 3482. https://doi.org/10.3390/polym14173482

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop