Next Article in Journal
Laser Heating Study of the High-Temperature Interactions in Nanograined Uranium Carbides
Previous Article in Journal
Processing and Thermal Conductivity of Bulk Nanocrystalline Aluminum Nitride
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Combined Experimental and Modeling Study for Pellet-Fed Extrusion-Based Additive Manufacturing to Evaluate the Impact of the Melting Efficiency

1
Centre for Polymer and Material Technologies (CPMT), Department of Materials, Textiles and Chemical Engineering, Ghent University, Technologiepark, 130, Zwijnaarde 9052, 9000 Ghent, Belgium
2
Centre for Textiles Science and Engineering (CTSE), Department of Materials, Textiles and Chemical Engineering, Ghent University, Technologiepark, 70A, Zwijnaarde 9052, 9000 Ghent, Belgium
3
Laboratory for Chemical Technology (LCT), Department of Materials, Textiles and Chemical Engineering, Ghent University, Technologiepark, 125, Zwijnaarde 9052, 9000 Ghent, Belgium
*
Author to whom correspondence should be addressed.
Materials 2021, 14(19), 5566; https://doi.org/10.3390/ma14195566
Submission received: 18 August 2021 / Revised: 8 September 2021 / Accepted: 21 September 2021 / Published: 25 September 2021

Abstract

:
To improve the product quality of polymeric parts realized through extrusion-based additive manufacturing (EAM) utilizing pellets, a good control of the melting is required. In the present work, we demonstrate the strength of a previously developed melt removal using a drag framework to support such improvement. This model, downscaled from conventional extrusion, is successfully validated for pellet-based EAM—hence, micro-extrusion—employing three material types with different measured rheological behavior, i.e., acrylonitrile-butadiene-styrene (ABS), polylactic acid (PLA) and styrene-ethylene-butylene-styrene polymer (SEBS). The model’s validation is made possible by conducting for the first time dedicated EAM screw-freezing experiments combined with appropriate image/data analysis and inputting rheological data. It is showcased that the (overall) processing temperature is crucial to enable similar melting efficiencies. The melting mechanism can vary with the material type. For ABS, an initially large contribution of viscous heat dissipation is observed, while for PLA and SEBS thermal conduction is always more relevant. It is highlighted based on scanning electron microscopy (SEM) analysis that upon properly tuning the finalization of the melting point within the envisaged melting zone, better final material properties are achieved. The model can be further used to find an optimal balance between processing time (e.g., by variation of the screw frequency) and material product performance (e.g., strength of the printed polymeric part).

1. Introduction

Over the past few years, the use of additive manufacturing (AM) has significantly increased in the production of polymeric products and material shapes [1,2,3,4,5,6,7,8,9,10]. Most polymeric AM applications start with a filament [9,10,11,12], which is forced through a heated nozzle, promoting its melting to then deposit consecutive material layers. This is performed according to geometric and process parameters that can have a strong influence on the properties of the final parts [13,14,15,16]. In recent years, however, a growing number of AM applications introduced the use of single-screw extruders (SSEs) fed with polymer pellets [17,18,19,20,21]. The global term for such AM modifications is extrusion-based AM (EAM), complementary with conventional large-scale extrusion [22,23,24,25,26].
EAM allows us to overcome several limitations of filament-based AM, which is also known as fused filament fabrication (FFF). Specifically, soft polymers can be manufactured that are difficult to process through traditional 3D printing [19,20]. Printing employing polymer pellets leads to a lower cost of the final material [20] and, due to the higher possible melting rate [19], a faster processing speed for, in principle, a wider range of materials is realized [20]. It allows us to bypass the filament extrusion process that normally precedes the AM process, relieving the material of one extrusion cycle, leading to less degradation [27].
One of the remaining EAM challenges is to select, for a given polymeric material, the appropriate operating conditions to avoid/minimize both thermal degradation and solid slippage. Here, control of the melting zone is crucial. However, this is far from trivial, as the melting process depends on many interdependent parameters, ranging from process parameters such as the barrel temperature and screw frequency to the shear dependent viscosity and stress, and eventually to material properties such as density, thermal conductivity, and heat capacity [28,29,30,31,32,33,34,35,36,37]. Unfortunately, only a limited number of studies have investigated EAM melting, as (i) typical AM research is devoted to the macroscopic analysis of the printed product as such and (ii) EAM has only emerged as a key production technique this century. If melting is investigated, one needs to refer to studies on large scale extrusion, in which the dimensions are of a totally different order of magnitude [28,29,30,31,32,33,34,35,36,37].
Looking back at the long history of conventional extrusion, Maddock [28] was the first to experimentally investigate the melting mechanism in conventional SSE. Later on, Tadmor developed [29] and experimentally validated [30,31] an analytical model by considering a heat and mass balance for the melt film between the solid bed and the barrel. This model was then modified to include the effects of shear rate and temperature on the viscosity [24]. Upon accounting for shear thinning with, for instance, a power law index and assuming melt removal by drag, the number of screw rotations or, equivalently, the unwounded screw length at which melting stops for a given set of processing conditions could be more reliably predicted.
In parallel with the work and modifications of Tadmor, other large-scale extrusion melting models have been developed and applied. Donovan [32], for example, focused on the melting efficiency, assuming a constant melt film thickness and an identical down-channel velocity for both the solid bed and the melt pool at any cross section. This model allowed, at first glance, a good prediction of the melting behavior, but was later criticized for the poor prediction of the overall melting profiles [33]. Edmondson and Fenner [34] introduced a model that accounts for the presence of a molten film between the screw and the solid bed, which is thus an extra drag contribution compared to the bed and barrel. The predictions of the model appeared to be satisfying for both the solid bed and pressure profiles. Furthermore, Mount et al. [35] proposed and satisfyingly validated analytical equations for the melting rate prediction that did not require an iterative calculation. This was done by reducing part of the constitutive equations to non-dimensional forms, employing characteristic scales and non-dimensional numbers.
Syrjala [36] in turn evaluated the melting profile for SSEs by only using relevant flow field equations, without any focus on the mechanism of melting. Notably, in a more recent work, Altınkaynak [37] performed three-dimensional finite element simulations of the melting process by solving the conservation equation for mass, momentum, and energy along with a generalized Newtonian constitutive equation. Such model efforts thus go beyond analytical approximate solutions, but the solving of the related equations is too slow in the context of online process control.
In our recent work [38], we theoretically studied the efficiency of EAM melting by downsizing Tadmor’s pioneering melting mechanism from conventional large-scale extrusion to micro-extrusion, selecting acrylonitrile–butadiene–styrene (ABS) as a viscoelastic reference polymer. We highlighted that the variation of the Brinkman (Br) number [39,40,41,42,43], which is the ratio between the melting contribution through heat dissipation and that through conduction, along the unwounded screw length is different for EAM compared to conventional SSE. This Br difference is due to a different balance of shear rates, melt layer thicknesses and screw geometry, consistent with speculations in previous studies [44,45,46]. In the process of scaling down SSE, a reduction in the viscous dissipation was thus witnessed in our previous work. This could be made clearer upon comparing the melting profiles for micro-extrusion and conventional extrusion under typical operation conditions for each technique on its own [38], consistent with preliminary results [47]. We also put forward in our previous theoretical study that screw parameters such as the screw frequency and pitch angle and material parameters such as the power-law index influence the relative position of the final EAM melting point.
One can thus not generalize the melting behavior of one polymer directly to another, highlighting the relevance of the development of more generic theoretical frameworks to support the identification of ideal EAM conditions. Furthermore, it can be expected that good printed macroscopic properties require a good melting. Consistent with this, the relevance of an adequate melting was also presented in previous studies for PLA [48] and ABS [49], showing how low printing temperatures can cause a steep decrement in the properties of the final parts. Hence, an understanding of the melting can help to better identify better printing conditions.
In the present work, we therefore further increase the impact of our modeling framework by benchmarking the EAM melt model to in-house determined experimental data. To record these data, we perform a series of so-called screw-freeze experiments [28] to visually inspect the positioning of the solid bed and melt pool at distinct screw positions. The material portfolio is also expanded beyond simply ABS to include several rheological behaviors. The first polymeric material ABS is characterized by an amorphous structure, the second polymeric material is polylactic acid (PLA) which has, in contrast, a semicrystalline structure, and the third polymeric material is styrene-ethylene-butylene-styrene polymer (SEBS), which falls into the class of thermoplastic elastomers. A comparison is included regarding the melting mechanism for the three different materials, through a comparison of the Brinkman number variation [39,40,41,42,43]. Moreover, the quality of the printed parts is related to the melting efficiency, as accessible for any set of operating conditions through application of the modeling framework. In addition, a first step toward a full macroscopic property analysis is included by including imaging and tensile results.

2. Materials and Micro-Extruder Machine Details

In this study, three commercially available extrudable polymeric materials were used: (i) acrylonitrile butadiene styrene polymer (ABS) MagnumTM M3404 from Dow Chemical Company (Midland, MI, USA), (ii) polylactic acid (PLA) Ingeo ™ 3D850 from NatureWorks (Minnetonka, MN, USA), and (iii) styrene-ethylene-butylene-styrene polymer (SEBS) KRATONTM G1657M from Kraton (Huston, TX, USA). Table 1 summarizes the most important materials characteristics. Additionally, a small amount of Clariant Plastiflac polyethylene (PE)-based red pigment was used for the screw-freezing experiments.
The micro-extrusion machine geometry is the same for all of the material processing and simulations, and its details are reported in Table 2. The dimensions of the geometry are consistent with our previous work [38] and resemble a typical micro-extruder design.

3. Experimental and Theoretical Methods

3.1. Melting Model

The melting profile was simulated using the melting model and the micro-extruder geometry described in our previous theoretical study [38]. For completeness, the main model parameters are reported in Table 3, making a differentiation between the three polymer types and extracting values from the previous Tables (Table 1, Table 2, Table 3, Table 4, Table 5 and Table 6). Essential in the model is the calculation of the melting rate per unit of channel width q:
q = X   ρ m     | V bx |   λ + c s   ( T m T 0 ) + c m   Θ   ( T b T m ) ( T b T m ) k ( 1 + B r 2 )
In this equation, reported from our previous work [38], X is the width of the solid bed (initially X is equal to the channel width W),   ρ m   is the density of the molten phase, Vbx is the velocity of the molten phase at the barrel in the x direction (as reported in Figure S1 of the Supplementary Materials), λ is the heat of fusion, cm is the specific heat capacity of the molten phase, cs is the specific heat capacity of the solid phase, Θ is the average temperature of the melt that is related to the Brinkman number Br [39,40,41,42,43], Tm is the melting temperature, k is the thermal conductivity, and Tb is the barrel temperature.
Br is a dimensionless number that indicates the ratio between the viscous heating and the (molecular) conduction heating [39,40,41,42,43]. The viscous heating is generated by the friction between the flowing elements at different velocities and is finally supplied by the engine that turns the extrusion screw. The conduction heat, on the other hand, is supplied by the external heater bands. Br is defined in Equation (2), reported from our previous work [38], and was used to evaluate differences in the melting mechanism between the materials tested.
Br = η   V j n + 1 δ ¯ n 1   k ( T b T m )  
In Equation (2), η is the shear rate dependent melt viscosity, Vj is the velocity difference between the velocity of the solid bed in the z direction Vsz and the velocity of the melt layer at the barrel surface Vb, n is the pseudo-plasticity index, and δ ¯ is the average thickness of the molten layer that can be assessed and updated along the unwounded screw as explained in our previous work [38]. Note that Vbx, reported from our previous work [38], is a function of the screw frequency N, the inside diameter of the extruder barrel Db, and the helix pitch angle of the screw Ɵ:
V bx = V b   s i n Ɵ = π   D b   N
It should be stressed that Vsz is dependent on the volumetric flow rate at the inlet of the extruder, Q0. The value of Q0 was estimated by assuming stationary conditions in the extruder and weighing the extruded amount for defined intervals of time. Hence, no ball-park value was selected, but for each polymeric material, a true value based on actual steady-state experiments was considered.

3.2. Rheological Measurements and Cross Model Fitting

Small amplitude oscillatory shear (SAOS) tests were performed using a MCR 702 rheometer (Anton Paar, Graz, Austria). The storage modulus (G′) and loss modulus (G′′) were monitored as a function of frequency (1 × 10−1 to 6 × 102 rad s−1), using the parallel plate configuration with a 25 mm diameter, a gap of 1 mm and an amplitude of 1%. The tests were performed under nitrogen atmosphere at the same (target or average) processing temperatures used for the actual EAM experiments utilizing the operating settings as reported in Table 3. The latter table provides an overall summary of parameters directly relevant for the results and discussion. Specifically for ABS, the processing temperature is higher, as preliminary screening has shown that the viscosity is higher. The disks were obtained by compression molding of dried pellets at the same temperatures, considering a diameter of 25 mm and 1 mm thickness. The discs were also dried before the rheological measurements.
As shown by the symbols in Figure 1, the three polymeric materials, with specifications in Table 1, globally have a similar rheological behavior, with the viscosity η exhibiting a shear thinning behavior after the Newtonian region, as is most evident in the case of PLA, shown in Figure 1b. There are, however, keeping in mind the different temperatures at which the curves were obtained, some clear differences in the values for (η0), which is the Newtonian viscosity limit before the critical shear rate (γcr). It particularly follows that ABS is considerably more viscous than SEBS, with PLA exhibiting the lowest (η0) value among the three polymeric materials. A previous study [58] observed that the rheological behavior of ABS is strongly dependent on the (average) molar mass of the SAN (styrene-acrylonitrile) phase. In particular, an increase in the mass average molar mass of the SAN phase leads to the increase in viscosities under both shear and uniaxial extension. This increment in average molar mass leads to longer SAN chains being grafted onto the PB (poly-butadiene) particles. This results in the formation of a denser structural “network”, with a higher yield stress and storage modulus in the low frequency plateau region. For completeness, it is mentioned here that strong shear thinning effects are also indicative of a wide molar mass distribution for the SAN phase [59]. Furthermore, for SEBS, a more complex rheological behavior is recorded in Figure 1c. In a previous work [60], it has been observed that the shear thinning effect in the region between 0.1 and 1.0 rad s−1 is only present at temperatures at which the SEBS’ microphase-separated state changes into a single-phase state, which is called the order-disorder transition (ODT). This effect is a strong indication of the presence of several melt structures, and the extreme shear thinning behavior at low shear rates can be associated with a melt yield stress.
The experimental rheological data for each polymer were also used to formally fit a Cross flow model at the selected temperature:
η = η 0 ^ 1 + ( η 0   γ ˙ τ * ) 1 n
η 0 ^ = η 0 ( 1 + ( η 0   γ cr   ˙ τ * ) 1 n )
In Equations (4) and (5), reported from our previous work [38], η is the shear rate dependent melt viscosity, n is the pseudo-plasticity index (or shear thinning factor), τ* is the critical Cross model stress at the transition to shear thinning, and γ is the shear rate (formally assumed equal to the frequency), while η 0 ^ is a short notation to retrieve the Newtonian limit thus before the critical shear rate γ cr   ˙ .
The corresponding fits are highlighted as lines in Figure 1a–c, and the Cross-model parameters are reported in Table 3. It follows from Figure 1 that for ABS and PLA, good descriptions are obtained with the Cross model, but for SEBS, the curve’s fit is less accurate. A cross-model only grasps an average overall behavior or one structure, and thus not the ODT.

3.3. Thermal Properties

The thermal conductivity k of the polymeric samples was measured using the Transient Plane Source method with a Hot Disk TPS 2500S (Hot Disk, Göteborg, Sweden) according to ISO 22007-2. The conditions of the experiments are reported in Table 4. The k value was obtained by an average of 3 measurements and is assumed to be representative (during the processing) at all temperatures by taking the average over all data in Table 4.
The specific heat capacities cs (solid) and cm (melt), the heat of fusion λ, and the “melting” temperature Tm of the polymeric materials were evaluated through differential scanning calorimetry (DSC), using a Polyma DSC 214 from Netzsch. The measurements were obtained using the ASTM standard E1269, which requires three scans: (i) a baseline scan, (ii) a scan using a sapphire standard, and (iii) the sample scan. The specific heat capacity was measured at 300 K, which was the initial temperature of the three materials, and at 460 K, at which all materials were in the molten state. The experiments were conducted at a heating rate of 5 K min−1. The conditions and results are reported in Table 5. Again, constant (average) values were considered in what follows, further supported by the sensitivity analysis in our previous work [38].
To take into account a possible difference in the PLA crystallinity (special case in Table 1 regarding crystallinity) between the injection molded samples for the Hot Disk and the solid pellets for EAM, further DSC tests were conducted at 10 K min−1. A variation of only 1.70% was found for the area of the crystalline peak, as shown in Figure S2 in the Supporting Information. The effect of this variation in the PLA crystallinity during the sample preparation was neglected, as previous studies observed a marginal impact of the degree of crystallinity on the thermal conductivity [61]. Moreover, previous studies on the melting mechanisms in SSEs observed that small variations in the thermal conductivity have a neglectable impact on the solid bed profile (SBP) [21,38].

3.4. Solid and Melt Density

The density of the solid phase was measured according to the Immersion Method of ISO/DIS 1183-1, weighing the pellets in air and in ethanol as an auxiliary liquid. The weight scale used was a Precisa XR205SM-DR manufactured by Precisa Gravimetrics AG (Dietikon, Switzerland).
A Melt Indexer MPX 62.92 by GÖTTFERT ® (Buchen, Germany) was used to measure the Melt Volume Rate (MVR) of the different polymeric materials under the conditions described in Table 6, according to ISO 1133:2005 and ISO 1133:2005 Cor.1:2006. Under the same conditions, the mass flow rate (MFR) was estimated by weighing the output from the melt indexer for a given time period. The melt density was evaluated as the ratio between MFR and MVR. The results are also presented in Table 6. The final density values were then averaged and used at any temperature for the processing simulations.

3.5. Screw-Freezing Experiments

The EAM melting model predictions were validated through a series of screw-freezing experiments that were conducted through rapid water cooling of the operating micro-extruder (Ghent University, Gent, Belgium). This was performed through water-cooled aluminum bands that were positioned along the extruder barrel, as visible in Figure 2a. After a screw-freezing experiment was completed, the extruder was disassembled. Then, the solidified polymer channel was extracted from the barrel and cut into sections to validate the model predictions on SBP.
The screw freezing technique was first introduced by Maddock for conventional extruders [23]. This technique starts with the extruder running in steady state. The screw rotation is then stopped, and the barrel is air cooled to solidify the polymer in the screw channel. The screw with the solidified polymer is extracted and the cross-sections of the solidified polymer are polished and investigated. In the present work on micro-extruders, the same approach was applied, although some modifications were made. In greater detail to highlight the regions where the polymer was molten, a small amount of PE pigmented feedstock polymer was introduced at a mass ratio of 0.5:100. The solid bed still showed the natural color of the pellets, while the regions where the polymer had melted were strongly colored by the pigment, according to the same principle described in Maddock’s study [28]. Additional water-cooling rings were added to the compression and metering sections of the extruder barrel to improve the cooling behavior. The cooling rings were assembled to promote the cooling of the lower section first and then the upper section. The cooling curves for the two sections after EAM competition are depicted in Figure 2b, with the positions of the thermocouples visible in Figure 2d. It follows that a relatively efficient cooling is realized.
For most of the tests with ABS, the solidified polymer channel had a negligible deformation during the extraction process, which allowed for epoxy embedding to better showcase the melting mechanism and the SBP. The embedded samples were finally polished, and a high-resolution scanner was used for the image acquisition, with the results depicted in Figure 2 for three selected screw frequencies. For the tests with PLA, SEBS, and two specific ABS tests, due to the deformations or fractures of the channel during the extraction process, it was impossible to replicate the methodology used for the previous samples mentioned. After cutting and polishing of the sections, the image was acquired through a camera. Examples of sections for the non-embedded samples are visible in Figure S3 of the Supporting Information.
Similarly to what was observed by Maddock and Tadmor [28,29,30,31] for conventional extrusion, the polymer melt appears at first sight as a film in between the hot barrel and the compacted pellets, as conceptually shown in Figure 3a. Advancing along the channel, however, the recorded behavior departs from what Maddock and Tadmor observed [28,29,30,31], due to the clear presence of an outer recirculation zone, as conceptually depicted in Figure 3b. This can be explained by the different geometry of the micro-extruder, in particular the conical shape of the screw and barrel that leads to a substantial variation in the velocity profiles, as pointed out in our previous work [38].
Due to the presence of such outer recirculating flow, the width of the molten phase was not constant along the solid bed (in a given screw element), and the percentage of unmolten area over the total area was selected to represent the ratio of the solid bed width and screw width for comparison with the numerical model predictions. The area (A) analysis was executed using the IMAGEJ image processing and analysis software (Version 1.8.0, U. S. National Institutes of Health, Bethesda, MD, USA) for the different sections. For the non-embedded samples, a different image acquisition interpretation was found to be related to the separation of the two phases, due to the lack of a very clear separation in the color distributions in contrast to the very clear color contrast for the embedded samples. The threshold for these cases was selected based on visual feedback, in coherence with the melting mechanism. Too high or low values of the threshold would, for example, show the presence of molten or solid particles in positions not consistent with the expected melting mechanism.

3.6. SEM Analysis of Extruded and Solidified Samples

A Phenom Pro X Desktop Scanning Electron Microscope (SEM) from ThermoFisher Scientific (Waltham, MA, USA) was used to investigate the surfaces of extruded filament samples at different screw frequencies. The images were acquired using a backscatter electron detector (BSD) technique with a beam energy of 15 keV with a magnification of 490:1.
The scope of the analysis was to investigate the presence of surface defects caused by prolonged residence times, or by the presence of solid particles in the metering section. The selected frequencies were 0.5, 5 and 20 rpm, based on model screening, allowing us to identify the position of the point of melt finalization. Two filament samples were always collected from the extrusion under steady state conditions, at a 15 min interval. The process was repeated using 2 nozzles with 0.8 and 1.5 mm final opening diameter, for a total of 12 samples.

3.7. Tensile Testing of the Extruded Filaments

An Instron 5566 tensile machine (Instron, Norwood, MA, USA) was used to investigate the tensile properties of filaments of ABS extruded at different screw frequencies. The method used for the tensile test was based on a previous work [62]. A constant velocity of 10 mm·min−1 was used to evaluate the stress–strain relationship. The gauge length was set at 50 mm. Again, the selected screw frequencies were 0.5, 5 and 20 rpm.

4. Results and Discussion

4.1. Model Validation Regarding the Point of Melt Finalization for Three Polymeric Materials

Three screw-freezing experiments were first performed for the ABS polymer at a nominal screw speed of 5 rpm, and the results after data analysis are reported in Figure 4a (symbols with connections). The average (mean) values of these experimental results are compared in Figure 4b with the model predictions, inputting the rheological data from Figure 1a (parameters in Table 3). For completeness, these modeled data are also reported in Figure 4a. It follows from Figure 4b that the melting is performed in the compression section at ca. 250 mm and the decrease in the solid fraction is gradual and sufficiently captured by the model. The details of the image processing steps are provided in Figures S4–S6 of the Supplementary Materials, with, in the case of the first 2 Supplementary Materials figures related to experiment 1 and 2, non-embedded samples being considered, and for the last figure (Figure S6 of the Supplementary Materials), focus is on the embedded sample (experiment 3).
For ABS, two additional experiments were performed at 2 and 8 rpm to evaluate the effect of an increment and a decrement in the screw speed. The results are depicted in Figure 5a,c directly as average experimental results, also including the model predictions. For the sake of comparison, the ABS results from before at 5 rpm are also repeated in Figure 5b. The details of the image processing are visible in Figure S7 of the Supplementary Materials for the 2 rpm test and in Figure S8 of the Supplementary Materials for the 8 rpm test. In Figure 5, a good agreement is observed between the experimental results and the numerical model prediction for all three screw frequencies, also realizing that a model can wipe out experimental issues, in the present work due to non-embedding. Closer inspection reveals that the observed final differences between the 2 and 5 rpm results are still safely within the limits of the envisaged melting zone (but with 50 mm difference), while for 8 rpm, a shift in the melting profile to higher values is observed, even up to 300 mm, which is outside the envisaged melting zone.
A further interpretation of Figure 5a–c also enables a better assessment of the experimental errors and the relevance of the melting model. For example, for 2 rpm, the observed difference between the position of the point of melt finalization in the numerical model and the experiments is around 17.9% for the ninth experimental point. Focusing on the eighth experimental point at which the image processing delivers an estimate of 98.89% of molten material, the difference drops to about 8.2%. For 5 rpm, in turn, the observed difference is around 9.2% for the 9th experimental point. Considering the eighth experimental point, at which the molten fraction is estimated as 97.5%, the difference now increases to 18.8%. For 8 rpm, the observed variation is around 1.5% for the 10th experimental point, with the next point exhibiting about 96.4% of solid material and a 0.8% distance from the experimental profile. Hence, one can indeed state that the model filters out natural experimental errors and is a more reliable tool to fully extrapolate the final point of melting.
As shown in Figure S9a of the Supplementary Materials, the accuracy of the simulations’ results depends on the correctness of the input parameters. Employing the Q0 values from Table 3, the modeling outcome is in good agreement with the experimental results (full lines). Instead, the use of ballpark Q0 values, as considered in previous theoretical work [38], is too approximate (dashed lines). The need of benchmarked input parameters is further highlighted in Figure S11 of the Supplementary Materials, showing that for the ABS case, the use of the PLA and SEBS Q0 and n values from Table 3 are incorrect. Furthermore, in a previous study on SSEs [63], Mount and Chung also highlighted the importance of the method used to evaluate the solid bed velocity.
In view of extra model validation, a single screw freezing experiment has been additionally performed for the PLA and SEBS polymeric materials at the nominal screw speed of 5 rpm. The experimental (symbols) and modeling (lines) data are depicted in Figure 6a,b, again displaying the average experimental results. The corresponding results from the image processing are now visible in Figures S9 and S10 of the Supplementary Materials. It is worthwhile mentioning that while the color distribution for ABS and PLA clearly shows in most cases a bimodal distribution, allowing for a clear separation between the different phases, for SEBS, the threshold for separation between the two phases could only be obtained based on visual feedback. As shown in the Supplementary Materials, the SEBS distribution is not clearly bimodal, and thus this somewhat more arbitrary approach for melt efficiency determination is the most recommended.
It follows from Figure 6a,b that the numerical model predictions are again in good agreement with the experimental data, showing a good overall reliability, taking into account the different morphology of the materials in Table 1 and a wide range of viscosities, as expressed in Figure 1 and Table 3. For PLA, a difference in the position of the experimental and simulated point of melt finalization of about 12.3% is observed, while for SEBS, this difference is about 6.9%. Hence, the model is capable of tackling different material morphologies, highlighting the model’s efficacy in preventing experimental errors. Closer inspection shows that the melting for PLA is slightly slower and takes place for longer than the envisaged melting zone, while SEBS melting is still within the envisaged melting zone and is performed at a similar position as for ABS (Figure 5b). Overall, it can thus be concluded that the processing temperature variation in Table 3 allowed us to create similar global behavior in melting. Hence, this operating condition is crucial in the overall design.
A detailed model-driven approach is thus within reach if more rheological data and Arrhenius parameters regarding rheology are available, facilitating automated model screening. In this context, the potential of the model is further highlighted in Figure 6c by comparing at different rpm the melting behavior for the three polymeric materials. It follows that the relative trend (melting faster for the order PLA, SEBS and ABS) is always respected but a different rpm threshold can exist for each material type to obey the melting zone idea of the original micro-extruder design.

4.2. The Impact of the Polymer Material Type on the Melting Mechanism

Combining the modeled data in Figure 5 and Figure 6, the melting behavior can also be studied at a fixed rpm, varying the material type. This is conducted in the present subsection, specifically at 5 rpm. While the position of the point of melting finalization at 5 rpm moves, e.g., by 8% between ABS and PLA, as shown in Figure 7b, the underlying melting mechanism, as included in Figure 7a, shows a very different behavior. In this subplot, focus is on the variation of Br, with lower values of this number indicative of a higher share of the thermal energy being supplied by the external heater bands compared to the viscous heat generated by friction.
It follows from Figure 7a that the Br values show a very large difference in the heating mechanism for the different materials. While for the ABS material the viscous heat generated is initially still higher than the conduction heat (Br > 1), for the SEBS material, the viscous heat always represents a small fraction, and for the PLA material, it appears to be negligible. The reason for the difference in the melting mechanism can be explained by the differences in the viscosity profiles highlighted in Figure 1, and again highlights the relevance of the rheological curves.

4.3. Impact of Melting Efficiency on 3D Printed Properties

To evaluate the impact that bad melting can have on the final 3D printing applications, the surface of the extruded filament, selecting first a nozzle diameter of 0.8 mm, was analyzed with SEM. For a meaningful comparison, sufficiently different frequencies were selected so that the point of melt finalization is located early in the compression section (0.5 rpm case), at the end of the compression section (5 rpm case; reference case from before), and at the end of the metering section (20 rpm case). This is further confirmed by the model lines in Figure 8. For the 20 rpm case, in particular, the solid bed at the end of the metering section still accounts for about 0.9% of the section area. The SEM images, always recorded based on data in replicate format, are presented as subplots with the arrows highlighting the rpm value.
While the SEM images for both a screw frequency of 0.5 and 5 rpm show a similar surface roughness, the surface for the filament extruded at 20 rpm presents a higher number of defects and a higher roughness. As shown in Figures S15–S17 in the Supplementary Materials, the same differences in roughness are observed in case the nozzle diameter is equal to 1.5 mm instead of 0.8 mm. It is worthwhile mentioning that no substantial variation in Q0 was observed after the modification of the nozzle diameter. Hence, limitations exist regarding the melting efficiency to enable good material properties. In other words, the model can be used to identify suited operating conditions, saving experimental screening time and material usage.
Furthermore, the tensile curves in Figure S18 of the Supplementary Materials do not show substantial differences among the different frequencies. It is worthwhile mentioning that the amount of solid material predicted by the numerical model for the micro-extrusion at 20 rpm was under 1% and, due to the high temperature, the melting should have been completed in the filament after the extrusion. A slight increment in the yield strain and a decrement in the modulus with increasing frequency was observed, as shown in Figure 9a,b. These variations can be attributed to the variations in the average diameters of the filaments, caused by the different die swell, as visible in Figure 9c and coherent with previous extrusion work [64]. For higher frequencies, the larger diameters cause an increment in the cooling time, allowing the molecules to partially lose orientation, causing the drop in modulus. During the tensile test, the molecules can then be re-oriented, leading to a higher strain at yield for the larger diameters, and therefore higher frequencies.

5. Conclusions

Once rheological data are available, an EAM melting model, based on the drag removal principle and accounting for both conduction and heat dissipation, can be utilized to predict melting efficiencies. This has been confirmed by acceptably describing screw-freezing experiments for ABS under various screw frequencies and for two other polymeric materials, SEBS and PLA, for a fixed screw frequency. Notably, the model can wipe out natural experimental discrepancies, due to the impossibility to sometimes embed screw-freed samples. In any case, reliable model parameters on physicochemical parameters are needed, although sometimes well-measured averages suffice.
The EAM melting model can also be employed to perform the design, and here the material type plays a role. A sufficiently high processing temperature is needed to ensure sufficiently low viscosities. The model also allows us to establish the melting mechanism of the studied materials; for ABS only, a short period of a dominant contribution from heat dissipation was observed. The screw-freezing experiments further showcase that under several conditions, an outer circulation flow is present related to the geometry and dimensions of the micro-extruder.
The SEM analysis of the extruded filaments for a set of operating conditions that have a predicted point of melt finalization outside of the envisaged melting zone presents a higher amount of surface defects compared to the filaments obtained via operating conditions with a melting point well within the extruder’s limits. These defects can be attributed to the presence of unmolten material within the extruder, confirming the importance of the melting model in the process design of polymeric AM applications. Furthermore, the screw frequency also has an impact of the die swell and the tensile properties.
Hence, the current work highlights the relevance of controlling the melting in the overall process and the need for further material-oriented research to fully design the printing process.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/ma14195566/s1, Figure S1: Definition of x and z direction for the general equations for melt removal by drag, Figure S2: DSC for the PLA pellets and PLA injection molded hot disk sample. The melting temperature was defined as the onset temperature of the melting peak for pellets, Figure S3: Examples of sections of the solidified extruder channel in the micro extruder for a) PLA; b) SEBS, Figure S4: Image processing for ABS 5RPM (not embedded sample 1), Figure S5: Image processing for ABS 5RPM (not embedded sample 2), Figure S6: Image processing for ABS 5 RPM (embedded sample), Figure S7: Image processing for ABS 2RPM (embedded sample), Figure S8: Image processing for ABS 8RPM (embedded sample), Figure S9: Image processing for PLA 5 RPM (not embedded sample), Figure S10: Image processing for SEBS 5RPM (not embedded sample), Figure S11: (a) Relevance of the Q0 parameter; (b) Relevance of the n parameter, Figure S12: SEM Images for ABS extruded filament at 0.5 rpm, 0.8 mm nozzle, Figure S13: SEM Images for ABS extruded filament at 5 rpm, 0.8 mm nozzle, Figure S14: SEM Images for ABS extruded filament at 20 rpm, 0.8 mm nozzle, Figure S15: SEM Images for ABS at 0.5 rpm, 1.5 mm nozzle, Figure S16: SEM Images for ABS at 5 rpm, 1.5 mm nozzle, Figure S17: SEM Images for ABS at 20 rpm, 1.5 mm nozzle, Figure S18: Tensile curves for ABS filament extruded at a) 0.5 rpm; b) 5 rpm; c)20 rpm. Other operating conditions in Table 3 of the main document.

Author Contributions

Conceptualization, L.C. and D.R.D.; software, A.L.G., L.C. and D.R.D.; investigation, A.L.G., R.F., D.V.A.C., M.E., L.C. and D.R.D.; writing—original draft preparation A.L.G.; writing—review and editing, A.L.G., R.F., D.V.A.C., M.E., L.C. and D.R.D.; supervision, L.C. and D.R.D.; and funding acquisition, L.C. and D.R.D. All authors have read and agreed to the published version of the manuscript.

Funding

The authors are grateful to the funding of the Industrial Research Fund of Ghent University under grant agreement F2017/IOF-StarTT/090.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Huang, S.; Liu, P.; Mokasdar, A.; Hou, L. Additive Manufacturing and Its Societal Impact: A Literature Review. Int. J. Adv. Manuf. Technol. 2013, 67, 1191–1203. [Google Scholar] [CrossRef]
  2. Ben-Ner, A.; Siemsen, E. Decentralization and Localization of Production: The Organizational and Economic Consequences of Additive Manufacturing (3D Printing). Calif. Manag. Rev. 2017, 59, 5–23. [Google Scholar] [CrossRef]
  3. Calignano, F.; Manfredi, D.; Ambrosio, E.P.; Biamino, S.; Lombardi, M.; Atzeni, E.; Salmi, A.; Minetola, P.; Iuliano, L.; Fino, P. Overview on Additive Manufacturing Technologies. Proc. IEEE 2017, 105, 593–612. [Google Scholar] [CrossRef]
  4. Peng, T.; Kellens, K.; Tang, R.; Chen, C.; Chen, G. Sustainability of additive manufacturing: An overview on its energy demand and environmental impact. Addit. Manuf. 2018, 21, 694–704. [Google Scholar] [CrossRef]
  5. Majeed, A.; Zhang, Y.; Ren, S.; Lv, J.; Peng, T.; Waqar, S.; Yin, E. A Big Data-Driven Framework for Sustainable and Smart Additive Manufacturing. Robot. Comput.-Integr. Manuf. 2021, 67, 102026. [Google Scholar] [CrossRef]
  6. Dobrzański, L.A.; Dobrzański, L.B. Dentistry 4.0 Concept in the Design and manufacturing of Prosthetic Dental Restorations. Processes 2020, 8, 525. [Google Scholar] [CrossRef]
  7. González-Henríquez, C.M.; Sarabia-Vallejos, M.A.; Rodriguez-Hernandez, J. Polymers for additive manufacturing and 4D-printing: Materials, methodologies, and biomedical applications. Prog. Polym. Sci. 2019, 94, 57–116. [Google Scholar] [CrossRef]
  8. Franco, D.; Ganga, G.M.D.; de Santa-Eulalia, L.A.; Filho, M.G. Consolidated and inconclusive effects of additive manufacturing adoption: A systematic literature review. Comput. Ind. Eng. 2020, 148, 106713. [Google Scholar] [CrossRef]
  9. Attaran, M. The rise of 3-D printing: The advantages of additive manufacturing over traditional manufacturing. Bus. Horiz. 2017, 60, 677–688. [Google Scholar] [CrossRef]
  10. Khosravani, M.; Reinicke, T. On the Environmental Impacts Of 3D Printing Technology. Appl. Mater. Today 2020, 20, 100689. [Google Scholar] [CrossRef]
  11. 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, 20, 192–204. [Google Scholar] [CrossRef]
  12. Sisi, W.; Lore, C.; Cardon, L.; D’hooge, D.R. Can the melt flow index be used to predict the success of fused deposition modelling of commercial poly(lactic acid) filaments into 3D printed materials. Plast. Rubber Compos. 2018, 17, 9–16. [Google Scholar] [CrossRef]
  13. Boğa, C. Investigation of mechanical and fracture behavior of pure and carbon fiber reinforced ABS samples processed by fused filament fabrication process. Rapid Prototyp. J. 2021, 27, 1220–1229. [Google Scholar] [CrossRef]
  14. Vyavahare, S.; Kumar, S. Numerical and experimental investigation of FDM fabricated re-entrant auxetic structures of ABS and PLA materials under compressive loading. Rapid Prototyp. J. 2021, 27, 223–244. [Google Scholar] [CrossRef]
  15. Parpala, R.C.; Popescu, D.; Pupaza, C. Infill parameters influence over the natural frequencies of ABS specimens obtained by extrusion-based 3D printing. Rapid Prototyp. J. 2021, 27, 1273–1285. [Google Scholar] [CrossRef]
  16. Nathaphan, S.; Trutassanawin, W. Effects of process parameters on compressive property of FDM with ABS. Rapid Prototyp. J. 2021, 27, 905–917. [Google Scholar] [CrossRef]
  17. Duty, C.E.; Kunc, V.; Compton, B.; Post, B.; Erdman, D.; Smith, R.; Lind, R.; Lloyd, P.; Love, L. Structure and mechanical behavior of big area additive manufacturing (BAAM) materials. Rapid Prototyp. J. 2017, 23, 181–189. [Google Scholar] [CrossRef]
  18. Zhou, Z.; Salaoru, I.; Morris, P.; Gibbons, G.J. Additive manufacturing of heat-sensitive polymer melt using a pellet-fed material extrusion. Addit. Manuf. 2018, 24, 552–559. [Google Scholar] [CrossRef] [Green Version]
  19. Go, J.; Schiffres, S.N.; Stevens, A.G.; Hart, A.J. Rate limits of additive manufacturing by fused filament fabrication and guidelines for high-throughput system design. Addit. Manuf. 2017, 16, 1–11. [Google Scholar] [CrossRef]
  20. Liu, X.; Chi, B.; Jiao, Z.; Tan, J.; Liu, F.; Yang, W. A large-scale double-stage-screw 3D printer for fused deposition of plastic pellets. J. Appl. Polym. Sci. 2017, 134, 45147. [Google Scholar] [CrossRef]
  21. Liu, K.; Wang, M.; Li, H.; Zhao, D.; Jin, Y. Effect of characteristic scale on the extrudate swelling behavior of polypropylene melt in a micro-extrusion process. Polym. Eng. Sci. 2021, 61, 1864. [Google Scholar] [CrossRef]
  22. Wilczyński, K.; Buziak, K.; Lewandowski, A.; Nastaj, A.; Wilczyński, K.J. Rheological Basics for Modeling of Extrusion Process of Wood Polymer Composites. Polymers 2021, 13, 622. [Google Scholar] [CrossRef]
  23. Rauwendaal, C.; Gonzalez-Nunez, R.; Rodrigue, D. Polymer Processing: Extrusion. In Encyclopedia of Polymer Science and Technology; John Wiley & Sons: New York City, NY, USA, 2017. [Google Scholar]
  24. Tadmor, Z.; Gogos, C.G. Principles of Polymer Processing; John Wiley & Sons: New York City, NY, USA, 1979. [Google Scholar]
  25. Rauwendaal, C. Polymer Extrusion, 4th ed.; Hanser: Munich, Germany, 2001. [Google Scholar]
  26. Giles Harold, F.; Wagner John, R.; Mount Eldridge, M. Extrusion: The Definitive Processing Guide and Handbook. William Andrew: Norwich, NY, USA, 2004. [Google Scholar]
  27. Dal Fabbro, P.; La Gala, A.; Van De Steene, W.; D’hooge, D.R.; Lucchetta, G.; Cardon, L.; Fiorio, R. Influence of machine type and consecutive closed-loop recycling on macroscopic properties for fused filament fabrication of acrylonitrile-butadiene-styrene parts. Rapid Prototyp. J. 2020, 27. [Google Scholar] [CrossRef]
  28. Maddock, B. A visual analysis of flow and mixing in extruder screws. SPE ANTEC Tech. Pap. 1959, 15, 383–389. [Google Scholar]
  29. Tadmor, Z. Fundamentals of Plasticating Extrusion. I. A Theoretical Model for Melting. Polym. Eng. Sci. 1966, 6, 185–190. [Google Scholar] [CrossRef]
  30. Tadmor, Z.; Duvdevani, I.; Klein, I. Melting in plasticating extuders theory and experiments. Polym. Eng. Sci. 1967, 7, 198–217. [Google Scholar] [CrossRef]
  31. Tadmor, Z.; Klein, I. Engineering Principles of Plasticating Extrusion; Van Nost. Reinhold: New York City, NY, USA, 1970. [Google Scholar]
  32. Donovan, R.C. A theoretical melting model for plasticating extruders. Polym. Eng. Sci. 1971, 11, 247–257. [Google Scholar] [CrossRef]
  33. Chung, C.I. Plasticating single-screw extrusion theory. Polym. Eng. Sci. 1971, 11, 93–98. [Google Scholar] [CrossRef]
  34. Edmondson, I.; Fenner, R. Melting of thermoplastics in single screw extruders. Polymer 1975, 16, 49–56. [Google Scholar] [CrossRef]
  35. Mount, E.M.; Watson, J.G.; Chung, C.I. Analytical melting model for extrusion: Melting rate of fully compacted solid polymers. Polym. Eng. Sci. 1982, 22, 729–737. [Google Scholar] [CrossRef]
  36. Syrjala, S. A New Approach for the Simulation of Melting in Extruders. Int. Commun. Heat Mass Transf. 2000, 27, 623–634. [Google Scholar] [CrossRef]
  37. Altınkaynak, A.; Gupta, M.; Spalding, M.A.; Crabtree, S.L. Melting in a Single Screw Extruder: Experiments and 3D Finite Element Simulations. Int. Polym. Process. 2011, 26, 182–196. [Google Scholar] [CrossRef]
  38. La Gala, A.; Fiorio, R.; Erkoç, M.; Cardon, L.; D’hooge, D.R. Theoretical Evaluation of the Melting Efficiency for the Single-Screw Micro-Extrusion Process: The Case of 3D Printing of ABS. Processes 2020, 8, 1522. [Google Scholar] [CrossRef]
  39. Chung, C. Extrusion of Polymers; Hanser Publications: Cincinnati, OH, USA, 2020; Chapter 2. [Google Scholar]
  40. Pearson, J.R.A. Heat transfer in flowing polymers Progress in Heat and Mass Transfer. In Selected Papers of the 1970 International Seminar; Elsevier: Amsterdam, The Netherlands, 1972; pp. 73–87. [Google Scholar]
  41. Coelho, P.; Pinho, F. A generalized brinkman number for non-Newtonian duct flows. J. Non-Newton. Fluid Mech. 2000, 156, 202–206. [Google Scholar] [CrossRef]
  42. Winter, H. Viscous Dissipation in Shear Flows of Molten Polymers. Adv. Heat Transf. 1977, 13, 205–267. [Google Scholar] [CrossRef]
  43. Tso, C.; Mahulikar, S. The role of the Brinkman number in analyzing flow transitions in microchannels. Int. J. Heat Mass Transf. 1999, 42, 1813–1833. [Google Scholar] [CrossRef]
  44. Kelly, A.; Brown, E.; Coates, P. The effect of screw geometry on melt temperature profile in single screw extrusion. Polym. Eng. Sci. 2006, 46, 1706–1714. [Google Scholar] [CrossRef]
  45. Vera-Sorroche, J.; Kelly, A.; Brown, E.; Coates, P.; Karnachi, N.; Harkin-Jones, E.; Li, K.; Deng, J. Thermal optimisation of polymer extrusion using in-process monitoring techniques. Appl. Therm. Eng. 2013, 53, 405–413. [Google Scholar] [CrossRef]
  46. Abeykoon, C.; Kelly, A.; Vera-Sorroche, J.; Brown, E.; Coates, P.; Deng, J.; Li, K.; Harkin-Jones, E.; Price, M. Process efficiency in polymer extrusion: Correlation between the energy demand and melt thermal stability. Appl. Energy 2014, 135, 560–571. [Google Scholar] [CrossRef] [Green Version]
  47. Nastaj, A.; Wilczyński, K. Optimization and Scale-Up for Polymer Extrusion. Polymers 2021, 13, 1547. [Google Scholar] [CrossRef]
  48. Behzadnasab, M.; Yousefi, A.A.; Ebrahimibagha, D.; Nasiri, F. Effects of processing conditions on mechanical properties of PLA printed parts. Rapid Prototyp. J. 2019, 26, 381–389. [Google Scholar] [CrossRef]
  49. Maloch, J.; Hnátková, E.; Žaludek, M.; Krátký, P. Effect of Processing Parameters on Mechanical Properties of 3D Printed Samples. Mater. Sci. Forum 2018, 919, 230–235. [Google Scholar] [CrossRef]
  50. Acrylonitrile Butadiene Styrene (ABS) Typical Properties Generic ABS. Available online: https://plastics.ulprospector.com/generics/1/c/t/acrylonitrile-butadiene-styrene-abs-properties-processing (accessed on 15 June 2021).
  51. Overview of Materials for Acrylonitrile Butadiene Styrene (ABS), Extruded. Available online: http://www.matweb.com/search/DataSheet.aspx?MatGUID=3a8afcddac864d4b8f58d40570d2e5aa (accessed on 15 June 2021).
  52. Polylactic Acid (PLA) Typical Properties | UL Prospector. Available online: https://plastics.ulprospector.com/generics/34/c/t/polylactic-acid-pla-properties-processing (accessed on 15 June 2021).
  53. Polymer Data Handbook, 2nd ed.; Oxford University Press: Oxford, UK, 2009; pp. 627–635. [CrossRef]
  54. Farah, S.; Anderson, D.G.; Langer, R. Physical and mechanical properties of PLA, and their functions in widespread applications—A comprehensive review. Adv. Drug Deliv. Rev. 2016, 107, 367–392. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  55. SEBS. Available online: https://www.demgy.com/en/datacomp/259 (accessed on 15 June 2021).
  56. SEBS. SEBS-Thermoplastic Polymers, Elastomers and Additives. Available online: https://mexpolimeros.com (accessed on 15 June 2021).
  57. Thermoplastic Elastomer (TPE) Typical Properties Generic SEBS. Available online: https://plastics.ulprospector.com/generics/53/c/t/thermoplastic-elastomer-tpe-properties-processing/sp/18 (accessed on 15 June 2021).
  58. Dreval, V.; Vasil’ev, G.; Borisenkova, E.; Semakov, A.; Kulichikhin, V. Influence of Molecular Weight of SAN on Rheological and Mechanical Properties of ABS-plastics. Appl. Rheol. 2007, 17, 44745-1–44745-11. [Google Scholar] [CrossRef]
  59. Kubota, H. Flow properties of ABS (acrylonitrile–butadiene–styrene) terpolymer. J. Appl. Polym. Sci. 1975, 19, 2299–2308. [Google Scholar] [CrossRef]
  60. Veenstra, H.; van Lent, B.; van Dam, J.; Posthuma de Boer, A. Co-continuous morphologies in polymer blends with SEBS block copolymers. Polymer 1999, 40, 6661–6672. [Google Scholar] [CrossRef]
  61. Wieme, T.; Duan, L.; Mys, N.; Cardon, L.; D’hooge, D. Effect of Matrix and Graphite Filler on Thermal Conductivity of Industrially Feasible Injection Molded Thermoplastic Composites. Polymers 2019, 11, 87. [Google Scholar] [CrossRef] [Green Version]
  62. Duan, L. Extrusion-Based Manufacturing and Structural Engineering for Highly Sensitive Piezoresistive Strain Sensors: A Fundamental Study. Ph.D. Thesis, Ghent University, Faculty of Engineering and Architecture, Ghent, Belgium, 2019. [Google Scholar]
  63. Mount, E.M.; Chung, C.I. Melting behavior of solid polymers on a metal surface at processing conditions. Polym. Eng. Sci. 1978, 18, 711–720. [Google Scholar] [CrossRef]
  64. Tang, D.; Marchesini, F.H.; Cardon, L.; D’hooge, D. State of the-art for extrudate swell of molten polymers: From fundamental understanding at molecular scale toward optimal die design at final product scale. Macromol. Mater. Eng. 2020, 305, 2000340. [Google Scholar] [CrossRef]
Figure 1. Viscosity as a function of the frequency for (a) ABS, (b) PLA, (c) and SEBS (specifications in Table 1), taking the same temperatures as the (target/average) EAM processing parameters. Fitting parameters via the Cross model listed in Table 3.
Figure 1. Viscosity as a function of the frequency for (a) ABS, (b) PLA, (c) and SEBS (specifications in Table 1), taking the same temperatures as the (target/average) EAM processing parameters. Fitting parameters via the Cross model listed in Table 3.
Materials 14 05566 g001
Figure 2. (a) Micro-extruder cooling setup for screw-freezing experiments and (b) micro-extruder barrel cooling chart after stopping the extrusion with active water-cooling (T1: thermocouple 1), (c1c3) embedded samples for ABS at screw frequency of 2, 5, and 8 rpm, respectively, from left to right, (d) position of the thermocouples, units in mm.
Figure 2. (a) Micro-extruder cooling setup for screw-freezing experiments and (b) micro-extruder barrel cooling chart after stopping the extrusion with active water-cooling (T1: thermocouple 1), (c1c3) embedded samples for ABS at screw frequency of 2, 5, and 8 rpm, respectively, from left to right, (d) position of the thermocouples, units in mm.
Materials 14 05566 g002
Figure 3. Melt removal by the drag mechanism: (a) reported by Maddock and Tadmor [28,29,30,31] for large scale SSEs and initially observed in the present work; (b) outer recirculation flow observed in several cases in the conical screw micro-extruder. In the end, the overall measuring/modeling result remains the same for subplot a and b through quantification of the area of the molten part vs. the total available space.
Figure 3. Melt removal by the drag mechanism: (a) reported by Maddock and Tadmor [28,29,30,31] for large scale SSEs and initially observed in the present work; (b) outer recirculation flow observed in several cases in the conical screw micro-extruder. In the end, the overall measuring/modeling result remains the same for subplot a and b through quantification of the area of the molten part vs. the total available space.
Materials 14 05566 g003
Figure 4. (a) Predicted and experimental individual melting profiles for ABS for three runs at 5 rpm; (b) predicted profile and average experimental points for ABS processed at 5 rpm (Table 3).
Figure 4. (a) Predicted and experimental individual melting profiles for ABS for three runs at 5 rpm; (b) predicted profile and average experimental points for ABS processed at 5 rpm (Table 3).
Materials 14 05566 g004
Figure 5. (a) Predicted and experimental melting profiles for embedded ABS samples at: (a) 2 rpm; (b) 5 rpm; (c) 8 rpm (Table 3).
Figure 5. (a) Predicted and experimental melting profiles for embedded ABS samples at: (a) 2 rpm; (b) 5 rpm; (c) 8 rpm (Table 3).
Materials 14 05566 g005
Figure 6. (a) Predicted and experimental melting profiles for PLA; (b) predicted and experimental melting profiles for SEBS; (c) predicted melting profiles for ABS, PLA, SEBS for different screw frequencies (Table 3).
Figure 6. (a) Predicted and experimental melting profiles for PLA; (b) predicted and experimental melting profiles for SEBS; (c) predicted melting profiles for ABS, PLA, SEBS for different screw frequencies (Table 3).
Materials 14 05566 g006
Figure 7. Predicted Brinkman number (a) and predicted melting profiles (b) for ABS (blue), SEBS (red) and PLA (black) at 5 rpm. The lines in (b) are taken from Figure 5 and Figure 6 (Table 3).
Figure 7. Predicted Brinkman number (a) and predicted melting profiles (b) for ABS (blue), SEBS (red) and PLA (black) at 5 rpm. The lines in (b) are taken from Figure 5 and Figure 6 (Table 3).
Materials 14 05566 g007
Figure 8. Left: Model predictions for melting profiles for ABS at screw frequencies of 0.5, 5, 20 rpm. Right: SEM images of printed filaments; enlarged and extra SEM images are included in Supplementary Materials Figures S12–S17.
Figure 8. Left: Model predictions for melting profiles for ABS at screw frequencies of 0.5, 5, 20 rpm. Right: SEM images of printed filaments; enlarged and extra SEM images are included in Supplementary Materials Figures S12–S17.
Materials 14 05566 g008
Figure 9. (a) Tensile strain; (b) modulus; (c) final diameter for ABS filament extruded at 0.5, 5, 20 rpm. The other operating conditions are given in Table 3.
Figure 9. (a) Tensile strain; (b) modulus; (c) final diameter for ABS filament extruded at 0.5, 5, 20 rpm. The other operating conditions are given in Table 3.
Materials 14 05566 g009
Table 1. Material characteristics for ABS (column 2), PLA (column 3), and SEBS (column 4); “melting” temperature (Tm) based on DSC measurements (see Section 3.2); glass transition temperature (Tg). Data from Refs. [50,51,52,53,54,55,56,57].
Table 1. Material characteristics for ABS (column 2), PLA (column 3), and SEBS (column 4); “melting” temperature (Tm) based on DSC measurements (see Section 3.2); glass transition temperature (Tg). Data from Refs. [50,51,52,53,54,55,56,57].
ABSPLASEBSUnit
Density (at 300 K)1.01–1.081.22–1.300.87[g/cm3]
Tg381-382326–337183 (Butadiene)–373 (Styrene)[K]
Tm443–593418–459455–553[K]
Tensile strength, break22–59.3591.80–12.00[MPa]
Elongation at break2.5–40<10%310–860[%]
Modulus of Elasticity1.5–2.63.50.6–4.68[GPa]
MorphologyAmorphousSemicrystallineThermoplastic elastomer
Table 2. Machine geometry and dimensions for extrusion-based additive manufacturing in the present work.
Table 2. Machine geometry and dimensions for extrusion-based additive manufacturing in the present work.
Characteristic Unit
Axial lengthLa662[mm]
- FeedingLf280[mm]
- CompressionLc274[mm]
- MeteringLm108[mm]
Screw diameter-initialDi28[mm]
Screw diameter-finalDf16[mm]
Helical lengthLh690[mm]
Channel depth-initialHi6[mm]
Channel depth-finalHf2[mm]
Channel widthW16[mm]
Pitch angleΘ0.308[rad]
Table 3. Summary of materials properties and operating conditions.
Table 3. Summary of materials properties and operating conditions.
Material Properties ABSPLASEBSUnitSection
Cross model parameterη0115,126175510,700[Pa s]Section 3.1
Cross model parameterτ*49,521265,95479,318[Pa]Section 3.1
Cross model temperatureTcr483453463[K]Section 3.1
Pseudo-plasticity indexN0.240.18560.39[/]Section 3.1
Specific heat capacity solidcs127313151993[J kg−1 K−1]Section 3.2
Specific heat capacity meltcm227724252607[J kg−1 K−1]Section 3.2
Specific heat capacity (avg.)c177518702300[J kg−1 K−1]Section 3.2
Thermal conductivity (avg.)K0.1910.2270.331[W m−1K−1]
Heat of fusionλ047,0000[J kg−1]Section 3.2
“Melt” temperatureTm423442423[K]Section 3.2
Glass transition temperatureTg373331238[K]Section 3.2
Density (solid)ρs10501240910[kg m−3]Section 3.2
Density (melt)ρm9791119796[kg m−3]Section 3.3
Operating Conditions ABSPLASEBSUnitSection
Initial TemperatureT0300300300[K]/
Barrel TemperatureTb483453463[K]/
Screw frequencies (for model validation)N2,5,855[rpm]/
Volumetric inlet flow at 2 rpmQ0141.3--[mm3 s−1]/
Volumetric inlet flow at 5 rpmQ0104.398.1109.4[mm3 s−1]/
Volumetric inlet flow at 8 rpmQ047.1--[mm3 s−1]/
Table 4. Thermal conductivity measurements for ABS, PLA, and SEBS from Hot disk tests.
Table 4. Thermal conductivity measurements for ABS, PLA, and SEBS from Hot disk tests.
Material 1Temperature [K]Measured k [W m−1 K−1]Error [W m−1 K−1]
ABS2960.181±0.0003
ABS3730.200±0.0010
PLA2960.224±0.0020
PLA3830.229±0.0030
SEBS2960.331±0.001
1 The experiments were conducted according to ASTM standard E1269 at 50% relative humidity.
Table 5. Specific heat capacity for ABS, PLA, and SEBS from DSC at a heating rate of 5 K min−1.
Table 5. Specific heat capacity for ABS, PLA, and SEBS from DSC at a heating rate of 5 K min−1.
Temperature [K]ABSPLASEBSUnit
Specific heat capacitycs300127313151993[J kg−1 K−1]
Specific heat capacitycm460227724252607[J kg−1 K−1]
Table 6. Measurements toward melt density for ABS, PLA, and SEBS.
Table 6. Measurements toward melt density for ABS, PLA, and SEBS.
MaterialMVR
[cm3/10 min]
Error [cm3/10 min]Extruded Weight
[g/10 min]
Temperature [K]External Weight [kg]Melt Density [g/cm3]
ABS8.187±0.4017.910493100.966
ABS8.809±0.1388.727493100.991
PLA6.975±0.1507.8524732.161.126
PLA7.365±0.2868.1864732.161.111
SEBS3.701±0.0272.9304732.160.792
SEBS3.639±0.0322.9084732.160.799
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

La Gala, A.; Fiorio, R.; Ceretti, D.V.A.; Erkoç, M.; Cardon, L.; D’hooge, D.R. A Combined Experimental and Modeling Study for Pellet-Fed Extrusion-Based Additive Manufacturing to Evaluate the Impact of the Melting Efficiency. Materials 2021, 14, 5566. https://doi.org/10.3390/ma14195566

AMA Style

La Gala A, Fiorio R, Ceretti DVA, Erkoç M, Cardon L, D’hooge DR. A Combined Experimental and Modeling Study for Pellet-Fed Extrusion-Based Additive Manufacturing to Evaluate the Impact of the Melting Efficiency. Materials. 2021; 14(19):5566. https://doi.org/10.3390/ma14195566

Chicago/Turabian Style

La Gala, Andrea, Rudinei Fiorio, Daniel V. A. Ceretti, Mustafa Erkoç, Ludwig Cardon, and Dagmar R. D’hooge. 2021. "A Combined Experimental and Modeling Study for Pellet-Fed Extrusion-Based Additive Manufacturing to Evaluate the Impact of the Melting Efficiency" Materials 14, no. 19: 5566. https://doi.org/10.3390/ma14195566

APA Style

La Gala, A., Fiorio, R., Ceretti, D. V. A., Erkoç, M., Cardon, L., & D’hooge, D. R. (2021). A Combined Experimental and Modeling Study for Pellet-Fed Extrusion-Based Additive Manufacturing to Evaluate the Impact of the Melting Efficiency. Materials, 14(19), 5566. https://doi.org/10.3390/ma14195566

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